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THE TOWNSEND DISCHARGE IN A COAXIAL DIODE 
WITH AXIAL MAGNETIC FIELD! 


P. A. REDHEAD 


ABSTRACT 


An approximate theory is developed of the breakdown characteristics of a 
coaxial diode in an axial magnetic field, taking into account the effects of elastic 
collisions. It is assumed that the electron moves in a constant electric field 
between collisions and thus the theory is valid only in the appropriate range 
of magnetic field and voltage. Estimates of transit time and of space-charge 
effects are also made. Measurements in the pressure range 10-% to 10-* mm. Hg 
are in general agreement with the theory. 


1, INTRODUCTION 

A self-sustained Townsend discharge in a coaxial diode with an axial 
magnetic field can be used very effectively as a means for measuring extremely 
low gas pressures. The theory of the striking characteristics of such discharges 
has been studied by Somerville (1952) and by Haefer (1953a and 6, 1954), and 
the operation of vacuum gauges using this type of discharge has been de- 
scribed by Beck and Brisbane (1952) and Haefer (19530, 1954). 

This paper presents a theory of the striking characteristics of coaxial 
cylinders in an axial magnetic field which takes into account the effects of 
elastic collisions of the electrons. Measurements of the striking characteristics 
are also presented, which were made on a number of magnetron diodes designed 
for the measurement of extremely low pressures (below 10-!® mm. Hg). The 
operation of these diodes as pressure gauges is described in a companion paper 
(Hobson and Redhead 1958). 

A schematic diagram of the electrode arrangement used for the measure- 
ments is shown in Fig. 1. The inner cylinder was used as the anode (radius 7,4) 
and the outer cylinder as the cathode or ion-collector (radius r;,). The ends 
of the cathode cylinder were partially closed to prevent migration of electrons 
from the discharge in an axial direction. The auxiliary cathode is a box-like 
structure enclosing the ion-collector and acting as an electrostatic shield; for 

‘arity this is not shown in Fig. 1. Connected to the auxiliary cathode are 
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two short cylindrical shields which project through holes in the end-plates of 
the cathode. These cylindrical shields reduce the electric field at the edges 
of the end-plates and suppress any field-emission from the cathode. The 
structure of this tube suggested that it be christened the ‘‘inverted-magnetron 
gauge’ (abbreviated to I.M.G. hereafter). 


CATHODE 
a 









CATHODE 
END-PLATES 
ANODE 


AUXILIARY 
CATHODE 


Fic. 1. Schematic diagram of diode. 


The theory of the striking characteristics will be developed for infinitely 
long cylinders, i.e. neglecting the effects of the end-plates. This approximation 
is reasonable on the upper branch of the striking characteristic since, as will 
be seen later, the breakdown condition is first satisfied at the mid-point of 
the cathode cylinder where the probability of electron escape is greatest and 
the distortion of the electric field by the end-plates is least. 


2. THEORY OF STRIKING CHARACTERISTICS 

(a) Electron and Ion Trajectories 

Fig. 2 shows a typical potential distribution measured with an analogue 
field plotter for I.M.G. No. 16 (see Table II for dimensions). It will be seen 
that an electron liberated within the ion-collector cylinder is trapped axially 
by the electric field. Provided the magnetic field is in excess of the cutoff 
field, which is the only case considered, the electron is trapped radially by 
the magnetic field. Electrons will follow a hypocycloidal path around the anode, 
while simultaneously oscillating back and forth in the axial (Z) direction. 
Approximate electron trajectories are indicated in Fig. 8 where the crosses 
indicate inelastic collisions. Electrons leaving the cylindrical portion of the 
cathode surface perform one arch of a cycloid before striking the cathode 
surface again; if they suffer an inelastic collision in this first arch they are 
assumed to return to thermal energy and to start their cycloidal motion 
once more from the point of collision, and thus to escape from the cathode. 
Electrons which have escaped from the cathode and which do not suffer 
further collisions will circle the anode in a hypocycloid generated by a circle 
of diameter D rolling on a circle of constant radius. The mean 6-directed 
velocity of the electrons is E/B (EK, electric field; B, magnetic field). Thus 
electrons can only advance towards the anode by suffering ionizing collisions. 

The trajectories of the positive ions formed by inelastic collisions of the 
electrons are relatively unaffected by the magnetic field because of the large 
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in Table I. 


Fic. 2. Measured potential distribution in I.M.G. No. 16 with the dimensions indicated 
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Fic. 3. Approximate electron trajectories. 
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ionic mass. The positive ions will move to the cathode with a transit time 
which is very much shorter than the electron transit time. 

The above description of the trajectories is greatly simplified since the 
effects of space-charge and the distortion of the electric field by the end- 
plates have been ignored. The situation is further complicated by the presence 
of electronic oscillations in the gauge which are of considerable amplitude at 
higher pressures. 

The theory of the striking characteristics will be developed on the assump- 
tion that the electric field is that which would exist in a cylindrical condenser 
of infinite length, and the effects of space-charge and of oscillations will be 
ignored. 


(6b) Derivation of a 

Previous derivations of the first Townsend coefficient (a) in crossed electric 
and magnetic fields (Kugler et al. 1933; Haefer 1953a and 6; Somerville 1952 
have ignored the contribution to a made by those electrons which have been 
elastically scattered. The following analysis will take into account the elasti- 
cally scattered electrons by an approximate method. We will first consider 
the case of crossed, uniform electric and magnetic fields. The electron tra- 
jectories in such a system are cycloidal, as indicated in Fig. 4. The diameter 
of the rolling circle generating the cycloid is 


(1) D = 2E/nB?, 


where 7 is the charge-to-mass ratio of the electron. 


B® 


Y 





Fic. 4. Electron trajectories in crossed, uniform electric and magnetic fields. 


The arc length of the cycloid from a point of height xo to a point of height 
x is given by 
S, = 2D[1—(1—x/D)!”]—So, 
= 2D[((1—x0/D)!?”—(1—x/D)'"], 
and similarly 


Sz = 2D[(1—x0/D)'?+(1—x/D)'"], 


where So is the path length to a point of height xo. xo is the value of x at 
which electrons have sufficient energy to enter into inelastic collisions; thus 
xo = V,/FE, where V; is the ionization potential of the gas. It will be assumed 
that for x < xo only elastic collisions are possible, while for x > xo inelastic 
collisions can occur. 
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Electrons that are elastically scattered in the region x < Xo will perform 
orbits that do not reintersect the plane x = 0, and will in almost all cases 
eventually pass through the plane x = xo and be capable of entering into 
inelastic collisions, thus contributing to the avalanche process. In previous 
treatments of this problem it was assumed that the elastically scattered 
electrons do not enter into the avalanche process, thus underestimating a. 
We shall assume that all the elastically scattered electrons are capable of 
suffering inelastic collisions, at some later time, and contributing to the 
avalanche process. 

The contribution of the elastically scattered electrons will be taken into 
account in an approximate way by assuming that the number of electrons 
available at the plane x = x» is the same as that which would exist if no 
elastic collisions occurred in the region x < xo. Thus, referring to Fig. 4, the 
number of electrons at point A is mo, and all the electrons at point B are still 
available at point C. For this assumption to be valid, the average path length 
above xo of the electron elastically scattered below x» must approximate the 
path length above xo of the unscattered electron. Calculations have been 
made of the path length, above x», of the electron which is elastically scattered 
into the negative y-direction at x = xo. This electron has the maximum path- 
length above x» of any electron and this path length can be approximated as 


S = SV (i+0/—*), 
where r= V x0/D, 





and S, is the path length, above xo, of the unscattered electron. 

For xo/D = 0.36, S/S, = 2. Thus for xo/D < 4, the maximum path-length 
of an elastically scattered electron is less than twice that of the unscattered 
electron. Thus the assumption that the average path-length, above Xo, of the 
elastically scattered electron approximates that of the unscattered electron is 
reasonable, provided that xo/D is not large. 

The number of electrons at A being mo, then the number at point 1 is 


no exp(—Si/A), 


Ns, 


and, similarly, at point 2 is 


Il 


Ns, = no exp(—S2/X), 


where \ is the electron mean-free-path for all collisions and is assumed 
independent of electron energy. 

The number of electrons which have suffered inelastic collisions with energies 
greater than e/x (e, electronic charge) in the first cycloidal arch is 


: 2D tev ts 2 x \'/? 
(2) n(x). = Ns,—Ns, = 2No exe] -22(1-#) sinh - I-h F 


The arc length from A to B is 4D(1—x0/D)'” and thus the number of electrons 
at B is given by 
(33) Np = no exp[—(4D/d)(1—x0/D)!'”]. 
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Combining equations (2) and (3) we can show that the total number of electrons 
with energies greater than ex which have suffered inelastic collisions is given 
by 
n(x) l< sinh[(2D/d)(1—x/D)"”” 
(4) “<2 = — Yon(x); = a (2 soe eae : 
No No ‘m1 sinh[(2D/d)(1—x0/D)’"] 


When 2D/r\ « 1, then 


(5) n(x) _ | 1—x/D i 
. mo L1—x0/D : 


which is the ratio of the path length above x to the path length above x». 
By differentiation of equation (4) we find that the number of electrons 
that have suffered inelastic collisions in the energy range eEx to eH(x+6x) is 


dn(x) _ cosh[(2D/d)(1—x/D)"”"](1—x/D)"" 


_- = \ sinh[(2D/) (1—x0/D)'"} 


Then the most probable value of x for an inelastic collision to occur is 


aD 
x = J (x /no)dn(x) 
2 *? » cosh{(2D/d) (1—x/D)*"] 


7 d sinh (2D/X) (1—Xo, D)'*}. es (1—x/D)” ax. 


Integration then leads to 


1/2 1/2 
xv a Xo __ Xo f ; 2D = | = 
(6) a ¢ ss) coth r (: D op ° 


When 2D/\X is small, then 
(7) ¥ = 4(x9+2D); 


Haefer (1953a) finds @ = 2D, which is meaningless when x9 > 3D. 

The first Townsend coefficient, a, defined as the number of ion pairs formed 
per unit distance of electron advance in the direction of the electric field, is 
given by 
(8) ee I n(x») fa 4 . 

gE No ge 
where g is a correction factor, greater than unity, to account for the fact 
that when the electron energy exceeds eV; the probability of a collision causing 
ionization is not unity. 

Combining equations (1), (7), and (8), we find that for 2D/A < 1, 


3E 
(9) QS Fa fee et 
g(4E"/nB°+ Vi) 
where we see that a@ is independent of the mean free path. This independency 


arises in the uniform field case because we have assumed the discharge to be 
of infinite extent in the lateral direction. In a coaxial system the electron 
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orbits are closed and the electron, being trapped in the system, will circle the 
anode indefinitely at a fixed mean radius until collision occurs; thus we would 
again expect @ to be independent of X. 

This method of deriving @ and a is only applicable to cylindrical geometry 
when the cycloid height D(r) is sufficiently small that the electric field can 
be considered constant over the cycloid. When D(r) becomes an appreciable 
fraction of the cathode radius, the uniform field assumption is no longer valid, 


and the above analysis is inapplicable. 


(c) The Breakdown Equation for Inhomogeneous Fields 

For a discharge in a cylindrical condenser with zero magnetic field the 
electron multiplication is not given by exp fa dr (Morton 1946) since the 
electron energy distribution is not the same as in a constant electric field 
where the distribution has had sufficient time to come to equilibrium with 
the field. However, when an axial magnetic field is applied, the electrons 
advance in a radial direction very slowly and the electron energy distribution 
can reach equilibrium with the field. Thus the Townsend breakdown relation 
for a coaxial diode in a magnetic field may be written 


(10) 1 exp ( { ata | = |, 


where ro is the maximum radius at which ionization can occur, and T is the 
effective third Townsend coefficient for all secondary processes. Combining 
equations (9) and (10), we obtain 


(11) log( 1+ 4 = | adr= 2 Ext lol 14-2! ey | 
I J ;, 2 gV; | Exr; 


when the upper limit of the integral is taken as zero (i.e., 7o/r%a >> 1) and E;, is 
the field at the cathode surface. 
The characteristic may be divided into two regions: 
(a) where the maximum energy of the electron in a cycloidal arch originating 
at the cathode surface is greater than the ionizing energy, 
Vi (max) = 2, Dy = 26 e/nb > Vi: 


and (b) where the maximum cycloidal energy at the cathode is less than the 
ionizing energy, 
V.(max) < V;. 


In region (a) the secondary electrons emitted from the cylindrical surface of 
the cathode can cause ionization in their first cycloidal arch and enter the 
discharge directly. In region (6) electrons emitted from the curved surface of 
the cathode cannot cause ionization in their first cycloidal arch and thus 
most of these electrons will be recaptured at the cathode. Some electrons can 
escape by making elastic collisions but this number is extremely small at 
low pressures. In this region the effective gamma is essentially controlled by 
the end-plates of the cathode cylinder and by processes within the discharge 


volume. 
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In region (a), 79 = r,, and thus equation (11) becomes 


1 3 Er, | nV; ay 
? y = a" —— a ae 
(12) joe( 1+) 2 eV, log} 1+ ie , 


For (B/E;,)? small, 
3y(Br,)* 
(13) ae. i : 
g log(1+1/T) 
Equation (13) represents the upper parabolic branch of the striking character- 
istic shown in Fig. 5, and agrees with the limiting form found by Haefer (19538, 
Eq. 16). 


Vs 


Ve B? 
— 
REGION 4) 
— VY, (max)= Vj 
- 
= 
oa 
a= 
a 
- REGION (b) 





Fic. 5. Theoretical striking characteristic. 


In region (6) ro = 71, where 7; is the maximum radius at which an electron 
can acquire ionizing energy in its cycloidal path; thus 7 is the radius at which 
V(max) = V; or 


= furs)! 
(14) ry -2( B ; 


Substituting this value of r; for ro in equation (11) leads to 
(15) Eyr, = (5/3)gV;logA+1/P). 


Thus in region (0) the striking voltage will be independent of magnetic field 
provided that I remains constant. 

The minimum anode striking voltage, on the lower branch, derived from 
equation (15) is 


(16a) V,(min) = (5/3)gV; log(r;/ra) log A+1/TP). 


The minimum magnetic field is found by substituting equation (14), with 
r, = rz, in (15), which gives 


(16d) Br,(min) = (5/3)g(2V;/n)!” log(1+1/P). 
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(d) The Gamma Processes 

In region (a) of the striking characteristic the electrons emitted from the 
curved surface of the cathode have a finite probability of escaping from the 
cathode and entering the discharge. An electron emitted from the curved 
surface which performs one complete cycloidal arch without suffering a 
collision is recaptured by the cathode and fails to enter the discharge. Elec- 
trons which suffer inelastic or elastic collisions follow orbits which do not 
intersect the cathode surface a second time. Thus the effective gamma of the 
curved surface (I';) is the intrinsic gamma (y:) modified by the probability 
of escape of the electron from the curved surface. If mo electrons are emitted 
from the curved surface, only mo exp(—4D,/A) electrons will return to the 
cathode. Thus the effective gamma is given by 


T; = vi{l—exp(—4D, /r)], 
which, for 4D,/\ « 1, reduces to 
(1) Tr; — 47D, ‘Xr. 


When the effective gamma of the curved surface is large compared to that 
of the end-plates, the upper branch of the characteristic is controlled by T), 
and by substituting equation (17) in (166) we find that 


(18) Br,(min) = a, log(d:/p), 


where p is the pressure, and a; and }; are constants. 

When the cathode escape probability (4D,/A) is large, and thus the I; 
process is important on the upper branch of the characteristic, then D,/r;, 
will also be large at low pressures (large \). Under these conditions the assump- 
tion that the electric field is constant over the cycloid height is no longer 
valid and the foregoing analysis cannot be applied to the upper branch of 
the characteristic. 

In region (6) of the striking characteristic, the cathode escape probability 
is essentially zero. Secondary electrons emitted from the end-plates can enter 
the discharge with their axial velocity component unaffected by the magnetic 
field. These electrons can contribute to the formation of a discharge, provided 
that they form ion pairs at a radius where the maximum cycloid energy 
exceeds the ionization potential (V(max) > V,), so that the slow electron from 
the ion pair is capable of causing further ionization. The orbits of the electrons 
emitted from the end-plates are extremely complex, and a complete formu- 
lation of the effective end-plate gamma, Ts, is not possible. However, an 
approximate calculation has been made of the longest path length of an 
electron emitted from an end-plate before it reaches the opposite end-plate. It 
was assumed that the electron remains at a constant mean radius during its 
transit of the gauge and that the cycloid heights were small in all cases. A 
maximum value of the 6-directed velocity was found as a function of the 
radius at which the electron was emitted by using the theoretical value of 
radial field. This velocity was compared with the minimum value of the 
mean axial velocity at the same radius obtained from the measured potential 
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distributions. For the standard operating conditions the path length of an 
electron which travels an axial distance / is given by 

Si < 1.221. 

Thus it is reasonable to assume that [, will vary in the same way as T;, with 
the axial length of the gauge / replacing 4D, as the characteristic length; thus 
lr, = yol/X, for 1//X <1. 

Substituting this value of IT, in equation (16a) gives 
(19) V,(min) = az log(bs/p), 
where a2 and bs are constants. 

When T[, is very small the effect of the end-plates on the discharge is 
negligible and the lower branch of the striking characteristic will rise and 
tend to follow the line V,(max) = V, (see Fig. 5). This rising characteristic 
of the lower branch has been observed by Penning (1936) in long discharge 
tubes at low pressures, where the influence of the end-plates on the discharge 
will be minimized. 

3. TRANSIT TIMES AND SPACE CHARGE 

The electron path length for one complete revolution of the electron around 

the anode is given by 


arc length of the cycloid 








L = 2ar.—— 


circumference of the rolling circle 
= 2nrr.4D/rD = 8r. 


The electron will then circle the anode \/8r times between collisions, and the 
time between collisions is given by 
\/8r.2a4r ; 
te = = Am /4v;, 


” 


Ct 
where v, is the average tangential velocity of the electrons given by v, = //B. 
The average radial distance moved between collisions is #, which, for 
2D/ « 1, can be approximated as (see equation (7)) 


I= A(xo+2D). 


The radial velocity of the electrons is then given by 


a AV; 4 Ets)’ 
eh es i mr BA ae SH Br . 


The maximum electron transit time is given by 


” {" dr Ee Be mo 2h if; tan ( Bry a) | 
. ai Pro Ur a 4V; ; B(V in)” Eegre 2 


if the upper limit of the integral is taken as zero; i.e€., 7a K rp. 
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In region (a), 7o = rx, and when B/E, < 1, we find that 


3 
‘107° ‘ x= rem | B 
(19a) t_(a) n 16 E, . 


In region (6), ro = 1, and substituting r; from equation (14), we find that 


, wrE, 1. 

1¢ ze => — 7s 

(19d) t_(b) 5On) V2 

For air, V; = 15 volts, \ = 5X10~-5/p (p in mm. Hg), and taking r;/7, = 30, 
we find that 


7t(b) = 8.04X10-" V4/p. 


The trajectories of the positive ions are only slightly affected by the magnetic 
field because of the large mass of the ions. The ionic transit-time will be 
estimated as the transit-time of an ion in a cylindrical diode with zero magnetic 
field and no space-charge. The radius of the diode will be arbitrarily taken 
as 2r,. Using the data of Ferris (1936) we find the transit-time, as defined 
above, for ions of mass 29 to be 
1.2710 *r, 


ie ae 
V/ 


for f2/%, = 30: 
Values of the electron and ion transit-times are shown in Table I for the 
following conditions: 


Va = ORv., t/t, = 30, 4 = Kicm., M—=2o Ve=— 15 volts, 


and A = 5X10°*/p cm. 


TABLE | 


2 
(mm. Hg) 7_(b) (sec.) t-(b)/r, Space-charge 





10-2 4.8X10-7 0.2 + 
10-4 18X10~5 20 — 
10-6 1.81073 2X 108 “ 
10-8 t.8X10"! 2x 105 _ 
10710 18 2X10? = 

2x 109 


10-* 


4.8X 108 

The tabulated values are very approximate since the distortion of the 
assumed electric fields by the end-plates will seriously affect the estimated 
transit-times. Also, at the higher pressures A/r, will approach unity and affect 
the estimate of r_/7,. However, these values may be used as a rough guide 
to indicate the order of magnitude of these parameters. 

For a pressure of about 2X 10-* mm. Hg r_/r, is unity; thus for pressures 
above this value the space-charge will be positive and for lower pressures the 
space-charge will be negative. Thus at pressures above 2X10~* mm. Hg, the 
discharge current will decrease with applied voltage, and at lower pressures 
it will increase with voltage. 
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4. MEASURED STRIKING CHARACTERISTICS 
Figs. 6 through 8 show the measured striking characteristics of tubes 19, 
16, and 4, respectively, with dry air as the test gas. The salient dimensions 


I.M.G. No.19_ 
=0°5 


"a 


Leoem 9-0 


4 =30°0 











0 200 400 600 800 1000 
B [Gauss] 


Fic. 6. Measured striking characteristics for I.M.G. No. 19. 


of these tubes are given in Table II. The upper dashed parabola in these 
three figures is the locus of the condition D;,/r, = 1/10. This value of D,/ 7; 
has been arbitrarily chosen as the maximum value for which the constant 
field assumption used in the foregoing analysis can be considered applicable. 
Thus we cannot expect the characteristics in the region above this line to be 
explainable on this simple theory. For I.M.G. Nos. 16 and 4 all the upper 
branch characteristics, in the pressure range used in the measurements, lie 
above the D,/r, = 1/10 line. These characteristics are almost lines of constant 
magnetic field with a slight negative slope at high pressures. The character- 
istics of I.M.G. No. 19 which lie below this line show the approximate parabolic 
form that the constant field theory predicts. 


TABLE II 
(All dimensions in millimeters) 


I.M.G. No. 





19 16 4 3 
Ya 0.5 0.5 0.125 0.5 
rk 9 15 15 15 
l 30 13.5 15 15 
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I.M.G. No.I6 


| V, (mox)= V; 
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Fic. 7. Measured striking characteristics for 1.M.G. No. 16. 
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Fic. 8. Measured striking characteristics for 1.M.G. No. 4. 
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The lower dashed line in the three figures is the locus of the condition that 
the maximum energy of an electron in a cycloid originating at the curved 
surface of the cathode is equal to the ionizing energy (V,;(max) = V;,) for 
V= 15 volts: 

The measurements for I.M.G. Nos. 4 and 16 were continued out to magnetic 
fields of 3X 10% gauss. In both cases the lower branch from 10° to 3X 108 gauss 
was exactly horizontal. For I.M.G. No. 19, the maximum field obtainable 
with the solenoid used with this tube was 10? gauss. 

Measurement of the striking voltage below 10-* mm. Hg was impossible 
to achieve with any accuracy because of the large and variable time-lag 
between the application of the voltage and the establishment of the dis- 
charge. This time-lag increased as pressure was decreased, and was as long 
as 10 minutes at pressures near 10-” mm. Hg. However, no tubes failed to 
strike even at the lowest pressures obtainable. 

kKig. 9 shows striking voltage as a function of pressure for a magnetic 
field well out on the lower branch of the characteristic (2060 gauss) for two 








P [mm. Hg] 


Fic. 9. Striking voltage vs. pressure for B = 2060 gauss. 


tubes, differing only in their anode diameters (see Table II for dimensions). 
This value of magnetic field was used when operating the tubes at very low 
pressures. A linear relation between V,, and log p was observed, in accordance 
with equation (19); this linear relation was closely obeyed by all the tubes 
at high magnetic fields (ca. 2X 10° gauss). 

I.M.G. No. 8 closely approximated the dimensions of the gauges used 
for pressure measurements below 10-'? mm. Hg. Extrapolation of the curve 
of Fig. 9 to pressures of 10- mm. Hg indicates a striking voltage of 2.2 kv., 
which is in general agreement with rough estimates made at that pressure. 
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From the intercepts of the curves of Fig. 9 on the p-axis and assuming 
that T2 = yol/A, we can estimate the value of y2, the intrinsic gamma of 
the end-plates. For I.M.G. No. 3, y. = 710-4, and for I.M.G. No. 4, 
vo = 1>€10-. 

Fig. 10 shows Br,(min) as a function of pressure for three tubes whose 
dimensions are given in Table II. These data are a reasonable fit to a straight 
line on a Br,(min) vs. log p plot, in agreement with equation (18). For I.M.G. 
No. 16 there is considerable deviation from the straight line for pressures 


above 10-4 mm. Hg. 


600 


200 
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Fic. 10. Minimum magnetic field vs. pressure. 


5. COMMENTS 

The simple theory developed on the assumption that the electric field is 
constant over any one cycloidal arch is in agreement with the measured 
striking characteristics in region (6) and is in qualitative agreement with the 
observed striking characteristics in region (a) when D,/r, is sufficiently small. 
This theory cannot be applied in that part of region (a) where D,/r, is large 
and steeply rising characteristics are observed. An explanation of the charac- 
teristics when D,/r, is large awaits the development of a theory which takes 
account of the more complex electron paths that occur in this region. 

The measured characteristics of this type of discharge are significantly 
affected by electronic oscillations in the discharge. Any theoretical treatment 
of the discharge characteristics that ignores these oscillatory effects can only 
give a first-order approximation to the measured characteristics; in particular, 
such a theory will be most in error at higher pressures. 
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OPERATION OF AN INVERTED-MAGNETRON GAUGE IN 
THE PRESSURE RANGE 10-* TO 10-" MM. Hg! 


J. P. Hopson AND P. A. REDHEAD 


ABSTRACT 


A cold-cathode ionization gauge with axial magnetic field and radial electric 
field is described which is useful in the pressure range 10~* to 10-"% mm. Hg and 
has a sensitivity of about 1 amp./mm. Hg. This gauge has the structure of an 
inverted magnetron with an auxiliary cathode which provides the initial field 
emission and allows the positive-ion current to be measured independently of 
the field-emission current. The operation and calibration of the gauge in a 
dynamic vacuum system is described for the pressure range 107’ to 5X 10-9 mm. 
Hg. Pressures ranging from 1078 to 107% mm. Hg were obtained in sealed-off 
systems. Liquid-helium trapping was used for pressures below 10~'® mm. Hg. It 
is shown that the gauge gives a reliable indication of pressure below 107! mm. 
Hg when the pressure is determined by helium diffusion through the envelope. 


1. INTRODUCTION 


Measurement of extremely low pressures is normally performed with a 
triode ionization gauge. A lower limit to the pressure measurable with such 
a gauge is reached when the ion current becomes so small that it is masked 
by the pressure-independent photocurrent to the ion collector. This photo- 
current is produced by X-rays released from the grid structure by the ionizing 
electron-current. The Bayard—Alpert gauge (Bayard and Alpert 1950) was 
designed to reduce this effect by decreasing the surface area of the ion collector 
to a minimum; the lowest measurable pressure with this type of gauge is 
about 107!° mm. Hg. 

The possibility of using a Penning type of cold-cathode gauge for the 
measurement of pressures below 10-!° mm. Hg was examined, since this type 
of gauge appeared to have two advantages: (a) there is no X-ray limit since 
the electron current producing the X-radiation is proportional to pressure; 
and (6) the lowest measurable pressure is not limited by the vapor pressure 
of a hot tungsten filament, which is 10-" mm. Hg at 2300° K. (Alpert and 
Buritz 1954). Previous designs of Penning ionization-gauges suffered from 
two defects which made them unsuitable for the measurement of extremely 
low pressures: (a2) most designs of Penning gauge fail to strike at pressures 
below about 10-> mm. Hg; and (6) the ion-collecting electrode is exposed to 
high electric fields and the resultant field-emission from this electrode is 
indistinguishable in the measuring circuit from the positive-ion current. This 
field-emission current establishes a lower pressure limit analogous to the X-ray 
limit in the triode ionization gauge. The first defect can be overcome by a 
suitable geometric arrangement of electrodes which produces efficient electron 
trapping in the discharge region; such designs have been described by Beck 
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and Brisbane (1952) and Haefer (1953 and 1954). To overcome the second 
defect, previous designs have been modified by the addition of auxiliary- 
cathode electrodes which shield the ion collector from high electric fields and 
from which the field emission, required to start the discharge, is drawn. 

This paper describes a gauge designed with the above factors in mind, which 
has proved useful in the range of pressures from 10~* mm. Hg to the lowest 
pressures obtainable in this laboratory (about 1O~" mm. Hg). The structure 
of this gauge suggested that it be christened the “‘inverted-magnetron gauge”’ 
(1.M.G. hereafter). 

The theory and measurement of the striking characteristics of this gauge 
are described in a companion paper (Redhead 1958). 


2. DESCRIPTION OF GAUGE 


Fig. 1 shows a cutaway view of the basic arrangement of electrodes in 
the I.M.G. The ion collector is a close-ended cylinder with holes in the center 
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Fic. 1. Cutaway diagram of inverted magnetron gauge and associated circuit. 


of the end plates through which the anode and the shielding tubes pass. The 
auxiliary cathode is a box-like structure enclosing the ion collector and acting 
as an electrostatic shield. Welded through the sides of the box are two short 
tubular shields which project into the ion collector and through which the 
anode passes. Dimensions which have proved to give satisfactory performance 
are given in Table I. The shield cylinders project 2 mm. inside the ion collector 
and protect the end plates of the ion collector from exposure to high electric 
fields; they also provide the field-emission current which initiates the dis- 
charge. The current to the auxiliary cathode, consisting partly of the field- 
emission current, does not flow through the ion-current amplifier. 
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The anode was made of tungsten, all other electrodes were made of Nichrome 
V. All parts were rigorously hydrogen-fired before assembly; in the later 
gauges the parts were also vacuum-fired before assembly. 


TABLE I 
DIMENSIONS OF I.M.G. No. 16 


Dimension 


Electrode Symbol (mm.) 
Anode radius Ye 0.5 
lon collector radius rk 15 
lon collector length 2.5 


Shield cylinder radius 


3. OPERATING CHARACTERISTICS ABOVE 5X10-9 MM. Hg 
(a) Static Characteristics 

Measurements in the pressure range 10-* to 5X10~-* mm. Hg were made in 
a continuously pumped system by introducing dry air or helium through a 
variable leak. Pressure calibration was done with a Bayard—Alpert gauge 
(B.A.G. hereafter) of standard design. The B.A.G. was operated with the 
grid at +300 v., the filament at +45 v., the collector at ground, and a grid 
current of 8 ma. The sensitivity of the gauge under these conditions was such 
that the equivalent nitrogen pressure (in mm. Hg) was equal to five times the 
ion current (in amperes). Pressures quoted in this paper are equivalent 
nitrogen pressures. 

A standard operating condition for the I.M.G. was chosen with a magnetic 
field of 2060 gauss and an anode voltage of 6 kv. Over a wide range of pressure 
the collector current obeys the relationship 
(1) ty = cp", 
where p is the pressure and ¢ is a constant, and where a lies between 1.10 
and 1.4 for various gauges. The value of the exponent 1 is essentially independ- 
ent of the anode voltage, and approaches unity as the magnetic field is in- 
creased. The dimensions of the gauge have little effect on the value of x. 
Fig. 2 shows the current~—pressure characteristic for 1.M.G. No. 4 at a magnetic 
field of 2060 gauss, with dry air as the calibrating gas. The data of Fig. 2 are 
fitted by n = 1.10. 

The measured current-voltage characteristics at a magnetic field of 2060 
gauss are all of similar shape at pressures below 10~*mm. Hg. At these pressures 
the voltage-current curve rises steeply and approaches a constant-current 
characteristic at high voltages. Above this pressure a radical change in shape 
occurs, and the curve approaches a constant-voltage characteristic. Calcula- 
tions of the ratio of the transit-time of positive ions to that for electrons 
show that the net space-charge is negative for pressures below about 10~* 
mm. Hg, and positive for higher pressures (see Redhead 1958, Table I). Thus 
the measured change in the behavior of the current-voltage characteristic 
occurs at a pressure where the approximate theory indicates that the space 
charge changes from negative to positive. 
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Fic. 2. Collector current versus pressure. 


Measurements were made of the magnetic field and voltage necessary to 
give a particular collector current at a fixed pressure, and these constant- 
current characteristics plotted on the V4—B plane. The shape of these curves 
followed the shape of the striking characteristic and suggested that the relation 


(2) n/p = k(Va—Vs) 


might be approximately correct. Here V4 is the anode voltage and V, the 
striking voltage. Fig. 3 shows V,—V, plotted against 1,/p for the standard 
magnetic field (2060 gauss). At high magnetic fields the above relationship 
is seen to hold over a pressure range of four decades. This relationship has 
been found to be approximately correct for all gauges over a very wide range 
of pressure. The non-linear relation between ¢& and p for fixed V4 and B 
(equation (1)) occurs because V, is a slowly varying function of pressure 
(Redhead 1958, equation (19)). Equation (2) indicates the importance of 
keeping V, as low as possible in order to obtain a near-linear relationship 
between ion current and pressure. 

The sensitivity ratio of the gauge to helium and air was measured in the 
pressure range 10-4 to 10-* mm. Hg, and is shown in Fig. 4 for a magnetic 
field of 2350 gauss, and an anode voltage of 6 kv. This ratio decreases with 
pressure, and is a function of anode voltage. 
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PRESSURE .[mm. Hg] 
rs 108 
2 107 
r 10° 
2 107° 
v2 10% 


H = 2060 GAUSS 
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4 
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Fic. 3. Va— Vs versus 1%/p. 
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Fic. 4. Sensitivity ratio for helium and air as a function of pressure. 


(b) Oscillatory Behavior 

Oscillations were observable in the gauges under almost all conditions; over 
certain limited ranges of the parameters these oscillations were remarkably 
stable, both in frequency and amplitude. Some general comments can be 
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made about these oscillations in the regions where they were sufficiently stable 
to allow meaningful measurements to be made. 

The observed frequencies of these oscillations (in I.M.G. No. 16) varied 
from 400 c.p.s. to 150 ke./s., in the pressure range 10-7 to 10-* mm. Hg, and 
the frequency was independent of the impedance of the external circuit. Dis- 
continuous jumps in frequency from one mode of oscillation to another were 
observed, which coincided with small changes in the collector current. These 
mode changes account for some of the small discontinuities in the pressure— 
current characteristics that were observed. 

For fixed magnetic field and voltage the frequency increased linearly with 
collector current and was independent of the nature of the gas present. Measure- 
ments of frequency were made with helium and air at a pressure of about 
5X10-® mm. Hg as a function of collector current by varying anode voltage. 
For equal collector currents the frequencies observed for helium and air were 
equal to within the experimental error. 

Fig. 5 shows frequency measured on I.M.G. No. 16 as a function of anode 
voltage, magnet current, and collector current. Magnetic field was related to 
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Fic. 5. Frequency of oscillations as a function of anode voltage, collector current, and 
magnet current for [1.M.G. No. 16. 
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magnetic current in the following way: B(gauss) = 10.8 imag(ma.). These data, 
taken in the pressure range 10-' mm. Hg to 10-' mm. Hg with both air and 
helium, are a reasonable fit to the relation 


(3) f = k(Vat,/B)"'*. 

When the amplitude of the oscillations was large, electrons could gain 
sufficient energy to reach the auxiliary cathode. The electron current to the 
auxiliary cathode was quite large at higher pressures, and changed discon- 
tinuously as the coherent oscillations started or stopped. 

4. OPERATING CHARACTERISTICS BELOW 107-§ MM. Hg 


(a) Apparatus 
The apparatus (Fig. 6) used in this pressure range consisted of a sealed-off 
pyrex envelope containing an I.M.G., a standard B.A.G., a screen between 


Styrofoam Cap 


Immersed Finger 








Moveable Dewar 


Platform 








Fic. 6. Experimental apparatus for measurements below 1078 mm. Hg. 


the two gauges, and a finger extending down into a simple liquid helium 
dewar which could be raised or lowered on a movable platform. The enclosed 
volume was usually about 900 cc. and the inner diameter of the finger was 
2.3 or 1.1 cm. The over-all finger length was 30cm. Normally about 170 cc. 
of liquid helium were transferred. The screen was operated about 700 v. below 
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ground potential to prevent interaction between the two gauges, which was 
often found to be present even at the very low pressures used, in contrast 
to the results of Leck and Riddoch (1956). Ten sealed-off systems of this 
type, designated by the letters A, B, C, etc., were constructed and similar 
experiments performed on each. The system letter and corresponding I.M.G. 
design number are here given for future reference: A7, B10, Cll, A’ = A 
repumped 7, D7, Ell, F13, G15, H15, J15. 

The systems were evacuated by a mercury diffusion pump with a solid- 
CO, trap in series with a liquid-nitrogen trap between pump and system. The 
systems were baked for about 15 hours at 475°-500° C. Pressures between 
5xX10-" and 5X10-'® mm. Hg were obtained after seal-off as measured by 


the B.A.G. 


(b) Experimental Procedure and General Results 

The results obtained on all systems were not quantitatively the same but 
the procedures followed are described with reference to Fig. 7, which gives 
the results of a run on system A. 

The gauges were switched off overnight to permit the pressure to rise in 
the system by diffusion of atmospheric helium through the pyrex walls. The 
gauges were then switched on simultaneously and their collector currents as 
a function of time recorded with the whole system at room temperature (7; 
refers to collector current of I.M.G. or B.A.G.). The conductance of the 
tubing joining the gauges was approximately 37 liters/second for helium and 
the conductance of the I.M.G. end openings was greater than 12 liters/second, 
or about 500 times the I.M.G. electronic pumping speed for helium. Hence 
corrections for pressure drops between gauges could be neglected. It may be 
noted from Fig. 7 that the I.M.G. took a short time to strike. After subtraction 
of the B.A.G. X-ray limit (see Section 4(d)) the I.M.G. calibration curve 
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Fic. 7. Summarized record of experiment. 
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from 10-* to 10-!° mm. Hg was obtained. The results of this calibration for 
three I.M.G. designs differing in detail are given in Fig. 8. I.M.G. pressures 
quoted later have been obtained from the observed currents by means of 
these curves. 
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Fic. 8. 1.M.G. calibration curves below 10-§ mm. Hg. Magnetic field, 2060 gauss. Anode 
volts, 6 kv. 


Before any refrigerants were added, i.e., during the 2-hour period of Fig. 7, 
two measurements were made: 

(i) The room-temperature leak-rate was found by measuring the increase 
in pressure in the system after the gauges had been switched off for a known 
time (usually about 15 minutes). A typical record of I.M.G. collected current 
after gauge turn-on is shown in Fig. 9(a). The leak-rate C is given by 


(4) C= (Poeax — eq) /t, 


where Pyeax is the pressure corresponding tO Zpeax in Fig. 9(a), and Peg is the 
equilibrium pressure corresponding to iq. ¢ is the elapsed time the gauge was 
inoperative. With a new system the leak-rate was about 5X10-" mm. Hg per 
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minute, which is greater than the 3X10-'! mm. Hg per minute obtained by 
Alpert and Buritz (1954), probably because of the higher surface area to volume 
ratio of the systems used in the present experiments. The observed C rose 
as much as a factor of 10 a few days after the first use of liquid helium, which 
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Fic. 9. (a) Leak-up transient of I.M.G. at room temperature. (6) Leak-up transient of 
B.A.G. at room temperature. 


has been attributed to the increase of helium in solution in the envelope. For 
comparison, a similar leak-up record for a B.A.G. is shown in Fig. 9(0). It 
may be noted that there is no time delay between the application of voltages 
and the striking of the gauge, as there is with the I.M.G. This time delay, 
while often a disadvantage, permits a clear separation between the switching 
transient and the start of I.M.G. operation (see Fig. 9(@)). The characteristic 
electronic pumping time, zz, of either gauge could be measured from the 
decay of the pressure following the peak in Figs. 9(a) and 9(6). Typical values 
for tz obtained in these experiments were for the I.M.G., 0.7 minute, and 
for the B.A.G., 4 minutes; corresponding to electronic pumping speeds of 
0.02 and 0.004 liter/second respectively for helium. The curves of currents 
versus time after gauge turn-on shown in Fig. 9 have been termed ‘‘leak-up 
transients at room temperature’’. 

(ii) A collector-current versus grid-voltage characteristic for the B.A.G. 
was measured after the manner of Alpert and Buritz (1954). The results are 
discussed in Section 4(d). 

Following these two measurements liquid nitrogen was added to both 
dewars (to precool the inner dewar), which caused a pressure drop of as 
much as 30°% in some systems. An unknown portion of this drop may be 
assigned to the elimination of helium diffusion through the finger (Rogers, 
Buritz, and Alpert 1954). When the liquid nitrogen was removed from the 
inner dewar two peaks in collector current were generally observed and were 
assigned to adsorbed gases. 

The transfer of liquid helium caused both collector currents to drop sharply. 
In the case shown in Fig. 7 the B.A.G. current fell to 3107!!! amp., corre- 
sponding to 1.5X10-!® mm. Hg, and taken as the X-ray limit of the B.A.G. 
(see Section 4(d)), while the I.M.G. (No. 7) current fell to 8X10~-" amp., 
corresponding to a pressure of 2.6X10-" mm. Hg (from Fig. 8). At I.M.G. 
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currents less than 10~-" amp. the signal-to-noise ratio was about 2:1. Fre- 
quently, but not in Fig. 7, the I.M.G. current with liquid helium present 
fell by about a factor of two during a 2} hour immersion period with no 
change in B.A.G. current. 

While liquid helium was present the leak-rate measured by switching off 
either gauge for a known time was always zero, demonstrating that the pump- 
ing action of the gauges was not sensibly influencing equilibrium pressures 
and that all gases, in particular helium, were being continuously adsorbed 
on the cold finger. 

A repetition of (ii) above with liquid helium present gave results discussed 
in Section 4(¢). When the liquid helium level fell below the tip of the finger 
a rise in currents, usually fast, was observed (Fig. 7). This pressure rise was 
caused by gas, mainly helium, desorbing from the finger, followed by electronic 
pumping of the desorbed gas. This desorption transient, termed the “d-tran- 
sient’, could be produced at will by lowering the dewar at any time during 
immersion. For a d-transient with a steep leading edge a d leak-rate could 
be found from the peak height less the equilibrium pressure divided by the 
elapsed time the finger had been immersed in liquid helium, in a manner 
analogous to equation (4). 

(c) Lack of Reproducibility of I.M.G. Collector Currents with Liquid Helium 
Present 

The collector currents with liquid helium present showed essential repro- 
ducibility in B.A.G. currents, but large divergencies in I.M.G. currents. Fig. 10 
shows the equilibrium values of I.M.G. and B.A.G. currents, measured 
simultaneously, first at room temperature and a short time later with liquid 
helium in the dewar, for all immersions for which these four currents were 
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Fic. 10. Equilibrium I.M.G. collector currents at room temperature, and later with liquid 
helium present, for systems A-G. 
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measured. Conversion of currents to pressures with the data of Fig. 8 reduces 
but does not remove the essential scatter. The 15 immersions shown were 
performed with eight different systems. The ratio of pressure at room tem- 
perature as measured by the I.M.G. ranged from 1.5 to 150 over all experiments 
performed, both these two extremes being obtained with the same system. 

Most of these scattered results were in the pressure range below the B.A.G. 
X-ray limit and hence the B.A.G. could not be used as an independent measure 
of pressure. The question naturally arises whether the 1.M.G. provides a 
reliable measure of pressure in this range. An investigation was therefore 
undertaken to provide more details about the observed variations in I.M.G. 
currents with liquid helium present. 

It was found that high I.M.G. currents with liquid helium present resulted 
when an immersion followed a heavy outgassing of the B.A.G. or severe heating 
of the system. Low I.M.G. currents with liquid helium present were observed 
after the gauges had been turned off or had operated normally for many 
hours. For example, the two immersions with system G shown in Fig. 10 
were performed on the same day and differed only in that a heavy outgassing 
of the B.A.G. (250 ma. d-c. bombardment current at 800 v. a-c.) for } hour 
was carried out between the two immersions. The second immersion took 
place 1.3 hours after the return of the B.A.G. to normal operation. In Fig. 11 
are given results illustrating the same point. 

Further detailed information on the processes involved in these experiments 
was provided by the three immersions labelled E1, E2, E3, of Fig. 11, and 
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Fic. 11. Equilibrium collector currents at room temperature, and later with liquid helium 
present, for systems E, H, J. 

1. 46 hours after heavy B.A.G. outgas. 

2. 7 days after 1 but with no B.A.G. outgas in the interim, but with B.A.G. operating 
intermittently. 

3. 4 hours after 2 with 1 hour heavy outgas of B.A.G. in interim. 

4, 22 hours after seal-off, no B.A.G. operation. 
5. 45 hours after seal-off, no B.A.G. operation. 
6. 5 days after seal-off, no B.A.G. operation. 
7. 47 days after seal-off, no B.A.G. operation. 
8. 23 hours after seal-off, no B.A.G. operation. 
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by a fourth immersion with the B.A.G. operated alone in the same system 
a few hours before E2. In particular the results of El and E3, which yielded 
high I.M.G. currents with liquid helium, are to be contrasted with those of E2, 
which yielded a low I.M.G. current in the presence of liquid helium. 

It has been shown by Alpert and Buritz (1954) that if the equilibrium 
between helium diffusion through pyrex and electronic pumping determines 
the equilibrium pressure P’.g (mm. Hg), then 


(5) Pe/ta'C’ = 1, 


where rz’ is the characteristic electronic pumping time for He (minutes) 
= V/S,’ (here V is the volume of the system in liters, Sz’ is the electronic 
pumping speed for helium in liters/minute), and C’ is the leak-rate of helium 
into the system in mm. Hg/minute. For the runs mentioned above Peg, Tz, 
and C were measured as described in Section 4(0)(i) and equation (5) checked 
with the results shown in Table II. 


TABLE II 
RESULTS ON SYSTEM E 





Prec TE 
Gauge (mm. Hg) (minutes) (mm. Hg/minute) Peq/tTEC 
1.M.G.(E1) LE x<ie-* 0.70 2.7X107- 15 
B.A.G. 2.4X10-9 4.4 4.6107" 1.2 
1.M.G.(E2) 3.1XK10- 0.80 4.2107 0.91 
1.M.G.(E3) 7.8X107- 0.77 5.5 1071 1.9 





For runs El and E3, equation (5) is not satisfied, while for the other two 
runs it is satisfied for both gauges within experimental error. It may be noted 
that C does not vary much from run to run. 

A comparison was also made between the room-temperature leak-up 
transients and the d-desorption transients. For each run one leak-up transient 
at room temperature and four d leak-up transients were observed. These five 
transients are arranged in horizontal rows across Figs. 12(a) and 12(d), those 
of Fig. 12(a) being room-temperature transients and those of Fig. 12(d) being 
d-transients. Each row corresponds to a different run, the top row with the 
B.A.G. operating alone and the second and third rows with the I.M.G. 
operating alone during runs E2 and E3 respectively. In Fig. 12(a) the bracketed 
figures represent the time in minutes the gauge was not operating. In the first 
and third d-transients of any run the gauge was operating when the dewar 
was lowered, while in the second and fourth d-transients the gauge was off 
when the dewar was lowered and was turned on a short time later. For the 
d-transients the first figure in brackets represents the time the finger was 
cold, and the second the time between the lowering of the dewar and the 
start of operation of the gauge (both in minutes). For either gauge all tran- 
sients have a characteristic shape except for two desorption peaks during 
run E3 when the dewar was lowered with the I.M.G. operating. In these 
cases a sharp spike is superimposed on the normal background. If the spikes 
are disregarded (dotted line in Fig. 12(0)) the results become the same as 
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those of E2. The tendency for the first d-transient to be lower than the others 
may be due to an initial irreversible adsorption of gas (Armbruster and Austin 
1944). If the spikes and the smaller desorption transients are disregarded, then 
the room-temperature and d leak-up transients are essentially the same. The 
presence of the spike causes the d leak-rate to exceed the room-temperature 
leak-rate, a result observed many times on other systems on both the I.M.G. 
and B.A.G. In particular, during run El the ratio of d leak-rate to room- 
temperature leak-rate was a factor of four. An interesting feature of the 
spike is that it is not seen if there is a delay of 2 minutes between the lowering 
of the dewar and the start of operation of the gauge, suggesting pumping 
of the gas in the spike by other than electronic means. 

The results of the B.A.G. run and E2 have been adopted as characteristic 
of a system with equilibrium pressures determined by helium diffusion (termed 
‘process I’’), and those of El and E3 as indicating the presence of a second 
and essentially different process (termed ‘process II’’). Process I has the 
properties: (1) that Peq = 7eC, and (2) that d-desorption transients are 
identical to room-temperature leak-up transients; while process II has the 
properties: (1) that Peg > 7eC, and (2) that d-desorption transients can have 
a sharp spike superimposed on the room-temperature leak-up transient. Pro- 
cess IT gives high I.M.G. currents with liquid-helium trapping and the presence 
of this process in varying degree is taken to account for the scattered results 
of Fig. 10. Further investigation is being carried out on process II, and it is 
hoped that results will soon be available. 

In the present paper discussion will henceforth be limited to those cases 
where no evidence of process II was observed. Such cases are now discussed. 


(d) Results for Systems with Equilibrium Pressures Determined by Helium 
Diffusion 

Fig. 13 shows two plots of the type described by Alpert and Buritz (1954) 
for B.A.G. collected current versus voltage between grid and filament. The 
upper curve was taken at room temperature and the lower with liquid helium 
present, a short time later. There is no doubt that the addition of the liquid- 
helium trap reduces the pressure, and the question arises whether any estimate 
of pressure with liquid helium present may be made from the B.A.G. Accord- 
ingly the dashed straight line was taken as the X-ray contribution, which was 
subtracted from the observed results to yield the dashed curves. The upper 
dashed curve has the familiar shape of an ionization-probability curve and 
the lower dashed curve does not. While the detailed shape of the lower dashed 
curve was not quantitatively reproducible either in the same gauge or in 
different gauges, its general shape was always different from the ionization- 
probability curve. It was concluded that no quantitative estimate of pressure 
in the presence of liquid-helium trapping could be made with the B.A.G. The 
current at the normal operating grid voltage of 255 appeared to be exactly 
on the X-ray line and was taken as the X-ray limit of the B.A.G. This corre- 
sponded to a pressure between 1.0X10-!° and 1.5X10~-'!’ mm. Hg for the 
gauges used, and the results of Fig. 13 indicate that the true pressure with 
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liquid helium present is at least an order of magnitude below this, placing an 
approximate upper bound on the true pressure at 10~!! mm. Hg. 
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Fic. 13. Collected current versus grid voltage for Bayard—Alpert gauge, with grid current 
held constant. 


A lower bound may be found by comparison between an ideally calculated 
pumping speed of the liquid-helium trap (S7) and a measured electronic 
pumping speed of the I.M.G. for helium (Sz). The expected pressure with 
liquid helium present (Pye) may then be obtained from the measured pressure 
with liquid nitrogen in the trap (Px) from: 


(6) Pue = PySz/(Sr+Sz). 


In the calculation of S; it has been assumed that the pressure inside the 
finger immersed in liquid helium is zero and hence Sz is just the conductance 
of the trap tubing. Seventeen centimeters of tubing were taken at room 
temperature, and 8cm. inside the dewar above the liquid helium level at 
liquid nitrogen temperature, with the remaining 5cm. immersed in helium 
and hence not contributing to the conductance. For A, A’, G, H, the tubing 
inner diameter was 2.3 cm., while for E it was 1.1 cm. 

The values of Py. obtained from (6) with S7 calculated as indicated above 
are shown in Fig. 14 along with the observed pressures as given by the I.M.G. 
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for nine immersions. With one exception, where the deviation is not serious, 
the I.M.G. observed pressure lies between the approximate upper bound of 
10-" mm. Hg and the calculated lower bound. 
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Fic. 14. 1.M.G. results with liquid helium trapping in the absence of process II. 


It is possible that the true pressure approaches the lower bound, and a 
residual current in the I.M.G. analogous to the X-ray current of the B.A.G. 
limits the use of the I.M.G. to a range above 10-" mm. Hg. Measurements 
rule out the possibillty that an ohmic leakage is responsible for this residual 
current. 

It is also possible that a true pressure of about 10~-” mm. Hg exists above 
10-* to 10-> monolayers of helium adsorbed on pyrex at 4.2° K., and that the 
I.M.G. is giving a true measure of this equilibrium pressure. The thermo- 
molecular pressure correction of 0.12 tends to be compensated by a factor 
of about 0.12 (Fig. 4, and see also Alpert 1953) for the ionization probability 
of helium as compared to nitrogen, and hence equivalent nitrogen pressures 
in the gauge at room temperature are numerically almost equivalent to true 
helium pressures in the finger at 4.2° K. While the helium adsorption isotherm 
has not been measured at coverages as low as those found in these experiments 
it is suggested by Meyer (1956) that more helium is adsorbed at low pressures 
than is predicted by an extrapolation of the Brunauer-Emmett-—Teller theory. 
Garfunkel and Wexler (1954) place an upper bound of 107!° mm. Hg for 
helium on brass at approximately 1.4° K., probably with higher coverages 
than those used here. Thus a true helium equilibrium pressure of 10~-” mm. 
Hg in the finger at liquid-helium temperature does not appear impossible. 

Thus the I.M.G. appears to give a valid measure of helium pressure down 
to 10-” mm. Hg. Experiments with liquid-helium trapping in pyrex apparatus 
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of the size used cannot be expected to produce pressures below 107!’ mm. 
Hg. Whether the lower gauge limit is instrumental or a true pressure is still 
to be decided. 


5. CONCLUSION 
The inverted-magnetron gauge has been shown to be useful in the pressure 
range 10-* to 10-" mm. Hg. The upper limit is set by the radical change in 
characteristics that occurs when space charge becomes positive. The lower 
limit has not been determined with certainty and may extend below 10~" mm. 
Hg. The sensitivity of this type of gauge is an order of magnitude greater 
than that of the Bayard—Alpert gauge. 
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THE SELF-CONSISTENT FIELD WITH EXCHANGE FOR NEON 
BY FERUT PROGRAM! 


BEATRICE H. WorSLEY 


ABSTRACT 


A general program for the computation of atomic self-consistent fields with 
exchange has been prepared for the University of Toronto electronic digital 
computer, FERUT. It was developed by application to neutral neon, for which 
the results are given, and tested against the independently calculated results 
for Fe*4 and Fet'®. 


INTRODUCTION 

The problem of defining atomic structures is a difficult one to reduce, even 
approximately, to practical terms. The usual treatment consists of applying 
the variational principle to the Schrédinger Hamiltonian, with an approxima- 
tion to the wave function for the many-electron system in the form of a 
determinant of one-electron wave functions. The problem can then be reduced 
to the solution of a system of ordinary differential equations of the second 
order. These equations have been named after their discoverers, Hartree and 
Fock. 


The Hartree-Fock equations are of the general form 
(1) P"' (nl; r) + g(r) P(nli r) +f(r) = 0, 


for the set of functions P(n/;r) defining the radial wave functions, there 
being one such function for each electron group in the configuration. The 
group is defined by the usual quantum numbers » and /. For each function 
P(nl; r), the corresponding function g(r) is known as the potential term, and 
it is in itself a function of all the wave functions in the set as well as of a 
Lagrange multiplier €,;,.; which defines the energy of the n, / electron group. 
For each function P(nl;7r) the corresponding function f(r) is known as the 
exchange term, and this is also a function of all the wave functions in the 
set as well as of a set of Lagrange multipliers e€,;,,,7 which define the exchange 
energies between electrons in the 2,/ group and electrons in all the other 
groups n’, /’. Each solution P(nl; 7) is subject to two-point boundary as well 
as orthonormalizing conditions. The non-linear nature of the Hartree-Fock 
equations as well as the above-stated conditions renders them difficult to 
solve. The only known method is a tedious numerical process involving 
successive approximations developed by Hartree (1927). It is based on the 
principle that the field of the average electron distribution derived from the 
wave functions must be the same as the field used in evaluating these wave 
functions. The atomic field so determined is therefore known as the “‘self- 
consistent field” (S.C.F.). It is, except for ions with very few electrons, prob- 
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ably the best theoretical description of atomic energy levels that can be 
obtained. 

The importance of the one-electron atomic wave functions has long been 
recognized. Up to the end of the Second World War, calculations had been 
carried out by hand for most of the light atoms up to argon, and for several 
atoms in the first long period. Applications continue to be calculated by 
hand wherever the results are urgently needed, although each configuration 
may require from three months to over a year. An early attempt at mechanizing 
the process was made by Hartree (1934), applying a differential analyzer to 
Hg and Hgt+. This machine gave graphical results of low accuracy, and 
required considerable human intervention and ingenuity. 

Recent interest in atomic structures arising in the fields of solid-state 
physics, crystallography, and astrophysics has coincided with the development 
of high-speed digital calculating machines. For these machines a program 
can be prepared once and for all, so that much of the tedium is removed 
from individual applications. 

Such programs have been prepared independently for the EDSAC I Com- 
puter of the University of Cambridge by the author (Worsley 1952), by 
Douglas (Douglas et al. 1955; Douglas 19560), by Ridley (1955), and by 
Froese (1957) ; for the Manchester University Computer Mark II by Altmann 
(1955); for the IBM 650 Computer by Piper (1956), and by Treanor of the 
University of Buffalo; and, as described here, for the University of Toronto 
computer FERUT, by the author. Considerable variety has been shown in 
these programs in many particulars. These include such aspects as the extent 
of applicability, the amount of judgment coded into the program, the for- 
mulation of the equations for evaluating the potential and exchange terms 
as well as of the Hartree-Fock wave equation, the choice of independent 
variable, the selection of the intervals of integration, the technique for satis- 
fying the conditions of orthonormality simultaneously with the boundary 
conditions in solving the wave equation, and finally such details as the numeri- 
cal representation, scaling, and decimal-point location. Some programs have 
not allowed for the inclusion of the exchange terms, and are therefore of more 
limited application. However, they all follow the basic Hartree procedure 
described in such papers as Hartree and Hartree (1936). A concise summary 
of the problem and of the applications to that date has been published by 
Hartree (1948). Many more aspects of the entire subject of calculating atomic 
structures have more recently appeared in book form (Hartree 1957). 


THE FERUT PROGRAM 


The entire procedure for applying the FERUT S.C.F. program to a given 
configuration is illustrated in Fig. 1. To begin with, it is necessary to write 
out the detailed Hartree-Fock equation for each wave function in the set. 
This amounts to listing the set of coefficients which define the potential and 
exchange functions in terms of the wave functions. This information must be 
transcribed into input tapes which FERUT can interpret through the main 
S.C.F. program. These are termed C tapes. A tape for the main program 
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has been prepared once and for all. Then it is necessary to prepare a set of 
numerical estimates to the solution. For each wave function this amounts to 
preparing a tape containing the wave function tabulated at the intervals of 
integration, together with estimates of the starting value [P(r)/r'*4Jo, = A 
say, the potential energy €,;,n:, and the set of exchange energies €;.,/). These 
are termed P tapes. 





of improved 
and print the 





Fic. 1. Flow diagram for the FERUT S.C.F. process. 


Auxiliary programs are provided to facilitate preparation of the C and P 
tapes. Also, since the S.C.F. process is essentially one of successive approxi- 
mation, the time required to refine a given set of estimates to some prescribed 
degree is very greatly dependent on the accuracy of the initial estimates. 
Methods for obtaining them are summarized in the Hartree text (1957) and 
much work has been done in this direction on the EDSAC I and FERUT 
computers by Froese (19576). 

The sequence in which the wave functions to be improved are selected 
is left to the operator’s judgment. It is usually best to refine the functions 
for the outer groups before disturbing the core. Also, it has been found that 
convergence is more rapid and direct for more highly ionized atoms. 

In its basic form, the main S.C.F. program with exchange has been devised 
so as to deal with all atoms in neutral, excited, or ionized states with certain 
restrictions. The total number of wave functions is limited to 14, of points 
in the range of integration to 117, of terms in the potential function* to 8, 





*Which are not of the form 2(2/-++1) Yo(nl, nl; p). Cf. pe 293. 
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of terms in the exchange function to 13, and of addends in each term of 
the exchange function to 5. Provision has been made for introducing certain 
variations into this basic form with ease. These include extension to more 
than 14 wave functions, neglect of exchange, selection of parameters to improve 
convergence, choice of a tractable approximation to © for the upper bound 
to the independent variable, and introduction of short cuts when testing, 
without improvement, the accuracy of a given set of estimates. 

The operations performed by a single run through the main S.C.F. program 
(boxes III to VI of Fig. 1) are illustrated in greater detail in Fig. 2. The 
subscript 7 is used to define the wave function and associated parameters 
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Fic. 2. Flow diagram of the main S.C.F. program. 


under improvement. The subscript j defines stages in the A-improvement 
cycle within the main program. Also, in practice, it is usually necessary to 
test the orthogonality of the wave function being improved with one other 
wave function. The subscript 0 is used to define this latter function. The 
independent variable has been chosen as p = log,cr, where ¢ and the interval 
6p are constants which may be treated as parameters of the program. 
With reference to Fig. 2, when tape input is complete (boxes I, I], and 
III), the program first calculates the potential function g;(p) and the exchange 
function f;(p) (boxes V and VI) for the wave function being improved. These 
functions are constructed out of functions Y;.(ml), melo; p) which are in turn 
solutions of differential equations of the second order satisfying two-point 
boundary conditions and having functional coefficients dependent on pairs of 
wave functions, P; and P», in the configuration. The k’s are small positive 
integers. The equations for Y, have been reduced, by means of an auxiliary 
function Z;, to a pair of first-order equations each having a one-point boundary 


condition, namely to 
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(2) dZ,/dp = rP\P2—kZ,, 
where the integration is carried outwards from Z,(0) = 0, and 
(3) dY,/dp = (kR+1) Y,—(2k+1)Z;, 


where the integration is carried inwards from Y,(©) = Z,(). As a refine- 
ment, all the Y functions of the form Yo(nl, nl; p) for nl ¥ n,l; are con- 
structed first and stored (box IV), so that they need not be recalculated in 
later recycling to improve the same wave function. 

The Hartree-Fock equation is next solved for P;, the wave function being 
improved (box VIII). This equation (1) has been rewritten in the form: 


(4) d°(P,r—'"*)/dp? = {r?(gitenini1) 1/4} (Pir) +r?/F;, 
where €.7,n1,; Must be chosen so that P;(p) ~ 0 as p— ~, and A; so that 
\ %> - 1/2 
(9) ) rP; dp = N 
ei 


approaches unity. 

The scheme for adjusting both e; and 4; will now be described (boxes VII, 
IX, X, XI, and XII). It has proved to be practicable and automatic, con- 
verging without operator intervention. It is well known that, if only an 
outwards integration is carried out, any solution of a Hartree-Fock equation 
diverges drastically to + after the asymptotic condition for large r has 
been satisfied. Also a value of ¢€, say e+, just above the true eigenvalue of 
(1) leads to a solution with (7—/) nodes, whereas a value of €, say €, just 
below the true eigenvalue leads to a solution with (7—/+1) nodes. The true 
solution has, of course, (7—/) nodes. Only outwards integrating runs are 
employed in this program. The e-adjustment to satisfy the condition P(p) > 0 
as p— ©, which takes place at the end of each run (box IX), is performed 
in two stages. In the first stage, the trial € is adjusted by an increment of 
preassigned magnitude +Ae, the sign being selected by the program according 
as divergence takes place to *@. As soon as values of € of each type, namely 
e* and e-, are discovered, the second stage of the process begins. This is a 
binary or averaging process. The improved e€ is taken to be (e++e~)/2. The 
two most recent values of «+ and e&~ are employed here, and their average is 
later used to replace e+ or e€~ according as it proves to be of the e* or & type. 
The stage II process continues until two successive values of e€ differ by less 
than a preassigned amount e,. When this occurs, an integrating run is made 
with the most recent value of e~, and the resultant P is then used to evaluate 
N= if, P2dr}?. It was found that with solutions of the e~ type, the best 
results were obtained when the outwards integration was carried out only 
as far as the (n—/+1)th node.* The trial values of P are used for the final 
values beyond this node, and checks are made that this approximation does 
not affect the accuracy of the orthonormalizing processes. 


*It is likely that an e* type solution would be equally satisfactory, provided that the out- 
wards integration were carried out only to the point beyond the (m—/)th node at which the 
slope of P(p) changed sign. However, this contingency was not so easy to detect automatically, 
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The A-adjustment is made by utilizing the close linear relationship which 
exists between good approximations to A, e, and V. After the first adjustment 
of € to determine the e; and V, which correspond to the trial value of A = Aj, 
a new trial value for 4, namely 4y, = A;/Njy, is applied and the e-adjustment 
process is repeated. Values «, and Vy are thus obtained. Ayy and ey are 
then calculated from the formulae 


Am = {(Ni—1) An —(Nu—-1) Ay} (Vi—Nn) 


and 


a= {(Mi- Leg —(Nu— Le} /(Vi-— Nin), 
which assume that the relationships 


Ay _ An _ Am a _ 1 _ an 

N,; Nu Min’ Ny Ny Nin’ 
hold. By a third application of the e-improvement process, this time with 
+1/100Ac, the converged value of ey, is determined. This value and Ay; 
are used in the final run to determine P, and a check is made that Vy is 
equal to 1 to the desired degree of accuracy. Note that the resultant P func- 
tion is never actually multiplied through by a normalizing factor. 

When the program has thus selected values of A; and ¢;, it calculates the 
corresponding value of P;, storing this as the “new”’ wave function in standard 
P tape form (box XIII), and differencing it against the estimate of the 
function input. Wherever necessary, a check is also made that this new 


and Nur = 1 


function is orthogonal to one other, say to Po, i.e. that 


(6) j rP}"P, dp=0 
v0 


tends to zero (box XIV). The operator must then decide whether or not to 
repeat the run with an improved estimate of the exchange energy €n,i;,n)%- If 
not improving €»; 7; 25%, the new wave function is next weighted with the trial 
wave function by some prescribed factor in a linear manner (box XV) and 
the resultant function stored and output as the improved estimate. The 
operator must then decide whether to recycle for the same wave function, 
starting now with the improved estimates of P;, €.:.n1;::, and A;, or to proceed 
to improve the next wave function (boxes XVI to XIX). 


NUMERICAL TECHNIQUES 
Integrating routines were specially developed for FERUT with a view to 
speed and accuracy. To integrate equations (2), (3), (5), and (6) the following 
process (Hartree 1956) was found to yield 5 or 6 steps per second, meanwhile 
ensuring 5 or 6 significant decimal figures and remaining stable with increasing p: 
To solve 


dy/dx+ky = f(x) 


introduce the integrating factor e*”. This gives 


de*y/dx = ef (x). 
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Upon integrating from x, to xj41 = xj+6x, one gets 


er 


j+1 
-kér ) kb k(z—-z;) 
Vor =e +e | e* 7) (x) dx. 
e 


zj 
By applying Simpson’s rule to the integral in the latter equation, and by 
putting 


Sing = Cfitfiar)/2— (OF) +8 f 41) /16 and 8fy = fj41— 2 Afi 
the following formula* can be derived: 
(7) vygn = 9 Oe ax 54.2/24+ {1/6496 YY /24) of 541 


+ {eh 24-4 6 /6) dxf — oY Oxf ,_1/24. 


To integrate (1), which is of the type y’’ = Fy+G, the Numerov formula 
(Numerov 1933) was applied in the form 


(8) Vier = (2yj—yy-1+ (6x) "[5F yy, /O+ Fwy 1/12 +G 541/12 
+5G ;/6+G;-1/12]} + {1—(6x)°Fy41/12}. 


Accuracy to five or six decimal figures was obtained without recourse to the 
refinement due to Olver (Hartree 1955). The speed is about two steps per 
second for the second-order equation. 

Starting values for equation (4) are calculated by the main program from 
the approximating series for P and equation (4) itself. They also require 
the current value 4 ;;, and constants generated by an earlier part of the main 
program and by the auxiliary program TC : PREP/B. Finally, whenever the 
independent variable appears explicitly, it is generated by repeated multi- 


plication. 


CODING NOTES 

In FERUT (Hume 1954) numbers are normally represented in binary 
form, each containing 40 bits and lying in absolute value between +1/2 
and —1/2. Furthermore, all arithmetic operations hold the binary point 
fixed. It has therefore been necessary to code in alternative numerical repre- 
sentations and operations such as scaled, truncated, and floating point, at 
various stages of the calculation. Were it not for limitations of speed and 
storage capacity of the machine, floating point could have been used through- 
out. Assignment of storage space to the various parameters, functions, and 
intermediate results has also been complicated by the physical structure of 


*Other methods investigated in this connection were: 

(i) The Runge Kutta Process modified by Gill (1951), which was very slow (two steps per 
second) and required tabulation of auxiliary functions at the half-interval. 

(ii) The Milne Process (1953), which developed a troublesome oscillating error as p increased. 

(iii) A one-dimensional Liebmann process to solve the two-point equation for Y, directly. 
This was much too slow, as of the order of 800 iterations were required over the entire range 
to obtain 4-decimal-figure convergence. Also convergence could take place to the wrong 
values when round-off errors were not kept to a minimum (Worsley 1952). 
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FERUT. The internal auxiliary magnetic drum store is of large capacity, 
holding the equivalent of over 16,000 standard numbers, or “words’’, but 
the high-speed store, which must share numbers as well as machine instruc- 
tions, holds only 256. To allow as many as 117 points in the tabulation of 
p it was not feasible to take advantage of the matching technique (Douglas 
1956a) for adjusting €,;,.,, in the integration of (4). 

All data, instructions, and intermediate results are stored on the magnetic 
drum at some stage of the calculation. Words are stored in tracks around the 
drum, each track containing 64 words. P tables require p tracks (Pb = number 
of wave functions in the set), C tables at most 33 tracks, machine instructions 
occupy 263 tracks,* and (22+2m) tracks are required for intermediate data, 
where m = number of terms in the exchange function. 

Several weeks are required to write out the wave equations, calculate 
estimates, prepare C and P tapes, and also to decide on the various runs 
necessary through the main program. However, the time to execute the main 
program on the machine is relatively short. One complete run to improve a 
specified wave function requires from 30 minutes to an hour, depending on 
the total number of functions in the system and the number of steps in the 
range of integration. If following up with an improvement run on the same 
wave function, about 10 minutes can be cut out of this time by by-passing 
read-in and evaluation of the unchanged Yo(n/, nl) tunctions. Also, if the 
value of V; happens to be within a preassigned tolerance of unity, a by-pass 
occurs automatically to cut out the next two .1-adjustment steps, thus saving 
a further 10 or 15 minutes. This latter is likely to occur eventually in the last 
one or two improvement cycles of each wave function. The main program 
is supplied with seven rescue points, so that in event of mechanical derail it 
is not necessary to restart the calculation at the beginning. Other console 
controls are coded in to facilitate such optional procedures as check printing, 
short cutting, and recycling. 

RESULTS 

The FERUT program results for neutral neon are presented, in atomic 
units, in Table I. The wave equations have been published in detail in the 
Hartree text (1957). The successive approximations were carried out until 
the maximum variations in the P’s, «l, and e€ were 0.09, 0.001, and 0.01% 
respectively. The maximum errors in the normalizing and orthogonalizing 
constants were then 0.02 and 0.05% respectively. It was not found necessary 
to adjust the value of e,.25 from its zero trial value. The results presented 
therefore satisfy the formulation to within two or three units in the least 
significant figure quoted. They differ by as much as 10° from the initial 
estimates, which were originally supplied by Prof. Hartree, and required 
some 15 hours of machine time, exclusive of check runs. 

Some early hand calculations on neutral neon were published by Brown 
(1933). These are of the self-consistent type, the equations having been 


*About 3400 machine instructions, since there are two instructions per word. 
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TABLE | 
NORMALIZED RADIAL WAVE FUNCTIONS WITH EXCHANGE FOR 
NEON (Is)*(2s)?(2p)® CALCULATED WITH INDEPENDENT VARIABLE 
p = log,.d500r, 5p = 0.1 





P(1s) P(2s) P(2p) 
Entnl 65.55 3.866 1.705 
A = (P/r'*1)o 60.77 14.27 27.87 
r P(1s) P(2s) P(2p) 
0.0020000 0.1191 0.0279 0.0000 
22103 0.1313 0.0308 0.0001 
24428 0.1448 0.0339 0.0001 
26997 0.1596 0.0374 0.0001 
29836 0.1759 0.0413 0.0002 
0.0032974 0.1938 0.0455 0.0002 
36442 0.2135 0.0501 0.0003 
10275 0.2350 0.0551 0.0004 
44511 0.2586 0.0607 0.0005 
49192 0.2845 0.0668 0.0006 
0.0054366 0.3128 0.0734 0.0007 
60083 0.3437 0.0807 0.0009 
66402 0.3775 0.0886 0.0011 
73386 0.4143 0.0972 0.0014 
81104 0.4544 0. 1066 0.0017 
0. 0089634 0.4979 0.1168 0.0021 
99061 0.5451 0.1279 0.0025 
0.010948 0.5963 0.1398 0.00381 
12099 0.6514 0.1527 0.0038 
13372 0.7109 0.1666 0.0046 
0.014778 0.7747 0.1815 0.0056 
16332 0.8431 0.1974 0.0068 
18050 0.9159 0.2144 0.0082 
19948 0.9933 0. 2323 0.0100 
22046 1.0751 0.2512 0.0121 
0.024365 1.1612 0.2710 0.0146 
26927 1.2511 0.2916 0.0176 
29759 1.3443 0.3128 0.0212 
32880 1.4404 0.3344 0.0256 
36348 1.5383 0.3561 0.0307 
0.040171 1.6370 0.3777 0.0369 
$4396 1.7352 0.3986 0.0441 
19065 1.8314 0.4185 0.0527 
54225 1.9237 0.4366 0.0629 
59928 2.0100 0.45238 0.0748 
0.066231 2.0881 0.4647 0.0887 
73196 2.1555 0.4730 0.1049 
80895 2.2094 0.4762 0.1237 
89402 2.2473 0.4731 0.14538 
Y8805 2.2664 0.4626 0.1701 
0. 10920 2.2644 0.4435 0.1984 
. 12068 2.2391 0.4147 0. 2303 
. 13337 2.1892 0.3750 0. 2662 
14740 2.1137 0.3236 0.3060 
16290 2.0130 0.2597 0.3500 
0. 18003 1.8883 0.1831 0.3979 
. 19897 1.7421 0.0939 0.4495 
21989 1.5781 —0.0073 0.5043 
. 24302 1.4011 —0.1192 0.5616 
. 26858 1.2167 —0.2398 0.6206 
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NORMALIZED 


TABLE I (Concluded) 


RADIAL WAVE 


FUNCTIONS WITH 





1958 


EXCHANGE FOR 
NEON (1s)%(2s)?(2p)® CALCULATED WITH INDEPENDENT VARIABLE 
p = log.d500r, 5p = 0.1 








P(1s) P(2s) P(2p) 
habs 65.55 3.866 1.705 
A =(P/r'*), 60.77 14.27 27 .87 
r P(\s) P(2s) P(2p) 
0.29683 1.0313 —0.3661 0.6799 
32804 0.8513 —0. 4948 0.7382 
36254 0.6825 ~0.6216 0.7937 
40067 0.5301 0.7419 0.8446 
44281 0.3977 —(). 8508 0.8890 
0.48938 0.2874 —0.9486 0.9249 
54085 0.1995 —1.0160 0.9506 
.59773 0.1325 — 1.0646 0.9646 
. 66060 0.0840 — 1.0870 0.9657 
. 73007 0.0506 — 1.0825 0.9535 
0. 80686 0.029 —1.0517 0.9280 
89172 0.014 —0.9969 0.8897 
98550 0.008 —0.9216 0.8398 
1.0891 0.003 —0.8304 0.7799 
1.2037 0.001 —0.7287 0.7122 
1.33803 —0.622 0.6389 
1.4702 —0.516 0.5624 
1.6248 —0.414 0.4852 
1.7957 —0.322 0.4096 
1.9845 —(0.242 0.3377 
2.1933 —0.174 0.271 
2.4239 —0. 120 0.212 
2.6789 —0.080 0.161 
2.9606 —0.050 0.118 
3.2720 —0.030 0.083 
3.6161 —0.016 0.060 
3.9964 —0.008 0.040 
$.4167 —0.008 0.025 
4.8812 —0.001 0.015 
5.3946 0.009 
5.9619 0.006 
6.5889 0.004 
7.2819 0.002 
8.0477 0.001 


modified to take account of the exchange between the 2p electrons only. The 


values of € (in cm.~') compare with Brown’s results as follows: 


FERUT 


Brown 


The author's preliminary EDSAC I results (1952) were used by Seaton 
(1951) to calculate photoionization cross sections. The present results have 
been circulated in advance of publication for calculation of various other 
physical properties for comparison with experiment. 
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As an independent check on the program, test runs were made on Fe*" 
and Fet!®, Accuracy similar to that stated above for neon was achieved, and 
the results obtained on EDSAC I were confirmed to this extent (Froese 1957a). 
The program might be capable of producing greater numerical accuracy by 
repeated application, but it was not deemed worth while to try to produce 
results unjustified by the physical approximation. A typescript ‘Manual of 
operation of the S.C.F. program’”’ (75 pp.) has been prepared as a guide to 
the use of the programs. 
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A SURVEY OF THE MATHEMATICS AVAILABLE 
FOR DESCRIBING FRACTURE! 


A. E. SCHEIDEGGER 


ABSTRACT 


An attempt is made to extend the mathematical finite strain theory of general 
rheological bodies to the inclusion of discontinuous displacements, i.e. to the 
description of fractures and related phenomena. A general analytical representa- 
tion of discontinuities is given. A particularly simple type of discontinuity has 
been singled out for special investigation, i.e. ‘‘dislocations”’. The notion of 
“dislocation”’ is investigated from the standpoint of finite strain theory and 
several theorems are proved. Finally, the requirements for the establishment of 
a dynamic theory of fractures are outlined. It is shown that the mathematical 
theory cannot be carried further without additional physical investigations. The 
direction in which these have to be sought is indicated. 


1. INTRODUCTION 

The fracture of materials is a very common occurrence. Although humans 
have been breaking things since the inception of civilization, it is an unfortu- 
nate fact that the whole subject of fracture is very incompletely understood. 
There are quite a number of rule-of-thumb criteria of fracture-——or better, of 
when. a structure is supposed to be safe so as not to fracture—but the basic 
problem of describing the progress of a fracture surface in a given body under 
given external stresses and boundary conditions has not yet been solved. This, 
of course, is owing to the fact that the engineer is not particularly interested 
in exactly how a certain structure fails as long as he can develop criteria for 
making it safe. However, in other sciences, the fracture process itself is of 
utmost interest. Thus, the present physiography of the Earth’s crust is the 
outcome of a giant deformation process in which fracture played a significant 
role. Similarly, earthquakes are certain to be the expression of a discontinuous 
displacement process within the Earth. If these subjects are to be more 
completely understood, it seems first necessary to obtain a thorough under- 
standing of the fracture processes involved. 

It is obvious that it would be a formidable undertaking to achieve a com- 
plete analytical formulation of the fracture problem. It is therefore proposed 
to investigate here, as a first step toward the final goal, at least the analytical 
means available to describe fractures in a proper, mathematical way. To the 
writer’s knowledge, not even this has ever been done. In the course of the 
investigation, it becomes obvious that the theory can only be carried up to a 
point where it is seen that further physical investigations are necessary in 
order to further elucidate the fracture process. It is hoped that the present 
survey will indicate the direction in which these have to be undertaken. 

In studying the mathematics of fracture, it becomes at once clear that it 
is necessary to invoke the formalism of finite strain theory. This is so because 

‘Manuscript received July 4, 1957. 
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the displacements, where fracture occurs, cannot be considered as infinitesimal. 
Unfortunately, this increases the mathematical difficulties tremendously. 

A particular type of discontinuous displacement, which has been called 
dislocations (Volterra 1907; Love 1927, p. 221; Burgers 1939), has been singled 
out for special attention. Such dislocations have already been applied to the 
explanation of earthquakes (Housner 1953; Vvedenskaya 1956) and it seemed 
therefore indicated to give the proper formulation thereof in finite strain 
theory. There is some hope that such dislocations may describe fractures in 
some types of rheological bodies. Unfortunately, it can be shown that dis- 
locations in general rheological bodies represent a rather restricted case of 
discontinuous displacement. In those cases where dislocations can be used to 
describe fracture, they effect a tremendous simplification of the mathematical 
formalism. 

The mathematical framework arrived at in this paper clearly shows the 
difficulties involved in any proper description of fracture phenomena. However, 
there are no reasons other than tediousness why the proper experiments could 
not be performed and the proper mathematical models could not be con- 


structed. 


2. GENERAL THEORY OF DISCONTINUITIES 
The treatment in the following pages is based upon an earlier summary of 
the writer (Scheidegger 1956), of the theory of finite displacements in rheo- 
logical bodies. We therefore again describe the displacements by the co- 
ordinates as functions of the parameters: 


(2. 1) X= x(a; t), 


where the x’s represent the co-ordinates of any ‘‘point’’ & of the body. Let 
us assume that the parameters £ are identical with the co-ordinates at time 
to as this will simplify matters. Then, if x; as a function of & is discontinuous, 
Eq. (2.1) describes a discontinuity in the body which has arisen in the time 
between ¢ and fo. 

It appears from physical experience that discontinuities are distributed in 
sheets through the material. Thus, let 


(2.2) f (é1’, fo", &3', t) = 0 


describe a surface of discontinuity S. Equation (2.2) implies that the surface 
S may change its position in the body with the passing of time. As indicated 
in Eq. (2.2), an arbitrary point on the surface S will be denoted by &,’. 

The equation of the surface of discontinuity has been given in parameter 
space, because fracture sheets are usually fixed in the medium, except for 
possible extension with time. It would theoretically be possible to substitute 
(2.1) into (2.2) in order to obtain the equation of the surface S in co-ordinate 
space. However, because of the assumed discontinuity, (2.2) expressed in 
co-ordinate space would have multivalued solutions. It thus turns out that 
the parameter space is much more adapted to describe discontinuities than 
the co-ordinate space. In the writer’s earlier paper, the formulas to transform 
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all dynamical quantities into parameter space have been put together. These 
formulas have to be used if certain quantities, in a particular problem, are 
not given in terms of the parameters. 

The rim of a surface of discontinuity is a line that must obviously be closed 
in itself or begin and end on the boundary of the body. The discontinuity at 
the surface S can be described by giving the vector of ‘“‘jump’’ X; in x; if one 
passes from one side of the surface to the other. The discontinuity is thus 
given as a vector field, defined on the surface S: 


. 9) r , 
(2.3) X, = Xi(ka’, t). 

The vector of jump can be calculated if a circuit around the edge of the 
discontinuity is drawn, as follows: 


. 
OX; 
2.4 X,= 9 — 
(2.4) =o 
ba ta 
encireling rim 


LE a. 


For physical reasons, only such displacement fields are allowed where this 
integral does not depend on the circuit (but it will naturally depend on é,’). 

Turning to the physical significance of the discontinuities discussed above, 
one notices that surfaces of discontinuity can describe, physically, various 
things. First of all, there is the possibility that they represent fracture sheets. 
The jump in displacement as one passes from one side of the surface of dis- 
continuity to the other may represent the picture of a physical disruption 
of the material, i.e. of a fracture. 

This is, however, not the only possibility. The discontinuity could be such 
that the material would have to penetrate itself in order to make the dis- 
placement possible. Physically, this corresponds to a case where a slice has 
been cut out of the body and the ends, then, have been pasted together. 

An important special case is obtained if the vector of jump X; is everywhere 
so behaved that its end point again lies on the surface of discontinuity; i.e. if 
the discontinuous displacement takes place entirely within that surface. This 
would correspond, physically, to a shear fracture. 


3. DISLOCATIONS 

A particular type of discontinuous displacement is called a ‘‘dislocation”’. 
Dislocations were encountered in the infinitesimal theory of elasticity as it 
was noted that Airy’s stress function in multiply-connected regions could 
have many-valued solutions if one tried to calculate the displacements corre- 
sponding to a given continuous strain field. It has been shown later more 
generally that the displacements in the infinitesimal strain theory of elasticity 
of multiply-connected regions are, quite apart from the method using Airy’s 
function, not necessarily one-valued although the strain components are. The 
proof, given by Love (1927, p. 222), was originally due to Cesaro. In order 
to restore one-valuedness of the displacement, one has to introduce surfaces 
of discontinuity (with regard to the displacement). Such discontinuities are 
called dislocations. 
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The surfaces of discontinuity are chosen in such a fashion that the body 
becomes singly-connected through them. It is thus evident that dislocations 
require the existence of non-evanescible circuits. In the limiting case, the 
multiple-connection of the body can be achieved by the assumption of singular 
lines upon which the strains are not continuous. These singular lines must 
either be closed in themselves or else begin and end on the external surface 
of the body. They are the rims of the surfaces of discontinuity. 

In order to prove the fundamental theorems about dislocations one usually 
uses the formulas of infinitesimal elasticity theory. An inspection of Cesaro’s 
proof, however, shows that it is based upon the compatibility relations only. 
Such compatibility relations obtain in various kinds of rheological theories 
and one would therefore expect that the notion of “dislocation’’ can be 
generalized for the finite strain theory of any type of rheological body. The 
particular assumptions made in elasticity theory may not be vital for the 
existence of dislocations. 

Thus, let us again consider a rheological body whose displacement is de- 
scribed by means of parameters and co-ordinates: 


(3.1) Xi = Xi(fa, t). 
The finite strain €g3 is defined as follows: 
(3.2) €ap = (Kap — Cas) 

where 
_ _ 9X; OXY 
a) oa 


and fag is the tensor of zero strain, which is a given and known tensor field 
(cf. Scheidegger 1956). 

Generalizing now from the elasticity theory of dislocations, we shall try 
to calculate the ‘‘displacement”’, i.e. the co-ordinates x, if the strains €,3 are 
given. We shall assume that the body is multiply-connected so that there 
exist non-evanescible circuits. In that case, we would expect “dislocations” to 
be possible. As mentioned earlier, the “thollows’’ causing the body to be 
multiply-connected can be contracted into singular lines upon which €gg is no 
longer regular; the lines, then, become the rims of surfaces of discontinuity 
through which the body can again be made singly-connected. Dislocations 
found in the above manner are, by their characterization, physically allowable 
discontinuities of a rheological body, belonging to the class discussed in 
Section 2. 

The condition that the €,3 be given continuous functions (except on the 
singular lines mentioned above) of the é’s also implies the same for kag. It is 
more convenient to base the following calculations on kag rather than on €gg. 
The problem, thus, is to calculate x; from a given x-field. The differential 
equation by which this has to be accomplished is (3.3). 

Let us assume that we have found the co-ordinates at a point A = (é1, &, &3). 
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Then, the co-ordinates at a point B = (&,’’, &2’’, &;’’) are given by the following 
expressions containing a line integral: 


Bax, , 

(3.4) Xi (Ex") = wilt) +f a diz = wile f Cia dEa 
A O&a e~J4 

with 

(3.5) Cie = Ox,/d€, 

and 

(3.6) Kap = Bpell sp. 


The basic problem in the present context is to decide what the conditions 
are that cause Eq. (3.4) to give a unique value for x;,(&’’). Obviously, this 
expression contains two possible causes for multivaluedness: the path of 
integration and the expression of dja in terms of kx. 

Let us consider first the question of the path of integration. Thus, let us 
assume that there exists a tensor field aj.(£) which is single-valued and satisfies 
Eq. (3.6). Furthermore, let us consider two paths of integration and any 
surface bounded by those two paths. The condition that the line integral in 
Eq. (3.4) is the same along the two paths is that 


(3.7) Aig = Via,B 


must be satisfied everywhere on the surface bounded by the two paths. Here 
we have set for abbreviation: 


(3.8) Qa = 0a/déa. 


We can therefore state the result of the previous argument as follows: The 
displacement cannot depend on the path of integration if x is so behaved 
that there exist single-valued solutions ajq of Eq. (3.6) which satisfy Eq. 
(3.7). It must be assumed that not al] imaginable «-fields will be so behaved 
as to permit this, but only those that satisfy a certain “‘regularity-condition” 
which may be represented symbolically as follows: 


(3.9) R(x) = 0. 


For the further discussion it does not matter at all what the regularity 
condition (3.9) is, as long as we restrict ourselves to x-fields that satisfy 
it.* Under such conditions, if the body is singly-connected so that all paths 


*An idea of the form of the condition represented symbolically by Eq. (3.9) can be gained 
as follows. If the condition is satisfied, then the transformation x = x(£) represents an actual 
displacement of a continuous medium. This can be described in such a fashion that, at the 
time ¢ = 0, € is assumed to represent the Cartesian co-ordinates of the medium (assuming a 
certain origin and three orthogonal fundamental directions). However, at time ¢, the tensor 
kag (3.3) can be regarded as representing a certain metric in a fictitious ‘strain space’’ (cf. 
Scheidegger 1956). This metric has been transformed, in the time 0 to ¢t, from the above- 
mentioned Cartesian metric dag into the metric kag. The necessary and sufficient conditions 
that such a transformation x = x(£) exists are (@) kag = kga and is analytically regular (or 
else at least twice differentiable), (6) cagdtadég is a positive definite form, and (c) the Riemann 
Christoffel curvature tensor Rixim referring to « as a metric is everywhere zero. The regularity 
condition (3.9) is equivalent to these three conditions. 
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can be continuously deformed into each other, then the co-ordinates x; are 
univalent functions if the strains are everywhere regular; and there cannot 
be any “‘dislocations’’. 

On the other hand, dislocations were defined as existing if the displacement 
turns out to be multivalued. It is obvious that this can only happen (if the 
strains are regular) if the body is multiply-connected, for then and only then 
cannot all paths be continuously deformed into each other. In other words, 
dislocations can be present only if there exist non-evanescible (also called 
non-reconcilable) circuits in the body. If there are no hollows in the body, 
then for dislocations to be possible, there must exist singular lines upon which 
the x-field is singular so as not to satisfy the regularity condition (3.9). It is 
also immediately obvious that these singular lines must be either closed in 
themselves or begin and end on the boundaries of the body; for if this were 
not the case, non-evanescible circuits would obviously not be possible. Fur- 
thermore, it is also obvious that the line integral in (3.4) taken over a non- 
evanescible circuit around a particular singular line cannot depend on the 
circuit by which one returns to the same point. This defines the strength of a 
dislocation at that point. 

In an actual body, displacements cannot of course be multivalued. In order 
to restore univaluedness, one has to introduce surfaces of discontinuity, 
bounded by the singular lines of the dislocations. For a particular dislocation, 
the surface of discontinuity can be any surface bounded by the singular 
lines. In a concrete case, of course, its position will be given by the physical 
conditions. 

The results found above may therefore be summarized as follows: 

(a) The defining characterization of a dislocation is the condition (i) that 
the finite strains are everywhere continuous and differentiable except on a 
singular line, and (ii) that the line integral in (3.4) around the singular line 
is different from zero. 

(6) The singular line engendering a dislocation is either closed in itself or 
else begins and ends on the boundaries of the body. 

(c) The line integral in (3.4) along any non-evanescible circuit around the 
singular line does not depend on the circuit. 

(d) Single-valuedness of the displacement can be obtained if a surface is 
drawn, completely bounded by the singular line (or possibly extending to 
infinity or the boundary of the body), upon which the displacement ‘‘jumps”’. 
The displacement is continuous everywhere else. 

(e) Dislocations form a special class of physically allowable discontinuities 
of a rheological body as they were generally discussed in Section 2. 

Turning now to the physical significance of dislocations, we note that, 
theoretically, discontinuities of displacement can be entirely arbitrary. It 
would be advantageous if one could confine oneself to dislocations which are 
a particularly simple type of discontinuity. In fact, it is customary to do 
this in infinitesimal elasticity theory. The justification is that, across a dis- 
continuity, the stress must be continuous (at least if inertia forces are neg- 
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lected). Since, in elasticity theory, the strains depend in a unique fashion 
on the stresses, this implies continuity of the strains also. Thus, the x-field is 
continuous, a condition characteristic for dislocations. 

In the finite strain theory of arbitrary rheological bodies, conditions are not 
so simple. Continuity of the stresses no longer implies continuity of the strains, 
not even in equilibrium problems. It is therefore undoubtedly true that dis- 
continuities of the general type discussed in Section 2 are possible. 


4. DYNAMICS OF DISCONTINUITIES 

After having discussed the geometrical description of discontinuities, it 
seems proper to investigate what can be said about the dynamical conditions 
necessary for their physical formation. 

Unfortunately, statements about the conditions for the formation of dis- 
continuities have to remain very vague. This is in line with the fact that not 
too much is known about the physics of formation of fractures in materials 
in general. 

Thus, let us assume that at time ¢; there is already a fracture surface present 
in the body 


(4.1) f(t’) = 0 
with- boundary (expressed in parameter space) 
(4.2) Fa = lZa(o), 


where o is the arc of the boundary (in parameter space). The jump vector 
at time f; is 


(4.3) X;, = 1X,(é’). 


What one would like to know is the boundaries and jump vectors at a later 
time. What can be the conditions to determine this? 

It has been mentioned that the fracture surface has a certain boundary. It 
will be necessary for physical reasons to assume that a fracture starts at a 
certain point &’ in parameter space and that the boundary grows from that 
point and sweeps through a certain region of the body, thereby engendering 
the fracture surface. At the same time, the jump vectors may change on all 
points of the surface. That part of the fracture surface which has been en- 
gendered would not change its position in parameter space. A particularly 
important case will be the one where there is slippage along the fracture 
shear” fracture where no hollows 


sc 


surface only, i.e. where there is complete 
are created inside the body. The condition for this is that the increment of 
the jump vector during an infinitesimal time dt is tangential to the fracture 
surface. 

The problem is thus to find the following two functions of the parameters 
for the time ¢: 
(4.4) Za = Za(o, ?), 


(4.5) X; = 'X4(Eq’, t). 
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We have already pointed out earlier that, in general finite strain rheology, 
the problem of finding the displacements as functions of time and parameters 
is determined by three types of conditions: (i) the law of motion, (ii) the law 
of continuity, and (iii) the rheological condition. In the last, the whole strain 
history of the body may be involved. With the presence of discontinuities, 
the conditions at the latter must be added. There is thus (i) a condition for 
the stresses across the surface of discontinuity (e.g. whether there is friction 
present) from which, by means of the rheological equation, a condition for 
the strains can be deduced (which again may depend on the strain-history), 
(ii) a continuity condition for the displacements across the surface of dis- 
continuity (e.g. stating whether the creation of cavities is allowed), and 
finally (iii) a condition for the spread of the rim. The latter again may depend 
on the whole strain history of the body. 

Experiments available to date do not permit one to present even one case 
of a complete solution of the above problem. However, the above enumeration 
of the conditions required might well suggest the direction in which further 
efforts might be expended. Perhaps it will be found possible, for instance, to 
connect the spread of a shear crack with the stresses present near the edge; 
or, possibly, the spread of a tension crack may depend only on the type of 
singularity of stress present near the edge. Such relationships could be found 
either by making the appropriate experiments or, possibly, by studying the 
physics of molecular cohesion theoretically. 


5. CONCLUSION 

If we make a survey of the mathematical theory of discontinuous displace- 
ments as far as it can be set up to date, it is inevitable that we become aware 
of its incompleteness. There is no difficulty in formulating the geometry of 
the problem and of devising a general scheme of describing the dynamics, 
although this, too, does not seem to have been done heretofore. However, 
when it comes to formulating the ‘‘law of spreading” of a crack through a 
body, we note that the understanding of the underlying physical principles 
is missing. One can state the types of physical laws that are required for the 
mathematical theory, and this might be helpful for the further elucidation of 
the problem. The missing laws are the dynamical conditions for the spread, 
speed, direction of the edge of a crack as a function of the conditions prevailing 
in its neighborhood. 

It might be argued that there is a large volume of literature on fracture 
and strength of materials (ably summarized, for example, by Orowan (1949)), 
but an inspection thereof shows that there are no instances, at least to the 
writer's knowledge, where there would be any such precise statements as are 
required for a mathematically complete formulation of a particular fracture 
problem. The closest to such a formulation may well be the Griffith (1920) 
theory of the spread of small tension cracks, in which the stresses at the tip 
of a crack have been calculated and a criterion has been formulated as to 
when the crack will spread. However, no information can be obtained as to 
the direction and speed of spreading of the crack during the next infinitesimal 
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time-interval. Perhaps it would be possible to extend the Griffith theory to 
provide this information for the case under consideration. 

Another attack to obtain the required basic physical information might 
well be by means of experimentation. If one were to combine modern tech- 
niques of high-speed photography (see, e.g., Jones 1952) with those of experi- 
mental stress analysis (see, e.g., Hetenyi 1950), it would appear expectable 
that the dynamic conditions near the edge of a crack could be established 
for various typical cases. 

Once the conditions at the edge of a crack can be formulated analytically, 
the problem of determining the spread of the crack through a body for given 
external conditions becomes a purely mathematical one. Even then, however, 
the solution of this problem is a very formidable undertaking. The problem 
is characterized by the fact that it contains a set of boundary conditions 
referring to a boundary which itself has to be determined by the solution. Such 
conditions are known as ‘‘floating boundary conditions’; they are notorious 
for the difficulties which they present if they occur in a problem. Nevertheless, 
it ought to be expected that solutions could be formed for a variety of typical 


cases. 
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ON THE SPIN-ORBIT INTERACTION IN 
DIATOMIC MOLECULES. I 


I. Kovacs 


ABSTRACT 


The purpose of the present paper is to give further details of an investigation 
of the spin-orbit interaction in diatomic molecules. The first part of the paper 
deals with perturbations between states of odd multiplicity in a two-electron 
system. With the aid of the matrix elements which have been obtained, a com- 
parison of the perturbations in the various components of a given state and in 
some cases a comparison of perturbations in two different states is possible. 
Also for certain cases it has been possible to establish relationships between 
the matrix elements which give the magnitude of the perturbations and experi- 
mentally measurable molecular constants of the perturbed states. 


INTRODUCTION 


The majority of the perturbations observable in the spectra of diatomic 
molecules—viz. those occurring between states of the same multiplicity — 
can be easily understood by taking into account the terms omitted at the 
separation of the wave equation with the aid of the perturbation method. 
As shown by many examples, however, perturbations can also occur between 
states of different multiplicity (Schmid and Ger6é 1935; Kovacs and Lagerqvist 
1950; Lagerqvist, Nilheden, and Barrow 1952; Nilheden 1955) or, in other 
words, anomalous phenomena have been observed which can be explained 
only by the occurrence of such ‘‘intercombination” perturbations. The anoma- 
lous behavior of the ‘II state of the O.+ molecule serves as an example for the 
latter (Bud6é and Kovacs 1954). An exact discussion of the intercombination 
perturbation could be carried out if an exact relativistic form of the wave 
equation were known for a many-electron system. Lacking this, one can 
start from the form of the Dirac equation as extended to an atom with two 
electrons, and the occurrence of terms in the Hamiltonian operator which 
are due to the magnetic interaction between the orbital and spin momenta 
of the individual electrons can be assumed. Taking these terms into con- 
sideration, attempts have previously been made to explain the intercombina- 
tion perturbations (Budé and Kovacs 1938, 1939, 1950; Premaswarup 1954a, 
b). These investigations, however, with the exception of one (Budé and 
Kovacs 1954), were of a fairly general character and thus were not adequate 
to draw conclusions concerning the magnitude of perturbing matrix elements 
for various electron configurations. The purpose of the present paper is to 
take into account individual electron configurations of two- and three-electron 
systems and to determine the values of matrix elements arising from the spin— 
orbit interaction. These matrix elements of various configurations are com- 
pared and related to the molecular constants known from the experimental 


data. 
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In order to attain the aims set forth above, the spin-orbit operator, as 
well as a detailed form of the wave-function of the individual electron con- 
figurations, will be needed. The former has, as is well known, the following 
form: 


(1) > a,(1,8;) = ray 4s ss ee (retin? Nee iSy;) 


-@ "(sie 9) bid) | , 


where ¢; is the azimuth angle of the 7th electron around the z axis which lies 
along the axis of the molecule and s;,, sy;, sz, are the well-known Pauli opera- 
tors. The operators f;, gi, and a; are independent of ¢; and their form depends 
on the nature of the model applied. 

To determine the electronic part of the wave-functions we shall start from 
the wave-functions of the individual electrons of the following form (Budé 
and Kovacs 1954): pie™#*ta,; and pe *iB,, where a; and @, are the wave- 
functions belonging tothe +4 and — 3 values of the spin quantum number o jy. The 
wave-functions independent of g; and corresponding to the state represented 
by the quantum numbers ”j, ip, Xip of the 7th electron are denoted by the 
symbol ~p; = Pn,,,4;,d;p- Similarly, according to the quantum numbers 1,,, 
leg, gg the wave-function gz = Qnig.tkg Akg Tepresents either for k ¥ 7 another 
state of another electron or for k = 7 another state of the same electron. 

Disregarding the mutual electrostatic interaction of the electrons, the 
wave-functions of the individual configurations can be given with the aid 
of the above-mentioned wave-functions in the form of determinants. These 
wave-functions of determinant form are in general the simultaneous eigen- 
functions of the operators L, = >}°;/., and S, = );s.,, whilst they are not 
eigenfunctions of the operators S? = S,?+5S,?+5S,°. It is, however, easy to 
form linear combinations of the above-mentioned eigenfunctions of determinant 
form which will be eigenfunctions of the operator S? with eigenvalues S(S+ 1) 
determining the multiplicity of the respective terms. 


TWO-ELECTRON SYSTEMS 

1. Eigenfunctions 

Assuming that all the electrons of a molecule, with the exception of two, 
are in closed shells, the states of the molecule might be characterized by the 
quantum numbers of these two valence electrons in the following manner: 
(1 ips Lips Nips F ip Nar legs Negr Fxq) OF, Since in Many cases the writing in full of the 
principal and azimuthal quantum numbers is unnecessary, simply (A,, A,). In 
this case the state (jy, lip, Nip, Tip) has to be understood as a state in which 
the ith electron is in a state denoted by p, where i = 1 or 2, whilst (mxq, dig, 
Arg» Teg) indicates that the kth electron is in a state denoted by g, where k = 2 
or 1. The possible types of states of the molecule are given by the different 
possible positions of the \’s relative to one another. Further on, the simplest 
cases will be dealt with, where |A,| = |A,| = A, or |A,| = A’ and [A,| = A’+1, 
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or |A,| = A’+1 and |A,| = X’. In the case of (A, A) it is well known that the 
following types of states are possible (Herzberg 1951): 'Z+, *+, 1(2d), #(2d), 
where for \ = 0 the terms =~ are lacking, whilst in the case of (X’, \’+1) 
and (\’+1, \’) the possible states are: 'II, *IT, !(2\’+1), #(24’+1). If \ or X’ 
= 0, the terms (2A) and (2A’+1) are not present, whilst if \ or \’ = 1,2,... 
they signify A, [...or ®,H... etc. types of states. 

Further on, the cases occurring most frequently, viz. interaction between 
> and II states, will be investigated. The electron part of the needed wave- 
functions assuming Hund’s case (a) is, on the basis of the formulation men- 
tioned above, the following: 


'e('S3) = 6('S*)(1/V2) (aiB:— Bia), 
(2) | o(*?;) = 6(°2*)a1a», 

| o(°*s2) = o(*=*) (1 /+/2)(a1B2+ Bia), 

(241) = o(°E*) 6182, 


where 


| ¢('=*) = (1 /2) (Pigo+ f2q1) ee Maia) 
(3) sin 


| Sw: j , COS 
| @(2*) = (1/24) (pige¥F fogs) sin \(g1— ¢2) 
for \ + 0. If \ = 0, the normalization factor will be 1/2/2.z instead of 1/2zx 
and as a matter of course owing to the sine factor only the wave-functions 
corresponding to the upper signs differ from zero. 

Similarly for II states we obtain: 


OC M41) = (11) (1/+/2) (ai82—Biaz), oC I’_1) = $('M’) (1/2) 

X (Bia2—a182), 
(4) | oCTy2) = 6M, $(°T’_2) = $(°T’)BiB2, 

CM) = $1) (1/2) (a182+Biae), o(T'1) = $(*T’) (1/2) 

X (Bia2+aiB2), 


Il 


o(*Mo) = o(°I1)B1B>, (1's) = (°M' ares, 
where 
( (etna) = (1/2/22) pee OO Og OOO), 
9) +4[d' (g1—¢2)—g2! 


i[\’ (¢1 v2)+91)) 


6 IW) = (1/22. 23) (pigne & poque 


for configuration (X’, \’+1) and 


Il 


(o(*'M) (1 24/2. 3) (Pigoe wortoil + pags 1X! (1—92)—#2] ) 


(6) a , é 9 + ’ —~2 2 
le? 1’) te (1 24/2.2) (pigze i[\’ (1 votes pa (¢1 2)—v2] ) 


for configuration (A’+1, \’). 

Taking into account that the energies corresponding to the two possible 
positions of the resultant 2 relative to the axis of the molecule (+2 and —2) 
are degenerate and completing the above-mentioned wave-functions by the 
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part arising from the rotations and vibrations of nuclei, the new wave-functions 
will be the following: 


re: >i = CER... oe 0s 
jy(z ex) = (1 V2) [oC S2))u, sit OC 271) z,-11Ra os 
y(°24)'= /V2)[oC ze us 41- OC EF) us, Ra, os 


(7) VC "e5) = CZa) RK. py Uy,0; 
mA i, r V2) [o(' IT,1) . )uy,- JR, vy 
| wCme) = (1/2) [o( Mig) ty,40+ (To) 7,-0]Rn. v» 


where R,,.,(p), wy.g, and v are the vibrational and rotational wave-functions 
and the vibrational quantum number of the molecule, respectively, and 1 
designates all the principal quantum numbers of the electrons. 

Since for Y states the approximation obtained from Hund’s case (a) has 
only a theoretical meaning, whilst for II states case (0) is a limiting case as 
important as case (a), it seems advisable to use the wave-functions of case 
(6) for the states considered above. These wave-functions are given as the 
linear combinations of the earlier ones where the coefficients of the linear 
combinations can be easily determined by taking into account, by the pertur- 
bation method, the terms omitted at the separation of the wave-equation. 


Hence 
y(‘z4) = y(‘=4), 
oa. i 748 an a : ; Sey l bn 
IWC2F41) = “W 3741" zi) -V/ a741 ¥ Si) = — lv CZ) 
| +¥(S4)], 
l¥(25) = ¥(25)’ 
= 
3 / a - J+1 a | 
yt za . Syt Ree by 
ieee i) = 4 2741 ¥O20-4 a7+i ¥6 =) = Alvi >i) 
| —y (24), 
y's) = y(‘m4), 
‘eer x 7 FEF -<8) J(J+2)_ 2b 
gs) (YOM = VW sernU+H WUD+W BIH as MOND 
l 357+ l Bes 1 ko 
+h sare Is) = 5¥( )+ a4 Ii) +54 C1), 
| / 
er re (J+2 i l os 
pesanee V “3+ ¥ CUE) — V FoEH VOD) 
| Il 2 3 i l r 
+—s¥(IIs) = —- = 1S) +54 (714), 
Ve Ve Ve 


| / / . 
ia ee a + / (J-1)(J+1) —_s 
Cera) ay OTH yy HCI) Qs+is7 Vy 


/ i Deg I ee, oa Ss 
YW 3 ant Ts) = 5¥ CMS) —~, 0 CS) +50), 


2 - 
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where the wave-functions of the left-hand side are characterized by the 
quantum number K as is usual in case (0). In those expressions, however, 
where the matrix elements between the terms of case (a) appear independent 
of the rotational quantum number J, the wave-functions of (8) should be 
taken in the simpler approximate form on the right-hand side. 

The perturbation matrix elements appearing between individual terms can 
be easily computed with the aid of the wave-functions obtained above and the 
operator (1). 


2. Perturbations AA = 0 

Primarily those cases will be dealt with where the terms perturbing each 
other belong to the same value of the resultant orbital angular momentum, 
viz. AA = 0. In such perturbations only the first member of the operator 
(1) vields non-vanishing matrix elements. 

(a) Perturbation }\2—*> 

It can easily be verified that between terms of the same symmetry in a 
first approximation there is no interaction (see Fig. 1), thus 


(9)° BOE Sh) = CTS) = BCS SS: = 6: 
E ay? 
~~ 





Fic. 1. 


Between terms of different symmetry we find the following matrix elements 


(see Fig. 2): 


ao) JHCE5a ED) = -HC2d4'53) = FW/2V2)\ay tay), 
H(*s* st) = 0 


The magnitude of the constants a, and a, will be dealt with later in the dis- 
1S (or 43%) 


a 


cussion of the selection rules. As is seen from (10), due to a 
state, perturbations of the same magnitude, and independent of J, will appear 
on the K = J+1 and K = J—1 components of the *2+ (or *2~) state. The 
order of magnitude of the perturbation might be corresponding to the values 
of the matrix elements appearing for heterogeneous perturbations of the same 
multiplicity. 


*The perturbations discussed in this paper are everywhere of the first order of magnitude. 
Therefore in some cases it might happen that here between certain components no inter- 
action exists in spite of results cited (Budé and Kovacs 1938, 1939, 1950; Premaswarup 1954a, 

These interactions are, however, in every case only second-order effects. 
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E ast 
oa—_—— 





Fic. 2. 


(b) Perturbation 1—*II 

Owing to the fact that in several cases for a *II state, at smaller rotational 
quantum numbers, Hund’s case (a) is a better approximation whilst for higher 
quantum numbers Hund’s case (0) is better, here both limiting cases are 
considered. As a matter of fact the best approximation would be obtained by 
an intermediate case; however, the formulae are fairly confused in this region. 

(1) '1—11 (a). 


(11) HCE M4) = —3[0('+1)a,+Xc,], HT m4) = HCH, *4) = 0 


(see Fig. 3). 





(2) I—I1(b).— 


fans *nt,,) = —H (nt Wt.) = —(1/22)[(\’+1)a,+)'apl, 
H('1i4 *1+) = 0 


(see Fig. 4). 

(c) Perturbation 1}>—'d* 

In case of !Y terms of the same symmetry (thus 'S+—!Z+*), operator (1) 
does not supply interaction matrix elements. The electrostatic interaction 
and certain terms omitted at the separation of the wave-equation do, how- 
ever, give matrix elements independent of J, which sometimes might be very 
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Fic. 4. 


large (Ger6 and Schmid 1943). Between terms of different symmetry, thus 
between '>+—!Z**, no interaction exists according to the aforesaid operators. 

(d) Perturbation *>—*X* 

Between components of *2 states of the same symmetry belonging to the 
same A the above-mentioned interactions similarly yield matrix elements 
independent of J, whilst owing to the operator (1) the perturbation between 
AK = 0, +2 components can be regarded as a second-order effect. 

On the other hand between *> states of different symmetry (thus *2+— 
8+) the operator (1) gives matrix elements of the first order of magnitude: 


SHC2541°E3) = HCZS S541) = HCZS*Z51) 
= H(*d4_,*S4) = +(i/2V2)aA(a,—a,) 


(see Fig. 5), whilst the other matrix elements are zero. 


(13) 


a* 





E 





Fic. 5. 


(e) Perturbation 'I1—'II* 

Operator (1) does not give matrix elements in this case, the perturbation 
being due only to the above-mentioned electrostatic interaction and to the 
terms omitted in the separation of the wave-equation. The matrix element 
is similarly independent of the rotational quantum number J. 

(f) Perturbation *1I—*II* 

Operator (1) and the above-mentioned operators give perturbations accord- 
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ing to the selection rule AQ = 0 if both terms *II belong to Hund’s case (a), 
or according to AK = 0, +1 if both terms belong to case (6). 


3. Selection Rules and the Values of the Constants a, and a, for Perturbations 
AA =0 

The perturbations dealt with previously do not occur in each possible case, 
but only if certain selection rules are satisfied. Such a selection rule, according 
to which operator (1) yields interactions only between »¥ states of different 
symmetry, has already been established. Other rules may be found by examin- 
ing the expressions of the a, and a, constants occurring in each perturbation. 
Thus 


ro a ok 
| Ap = | R, r(p)) Prip. lip Nip@iPn’ ip Vip r,,d7 (f Inky tkq: Ang dn’ ka U'kq Neg@Tx) 
XK Ry. (p)dp, 
il 


sige +k 
| Q, = J R,, »(p)f Inka. leq Agg@kQn' kq Ukq reg te (J Pnip. lip. dipPn’ ip. Vip deeeTa) 
X Rn,» (p)dp. 


(14) 


Those matrix elements from which these constants are derived, owing to 
the integrations carried out over the spin coordinates, are diagonal in o,, 
and o,,, thus Ao;, = Ao;,, = 0; also these expressions, owing to the integration 
already carried out over ¢; and ¢g,, are diagonal in \;, and A,;,, hence Ad;, 
= Ah,, = 0. This means that even if every other selection rule were satisfied, 
perturbations can take place only between Y terms belonging to the same 
electron configuration of the same \; e.g. if both Y states belong to the con- 
figurations (7, 7) or (6,6) and so on. (Assuming a configuration (¢, ¢) only 
‘y+ and *Z* states will exist and according to the selection rules already 
mentioned no perturbation occurs between these.) In case of I] states, per- 
turbations occur only if both II terms belong to the configuration (¢, 7) or 
(7, 6) or (7, ¢) and so on. 

Thus conclusions may be drawn concerning the electron configuration of 
the perturbing term from the occurrence or absence of the perturbations. 

Taking into account some assumptions dealt with in detail in Section 5, 
further separations might be carried out and accordingly it can be readily seen 
that a, and a, have to be diagonal in /;, and /,,. 

Finally owing to the orthogonality of the functions g and p, a, has to be 
diagonal in ,, and a, has to be diagonal even in 1;,. This statement vields 
the result that the values of the bracketed integral in the expressions of a, 
and a,, owing to the normalization of the wave-functions, are unity. 

According to the assumptions described in Section 5, the functions p and 
g as well as the coupling constant a do not depend on the internuclear distance 
p; hence (14) might be written simply in the following form: 


fay = (my, 1" p)IRy. (0) Ru, (p)dp, 


(15) - 
la, = a(Meg, 2x) Re. o(p)Rn’. 2 (p)dp, 
where 
‘ * 
ty a(Nip, 1" ip) = a ae 
(15a) ; > x 
a (Nx, n ka) = J Inka. Ieq. Akq@kQn' ka. pe Dag Fi 


and: 4= 1,2..k =—-2; 1. 
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If An,, ¥ 0, then according to (14) a, = Oand if at the same time An,, # 0, 
then also a, = 0. This means that in such cases no interaction at all exists 
between the two states since both constants vanish. 

If one of the requirements An;, = 0 and An,, = 0 is satisfied one of the 
two constants differs from zero. In such cases the value of the remaining 
constant can only with difficulty be related to the multiplet splitting constants, 
which can be determined experimentally. Fortunately in general these states 
are energetically far from one another; therefore a perturbation with inter- 
section between them will seldom occur. 

For An; = 0 and at the same time An,, = 0, viz. An = 0, the functions 
Pripstip dips Pn'ip tip dips AND neg, tea dq Yn'kq.U’kqdkqgr Which occur before and after 
a, and a, respectively in the wave-function of the two terms perturbing each 
other in equation (14), will also be the same. Under certain circumstances 
these constants can be related to the multiplet splitting constants 1, of the 
‘II term originating from the same electron configuration. 

In case of equivalent electrons all of the quantum numbers 7 jy, / ip, \ip agree 
with the values mx, Jig, Agg and the two electrons might differ only in the 
orientations of the \’s or o’s. In such cases evidently Ppy,,.1;, sp = Qnia.tig dia 
and Pniy.tep Akp = YUnkg. tka Akg aNd thus a, = a, = a. From the wave-functions 
(2) and (3) it follows at once that for a configuration \* with \ = 0 only a 
term 'X+ can exist, whilst with \ # 0 only the terms !Z+, *Y~, and 1(2A) can 
exist. In order to normalize the wave-functions occurring in equations (3) 
in the case of equivalent electrons they should be divided by 2. According 
to the present selection rules an interaction can only exist between !2+ and 
’Y~ states belonging to the same configuration \? 


~ 


(16) H (3541 '25) = —H(°25-1'25) = +(i/V2)da. 
If Agg = Awl, but nip = meq and |i, = ky, then the above-mentioned 


equation will hold also in (5) and (6) so that in this case also a, = a, = a. 


| H(t *m) = 4(2\+1)a 
(17) 2\+1 


2/2 4 


Summing up the foregoing briefly, for AA = 0 in Hund’s case (a) the follow- 


H(1F *14,:) = —H(‘15 T_,) = 


ing holds: 


(18) HH (nip, bss Nips o ip; Nkas bees Arg Orgs N,v, aA, 2); 
(1° ips Dips W’ tps O tps Beas Vege eqs Oki MV’, A’, Z')] HO 

for An, = 0 or An, = 0, Alip = Algg = 0, AXtp = Adee = 0, Avin = Ace, = O, 
AA = AZ = 0. 

According to what has been said, in Hund’s case (6) instead of the selection 
rule AS = 0 the rule AK = +1 will hold. 
}. Perturbations AA = +1 

For AA = +1, the part inside of the angular brackets of operator (1) 
supplies the perturbation matrix elements. Here the following cases are 


possible: 
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(a) Perturbation 12—*II 
As above, for the two limiting cases of *II states, perturbing matrix elements 
are also given here: 
(1) 's*—*II(a).— 
(19a) H('2s *Ts) = (—1)" él, 
where the exponent ¥ is equal to 0 or 1, according to whether we are dealing 


with a 'D+ or a !D~ state (see Fig. 6). The constant is expanded further in 
Section 5. Here it should be noted that the upper and lower signs on the 2 


Tia) 
c Ss 





Fic. 6. 


and II symbols are to be taken in the same order and hence no interaction 


exists between '2} and *II, states or between '> and *II* states. 


(2) '=+—"11(6) — 
st 3 441) = os I('> ‘= 76.,) = a E|, 


(19d) 


(see Fig. 7). 














KOVACS: SPIN-ORBIT INTERACTION. | 319 


(6) Perturbation *>—'Il 


i 


| dy - dy + —i) | 
a | C2541 05) = HCzS 5) = SS Ig, 
(<0) ‘ = 
H¢x$'13) = —S2 ig 


(see Fig. 8). 





Fic. 8. 


(c) Perturbation }2—'Il 

Since this perturbation has been dealt with by several authors already, this 
case is now mentioned only for the sake of completeness. Operator (1) does 
not supply matrix elements but they are supplied only by the terms omitted 
at the separation of the wave-equation. It is known that these matrix elements 
are as follows: 
(21) H (‘sé 4) = —(—i)'2nvV/2J(J4+1). 
The terms omitted at the separation of the wave-equation will supply matrix 
elements also for the case of the selection rules AQ = +1 or more explicitly 
in our case even for AA = +1 and AY = 0, which depend to a large extent 
on the rotational quantum number J. A more detailed consideration of constant 
n occurring there will be given in Section 5. 

(d) *X—II Perturbation 

Taking into account the interactions discussed above, then for a *2—*II 
perturbation, operator (1) as well as the terms omitted at the separation of the 
Wwave-equations will give perturbing matrix elements. Since in general the 
terms supplied by the two operators might be of the same order of magnitude, 
they should be taken into account simultaneously. In the earlier papers 
(Hebb 1936), operator (1) was more or less replaced by operator A(L.S). 
In the computation of terms diagonal in A and & these operators are equiva- 
lent but for non-diagonal elements (A, ©; A+1, 21); however, the two 
operators will give the same value of the matrix elements only if owing to the 
configuration neither of the two states perturbing each other will have inter- 
combination interactions with other states. A decision as to the correct 
operator can easily be made. For intercombination, perturbations cannot be 
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interpreted by the interaction A(L.S) and since their occurrence has been 
observed experimentally it is evident that the behavior of a molecule may 
be described more exactly by operator (1) than by the operator A(L.S). 
Thus results supplied by operator (1) are to be accepted as correct. 


(1) Perturbation *>—*II(a).— 


oe ge wd ae ---~—-——- 
H (2341 M2) = (—2) = vE-DOs) <a e 
| = (i) 9v/2(J—-1)(J+2), 


ae J 
(-—i) oA sc5 [2 +2n(J+2)] 


> 
~ 


| = (—1)"[2§+2n(J+2)], 

Fe ae ia 
a shea sz,/ JH 

| H(°2$41 It) = (—i84/ FE et ons+2) 


ll 


H(z%,; 1) 


~ (=) ep ons +2)), 
V2 


3 


be 
(22) | H(25 “T) 
| HCZ5 “14) 
| H(?=* *114) 


—(—i 2nr/(J—1) (J +2) , 
— (i) V2. (4E+2n), 
(—i)*29/J(J+1) , 


I 


ll 


Os) = (—I724/ Sein VDD) 
| = (nV20=) G2), 
Lnexs.0) = (—a8v2.4/ ZL eon 7-0) 

| es = (—i)°[48—-2n(J—1)], 
H(°E5-1 Tl) = — -0 fa [——2n(J—-1)] 

| = —(-i)1§-2n(J-1)1. 


II 
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As can be seen from (22), here only a single matrix-element is independent 
of the rotational quantum number J, i.e. 7(*2* *I+). Of the others 


A(*z+,, *Mt) and H(*=+ 


J~1 


3]IF) 


as well as H(*2+ *II+) and H(*=+ *I+) for the same J are approximately 
equal to one another and vary proportionately to J. (In Fig. 9, these are 
connected by a thin line.) Since one of the states, *2, belongs to case (4) whilst 
the other, *II, to case (a), no selection rule exists and at each intersection 
perturbations will appear. 


(2) Perturbation *2—*II(b).— 


H(2541 “W41) = (—4) -a4/ 5 5 a [E+4n(J+1)] 


= OP HV FD OF) 





} x= 2 J —1 P 
IC 2541 My) = (—i) ei eneat aT en R ONT EM te 
V/2(2J+1)(J+1) 
| I CZ5 41 i 1) = 0 
| “3 | =a “\z _ J(J+2) _ _(-if 
SC ee) = tk V5 QJ+1(J+1) = * 
CES MS) = (iP WV 2I +1) Ant —. 4 
wae ( IV 2I(I+1) 2n an a 
= (—i)2nV/2J(J+1) 
esis sz , /(J—-1)(J#1) -i) 
|H(C=s 3-1) = (-7) ae ~ ‘= ~§ 
| H(*2 25. 1 "Ss = 0 e 
\aCzg.3) = —(-)7—2tt_;--“ ; 
| we (2 J+J x 
| H(C=4_, 13-1) = —(-i? 4/ =— te dn J] 
= -‘SY etn G=T 
According to (23) matrix elements corresponding to AK = +2 are zero. 


Those corresponding to AK = +1 are constant and approximately equal 
to one another (in Fig. 10 these are connected by a thin line), whilst the 
elements corresponding to AK = 0 depend to a large extent on J and are the 
largest ones. 

This can be seen by noting that logically in the ideal case (a) » = 0, whilst 
in case (6) in (23) & = 0, since in (8) the coefficients of linear combinations 
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of the *II wave-functions are calculated similarly in the two limiting cases. 
On the other hand, we are of the opinion that (22) and (23), in the forms given, 
will provide a better approximation to the ideal limiting cases, in case (a) 
for a small J and in case (6) for a large J. 





Fic. 10. 


(3) The A-type splitting of the *II state in the limiting cases.—As is well known, 
the A-type splittings observed in II states are due to the perturbations by the 
more distant = states (Hebb 1936). According to what has been said in 
Sections (a) and (d) perturbations might be caused in *II states by 'Y and 
’¥ states.? The '2~ and *+ states on one side of the *II state (energetically 
above or below it) and these 'X*+ and *Z~ states lying on the other side of 
it will give rise to a shifting in the same direction. Taking into account all 
these interactions the following expression is obtained for the A-doubling: 


JH Cy >) |*—|H CMs 2.) |? 


» (3 = 22 J+1 
(24) hAv("M,) = (—1)""* De hv(I, ©.) ; 
where 7 = 2, 1,0 in Hund’s case (a) and 7 = J+1, J, J—1 in Hund’s case 


(6), whilst k is the index of the summation over each perturbing term 'Z*t, 


ty~, 82+, 82>. On the basis of Sections (a) and (d) expressions can be given 
for both Hund’s case (a) and case (6) coupling which hold rigorously not 
only for the two ideal limiting cases, but also hold well into the intermediate 
coupling range in their behavior with respect to variation with J. 

These expressions, in case (a), are the following: 


| hAv(*IIs) = 0, 
(25) ‘ hAv(*I;) = CoJ(J+1), 
hAv(?Mo) = Co+2C:+2C2 


3In case of two-electron systems *Y terms do not exist. 
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and in the neighborhood of case (8), 


Bele I J+2 
hAv(My41) = a 


l ¢ _ . 
| 500 2741 +C1(J+2)+C2(J+1)(J+2) 
(26) {hAv(Iy) = ahi —C,—C2J(J+1) 
P 1, J- ; 
hAv(*My-1) = 305 oT —Ci(J-1)+C.(J-1)J, 
where 
be cs (—1)"eC Zz, )1° (—1)° “|gC ML 
bas “as » hv (II *S, ~ hv (I'S, . 
(26a) 
io (—1)*|é ia TI *S,)| (1 
= 4D — 5 *s,) oe = 8) an 


From selection rules dealt with in Section 5, it can be seen over how many and 
what type of © terms the summation should be carried out.4 

Concerning the application of (26) we note, for instance, that it gives a 
fairly good approximation for the A-doubling of the B*II state of the N» 
molecule even for J = 30 (Hebb 1936) (Av(*II,,1:) = 0.98 cm~! instead of 
1.18 cm—!, Av(#IIz) = 0.89 cm —! instead of 0.86 cm—!, Av(*IIjz_1) = 0.65 
cm! instead of 0.48 cm~!) in spite of the fact that at this quantum number 
the splitting term hAv = C.A(K+1) (which is entirely valid in ideal case 
(6)) in (26) is much less important than the terms with coefficient Co. For 
higher rotational quantum numbers the approximation is even better. 


Selection Rules and the Values of the & and n Constants for Perturbations 
AA = +1 
As before, a closer examination of the constants £ and 7 leads to selection 
rules. The values of these constants are the following: 


_ . : 
(27a) é[(A, A)(A+1, AD] = 22 J R¥i (oe) J P* nip. pri fit rdgil 


x Pn’ ip. Uip At dt; | Poss lkq AVn’ kag. lkq dr; R, »(p)dp 


or 


II 


" 1 - 7 : 
(276) ef(\, A)(A, A+1)] 2724 R*,. (0) J Q* neg. tea XQel fe (LEA) ge] 


x Yn’ ko. Ug 2410 Ts SOPs. csi, Vip rdT Rn,» (p)dp, 


and two similar expressions for 7 except that the term B = h?/82*up* appears 

instead of a; and a,. This term does not depend on the integration coordinates 
; and 7;, but only on the internuclear distance p. 

In the matrix elements where ~ given by (27a) occurs, since the integration 

is carried out over the spin coordinates, o' jp = o ipl and o'¢g = ox, while 

where & from (276) occurs, o’ jp = ip and o'pg = ogg 1, and finally in matrix 


‘For a more accurate discussion of Co, the term arising from the spin-spin interaction 
should also be taken into consideration. 
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elements where n occurs, everywhere o’ ;, = oi) and o’., = oxg. In the integra- 
tions carried out previously over ¢; and ¢g, a closer examination of the constants 
leads to the fact that their value differs from zero only if in (27a) and the 
corresponding expression of 7, \’ ip = Nip] and Neg = Ag Whilst in (275) and 
the corresponding expression for 7, \’ ip = Nip aNd A’gqg = Aggel. Since in this 
investigation one of the terms is a ¥ term, it follows that |\;,| = |Agg] = ¥ 0. 
If \ = 0, then each constant has to be multiplied by +/2. These selection 
rules expressed by the resultant orbital and spin quantum numbers of the 
molecule give the well-known selection rules, namely, for the £ constant, 
AA = tland AY = ¥1, that is AQ = O, while for the 7 constant, AA = +1 
and A> =-0:(or AQ’ = :1). 

Assuming a special model for the molecule, simple expressions for the 
constants € and 7 can be obtained. Considering the two nuclei of the molecule 
united in the mass center (the so-called united atom approximation) let us 
take this united nucleus as the center of the coordinate system. With this 
model the following equations will hold: 


(28) fi = hO/dd,, gi = hcotd,, 

where #; is the polar angle of the 7th electron and where the orbital angular- 
momenta /,; and ], of the two electrons are being referred to the origin. The 
wave-functions of the individual electrons might be thus written: 


(29) Pay, tip hig CO v;) Ly ip. tip 11) O agp dip (Vi), 

(30) nig. tha deg Ter Ot) = Langa, tug (1) tag deg (8x) 

where the functions 9, (¢%) in which we are particularly interested are the terms 
of the derivatives of the Legendre polynomials. Using these, the expressions 
of — and 7 have the following form: 

(8la) é[(A, \)(A41, A)] = (1 2/2)[a(f+ U4 )g)]p f R* .(p) Rn’. » (p)dp, 
(316) (A, YA, ALD] = (1/2V2)[a(f+ UA) g)]¢ f RE o(p) Rv» (p)dp, 
(32a) f(A, (AHI, AD] = (1/272) [f+ 4 )g], f RF (0) B(p) Ru.» (p)dp, 
(32b) f(A, (A, ALD] = (1/2V2) [f+ 04 gl f RF »(0)B(p) Rn’. » (p)dp. 
After integrating the orthogonal and normalized functions g and p for equa- 
tions (3la) and (32a) values differing from zero will be obtained only if 


An;, = 0 and Al,, = 0, respectively, and for equations (316) and (328) only 
if An;, = 0 and Al;,, = 0. Moreover: 
(33) [a(f+(Q+A)g)], = J ae 8 Na(r [fit Ur)gi] 

X Pn'in. Vip rti(s, 8;) sin 8; dd; redr:, 


and a similar expression will also hold for [a(f+(1+))g)],. If we replace 
the values of » and q with those obtained from (30), we obtain the following 
(Condon and Shortley 1953): 


(34a) [a(f+ (1+ A)g)], = Fa(nip, 2’ yp) V(t —-AALL , 
(345) [a(f+ (ItA)g)]o = Fa(Meg, 2’ rg) V lala +l) —AALN) 
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where 
6 , r 7% 2 
(35a ) a(n ips n' ip) _ | Laity, lip (ria (ri) Ln’; >. Vip(Tai dr; ’ 
as r # 2 
(356) A(Meg, Mig) = J Larva. tag (Pe )@ (1h) Ln’ kg. tug (Ve) Ve are « 


Owing to the fact that the integrations are carried out over #, further selection 
rules occur as follows: in (35a), (84a), (31a), and (32a), quite apart from 
what has been said previously, A/;, = 0 and further in (358), (848), (818), 
and (320), Al, = 0. 

If An;, #0 and An,, #0 hold simultaneously, then according to this 
model no interaction exists, since the expressions (31) and (32) are (owing 
to the orthogonality of the functions Z(r)) zero. If An,;, 4 0 and Any, = 0, 
then expressions (310) and (3206) are zero, while if An;, = 0 and An, ¥ 0 
expressions (3la) and (32a) are zero. This means (assuming for instance that 
the S state belongs to a (z, 7) configuration), in the first case, that the terms 
in question can interact only with II states belonging to (67) and (oz) con- 
figurations whilst, in the second case only, they can interact only with II 
states belonging to (76) and (ac) configurations. 

For An;, = 0 and An,, = 0, i.e. An = 0, the expressions for — and 7 take 
the following forms: 


I —— : 
(36a) é[(A, A)(A#1, d)] = F 5 jaty Violet 1) — AOE) ; 
(360) E(A, A)(A, AEL)] = Fy jae Vise F1)— AON) . 
(37a) nl(A, A)(A41, d)] = 575m V1p(lp $1) — MAI) , 
(376) nl (A, A)(A, A#L)] = Fy jpBen Vil F1) — A(A1) , 
where 

| ay = (ip, Bip) J R* .(p) Rw.» (p)dp, 
(38) | = a(Mig, Meg) { Re o(p)Rn’,e' (p)dp, 


Ben = { Ry, .(p)B(p) Rn’. » (p)dp. 


The constants a, and a, appearing here will have the same meaning as in 
(15) if there m,, = m'y and my. = n’x,. Concerning the constant By, it can 
be said that it cannot be larger than the larger of the rotational constants of 
the ~ and TI states perturbing each other. 

If the Y states belong to equivalent electrons, then a, = a, = 0 and further- 
more /, = 1, and thus the two expressions (a) and (6) for — and 7 are every- 
where equal to one another. At the same time the constants should be multi- 
plied by +/2. This means for instance that for a = state belonging to a (z, 7) 
configuration it makes no difference whether the perturbing II state belongs 
to a (o, r) or a (x, o) configuration. 

Summing up what has been said above, for AA = +1, then in Hund’s 
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case (a) for the matrix elements given by operator (1) the following relation 
will hold (namely for those containing the constant &): 


(39), FEE WMet ies Lis Ng phous tec bess Naess Oker ey. Opa) 
, , , , , , , ’ , , ’ x7 
(n ips l ips r ip» o ips n kq bxes d kas o kas n ’ Vv, A ’ z= )] # 0 


for Anip = 0 or Any, = 0; furthermore A/;, = Al,g = 0, and at the same time 
either (1) ANip = 41 and Adz, = O or (2) AXy = O and Ady, = £1, viz. AA = 
+1 together with either (1) Aoi = 1 and Ao,, = 0 or (2) Acy = 0 and 
Aoxg = ¥F1 (i.e. AY = #¥1), then as in (18) the selection rule AQ = 0 is also 
valid here. 

In Hund’s case (b) mstead of selection rule AS = #1, AK = 0, +1 becomes 
valid. 

For matrix elements depending on the rotational quantum number J which 
arise from the terms omitted at the separation of the wave-equation in 
Hund’s case (a) (namely for those containing constant 7) the same selection 
rules are valid, with the difference that in case (a) we always have Agi = 
Aoxg = 0, i.e. AY = 0 and hence AQ = +1, whilst in Hund’s case (6) instead of 
AY = 0 the rule AK = +1 will be valid. 

To summarize the above, we present a table (Table 1) from which it can be 
readily seen which states belonging to what kind of electron configurations 
can interact with each other. 


TABLE | 
Electron Electron 
configuration Terms of molecules configuration 
oo IDt. spt oo 
oT I, 41 ro 
Tv 1 Dt, 82, “Saye te, Ay A 7 1 
7 6 I, *11, 1p, % 5x 
St, 2%, ; ’ 


6 6 Is + 3yt 1y- 35 se 37° 6 6 


The operator (1), as well as those arising from the terms omitted at the 
separation of the wave-equation and depending on rotational coordinates, 
allows for perturbations only between such terms (if permitted by the other 
selection rules) as are to be found in the same or directly neighboring rows of 
Table I. The other rules are that AA = 0, +1 and between the © states the 
+< — rules. No perturbations can appear between states belonging to 
different columns; thus for instance the II states belonging to (¢, 7) and 
(7, ¢) configurations do not perturb one another, neither do those belonging 
to (7, 6) and (6, 7) and so on. 


6. Relations between Matrix Elements and the Constants of Molecules 
Calculating the diagonal elements of operators (1) and A (L.S) and equating 

them we obtain a relationship between the experimentally determinable 

multiplet splitting constant 43, of the term II and the constants a: 


(40a) As = 43[(\’+1)4,—0'4)). 
TI 
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If apart from )’ the other quantum numbers of the two electrons agree, i.e. 
ad, = G, =a, then: 


(400) Asn = Sa. 


Here the vibrational wave-functions play no roles, since the wave-functions 
of the same state occur on both sides of the operator. 

Relations between d, and a,, and further between a, and a,, appearing in 
(40a) and (38) might be given for An;, = An,, = 0, that is if terms belonging 
to the same electron configurations are perturbing each other. In such cases 
the relation is supplied by (38) where a(ji,, mip) = Gp and a(Mxq, Meg) = og. 
In spite of the fact that in R,,,(p) and R,..(p) all the principal quantum 
number 2 of the electrons is the same, the two functions belong in general to 
two different sets of orthogonal functions. As a result, owing to the different 
values of the vibrational quantum numbers, the integral is in general not 
zero, but according to the magnitude of the ‘‘overlapping”’ region, it has some 
value less than unity. 

If, however, the potential curves of the two states perturbing each other 
lie in such a manner that the vibrational wave-functions of the two interacting 
states have practically no overlapping regions (i.e. if the difference between 
equilibrium internuclear distances is very large) then the value of the above- 
mentioned integral will be very small and, between the two terms, no observ- 
able perturbation occurs, not even in a case where every selection rule is 
satisfied. 

On the basis of what has been said, it can be stated that for the constants 
occurring in the matrix elements between the different components of two 
states perturbing each other the relations previously described will hold 
without question, and these perturbations might be compared with one another. 
On the other hand if two different perturbations are compared, for instance 
the (1S *Z) and (‘II *IL) perturbations or the (‘2 *I) and (2 4II) perturba- 
tions, and if the constants appearing in these perturbations are to be expressed 
by the multiplet splitting constant (Budé and Kovacs 1954), then, even if 
the electron part of the constants occurring in these is the same, it can be 
done only if the position of the potential curves is known and adequate values 
for the integrals of the vibrational wave-functions are taken into consideration. 

The second part of the paper will deal with the perturbations between 
terms of even multiplicity in a three-electron system. 
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ON THE SPIN-ORBIT INTERACTION IN 
DIATOMIC MOLECULES. IP 


I. Kovacs 


ABSTRACT 


This paper deals with perturbations between states of even multiplicity in a 
three-electron system. With the aid of the matrix elements which have been ob- 
tained, a comparison of the perturbations in the various components of a given 
state and in some cases a comparison of perturbations in two different states is 
possible. Also for certain cases it has been possible to establish relationships 
between the matrix elements which give the magnitude of the perturbations and 
experimentally measurable molecular constants of the perturbed states. 


The first part of this investigation (Kovacs 1958) dealt with the various 
perturbations between states of odd multiplicity of a two-electron system 
taking into consideration the spin-orbit interactions. The present paper treats 
the perturbations between states of even multiplicity in a three-electron 
system. 

1. EIGENFUNCTIONS 

If the electrons of the molecule (with the exception of three) form closed 
shells, the states of the molecule can be characterized by the quantum numbers 
of the three electrons in question in the following manner: (jy, lip, Ntp, @ ip; 
Nea legs Neqs Fkqi Mir» Liry Nar» Grr) OF briefly: (Ap, Ag, A,). The notations p, g, 7 
signify that the ith electron is in a state denoted by p, the kth in a state 
denoted by gq, and the /th in a state denoted by the index r, where (7, R, 2) 
can take up every possible permutation of the indices 1, 2, 3. Since here we 
intend primarily to deal with interactions between the ¥ and II states, we will 
deal with only two of the many possible electron configurations, and the 
wave-functions of these two only will be considered. These will be (A, A, ¢) 
and (A’, \’, ); this means in both cases |A,| = |A,| and in the first case A, = 0, 
whilst in the second, A, = 1. As is well known in the case of a (A, A, «) con- 
figuration the following states are possible: *2*, *2**, *Z+, 7(2A), 7(2d)*, 
4(2\), where for \ = 0 the =~ states are lacking whilst in the case (A’, \’, 7), 
the possible states are: 711, ?11*™, 211, 271*@), 2(2x’—1), 2(24’—1)*, 411, 
‘IE, 4(2)’—1), 2(20’+1), 2(2A’+1)*, 4(2X’ +1). If A = 0 or X’ = O the terms 
2\ or 2A\’—1 and 2\’+1 states are not present whilst if \ = 1,2... the 
symbols (2A) denote A, T,... states in case of \’ = 1,2..., the symbols 
(2\’—1) denote II, ¢,... states, and finally the symbols (2A’+1) denote the 
¢,... states. The states designated by the sign * arise from the various 
possible orientations of the spin momenta of individual electrons, whilst the 
states distinguished by bracketed numbers above on the right-hand side arise 
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from the different orientations of the \’s. It follows from the latter that for 

= 0 in case of a (6, o, 7) configuration the *II™ and *II® states as well as 
the further ones distinguished by “ and “ coincide; thus in such cases only 
*11, *1I*, and ‘II states are present. If \’ = 1, altogether six *II and three 
‘IT states occur, since in such cases the configurations (2\’—1) give also II 
states. In every other case four II and two ‘*II states occur. 

The electronic part of the wave-functions of the = and II states resulting 
from the above electron configurations can be given in Hund’s case (a) accord- 
ing to the method described in I. For the ?2* states of the (A, A, ) configura- 
tion they will be 


o(*2#)) = (1/V6){ d1(73*) (a1B2— Biae)a3— $2("E*) (a183— Bias)ax 

(1) + $3("E*) (a283— Bras)ar}, 
e(S*,) = (1/V6){ ¢1C3*) (B1a2—a182)B3— $2(°E*) (B1a3—a183) 8» 

+ 3("E*) (Boa3—a283) Bi}, 


where 
o1(S*) = (1/V82°) [pigetpoqilrs eee 
(la) 2(°*2*) = (1/V8n° ) [Pigs Paqilre sin * Mo— D3), 
$3(**) = (1 //8n°) [Peds pagelri — A( g2— $3), 
and for ?2*+ states 
o("2*4,) = (1//18)if o1(> ty * ) (Byara3+a1B2a3— 2a 10283) 
— ¢2(" ‘yet ) (Bax; — 2a) 8203+ a 10283 ) 
(2) =a ¢:(*2**) ( — 2Bia20;+018203+a1a283) ts 
o(°2**)) = 1//18) { o1(°2"*) (a18283+810283— 288203) 


— $2(°E**) (a18283;—2810283+ 618203) 

+ $3(°2**) (—2a18283+B10283+8 Boa) }, 
where all of the functions $;(?2*+), ¢2(?2**), $3(*2**) correspond to the 
functions occurring in (la) with the difference that the signs * occurring 
there have to be changed to * signs. 

For the wave-functions of 42 states we obtain 
fant tay 
o( 23/2) = o( >*)aya20;, 


o(‘=#;) = (1 //3)¢(3*) (a1a283;+a1B203+ Bia2a3), 


(3) co ok 
(223) = (1/V3) bC2*) (a18283+B10283+ 818205), 
@(*S*3)2) = o(**) 818283, 

where 

(3a) o('E*) = (1/V3)[b1C2"*) — b2°2"*) + bs7E**)]. 


If \ = 0 in equation (la) and (3a) the normalization factor will be 1/47 


instead of 1/+/8r* and in such cases only the expressions having the ‘‘cos’’ 


factor will differ from zero, i.e. the X~ states are lacking. 
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For the 711%, *1I®, and 2(2\’—1) states of (A’,\’, 7) configuration the 
forms of the wave-functions formally correspond to those of functions (1) 
with the difference that in place of the expressions $;(?2+), ¢2(?2+), @3(?2+) 
for *II,3/2 and *II,; states the functions ¢i:(71I*), d:2(?II*), ¢i3(7TI*) and for 
*11’_, and ?II’_3,2 states the functions ¢i:(?II-), ¢i2(?1I-), ¢is3(7II-) are to be 
written, where 7 can be equal to 1, 2, 3 corresponding to 7II™, *II®, and 
*(2\’—1) states. The form of the above-mentioned functions is the following: 


$1: (71*) (1) 4a/ x") [pagaet™ 9+ gape “Ole! 
(4a) o12(“II*) (1 /4n/n*)[pigae*™ © +-a:p:6° PO ght 
$13(TI*) _ ( /44/ x) [pigoe™™ @! PL gpge OP yg 


The @(7I*), d22(7TI*), ¢32(*11*+) functions can be obtained from the fore- 
going if we write —2’ instead of \’. For a (X’, \’, 7) configuration the functions 
1 = 3 will be the functions of the II state for \’ = 1. Thus 


$31(II*) (1 /4/ 2°) [pogstqops|net ere, 
(4b) 32(MI*) = (1/40 2") [pigs tqipslre® te, 
33(II*) 7 (1/4/ w°) [Pige+aipe)rset Oth. 
The wave-functions of the 7II*™, 211*@, and 2(2\’—1)* states can be obtained 
from the formulae (2) if for *1I*43,2 and 7II*,, the functions $;(?2**), @2(?2*+*), 
$3(?=**) are replaced by $i (?*1I**), $:2(?1**), dis(?T**+) and for *11*’_y and 
*T1*’ 3,2 by du (II*-), d2(I*-), o:3(71I*-) functions where 7 can be 1, 2, 3 
corresponding to the *1I*™, ?11*, and ?(2d’— 1)* states. These latter functions 
can be obtained from (4a) and (48) if we put within the bracket (but not in 
the exponent) the — sign instead of a + sign. 
Finally the wave-functions of the ‘II states obtained are from (3) if for 
*T145/2, M432, 4144, 4M, the (42+) functions are replaced by ¢;(*I*) whilst 
for *11’,3, 411’, #11’ 3,2, *1I’_5,2 they are replaced by ¢;(‘II-) functions, where 


2 


II 


Il 


5 a 2% 2% 2..% 
(4c) $:(M*) = (1/V3)[bn CT *)— b2(M*) + ¢0( *)). 
Taking into account the two possible orientations of the resultant 2 relative 
to the axis of the molecule as well as the vibrations and rotations of the nuclei, 
the complete wave-functions for Hund’s case (a) will be the following: 


ah ¥(2%) = (J ‘/2)[b(2F0)uy,40+ 6(2*0)u, ol R,, vy 
7 v(ZE)! = (1/4/2)[6(ZE0)t03,40— 6(Z4a)ty,—0] Rr. 
where for °?X states 2 is equal to 1/2 and for 42 states it is equal to 3/2 or 1/2. 
Also 
(5b) W(t) = (1/V2)[o(Mya)uy.40+ O(Mo 9) uy, o] Rn.» 
where for *II states 2 is equal to 3/2, 1/2 and for ‘II states 5/2, 3/2, 1/2, 
—1/2. 

Further the eigenfunctions for Hund’s case (8) are 
fv CZFy) = ¥CZP)’, 
lyC2#y) = vCxp); 


(6a ) 











(6c) 


(6d) 
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ae ee Ne gf ED ge | 
| WC 2743/2) = gee 741 ¥ (2372) + 741 2j 
| = — by ('Et2) — 4/3 v(O34) 
1 3(J— J+3/2 
I¥C2F4y) = jg (: D ycss, 2)’ + 4 = v(‘2F) / 
| = —3V3 ¥(' 22)’ — 40 2y )s 
fyi = ai /3(J+3/2) 4nt aes ys fot | 
| y( “J 1) = 2 V J+1 v( D572) T+1 v( >y) 
~ —3V3 C252) + CED), 


| : Hf FR og 3(J—1/2) 
| ¥CZFa2) = ly PE WBE) — i 3 ‘sty | 


| = — $0 (2h)! + V8 VEY: 
> ae FHP age 
¥ (TF, ) os i - 2( oe 1) v(t) — Vz J+} ) ¥(CTG ) 


= J+3/2 J—}— 
Lyme) = 4/ S88 vent. + sq Vem 


| 
| = /\/2) [WCIF2) — CIE), 
| = (1/+/2) [vCm#.) + ¥ Cif): 


ly pee et a 


/3(J—4) (J+5/2) " 3(J+5/2) 
+ WV (FE (T4372) ¥ CIES, oe Tu v(‘TIF) 


= —1/V8[y (ME 2) + 39 ‘Nt 2) +30 Cnp) + ¥(‘44)], 


Pree | /3(J—3/2)(J—3)(J+5/2) 

pre = a V so+na+say Vita 

| (J+9/2) (J—=4)_ = — 

| 7 V5 TEN 43/2) 1M) + T+ 4y VOD 

| +4/ aut ») cnr, pee [(--V3¥('ne 2) — v(m, 

| +4) + vi ‘II)), 
Ieee Ly -3/2)(J+3/2)(J+5/2) | 404 

ee ae G+DG-nuty — YM) 


| (ta PY T+SD ‘ 7B (J45/2) cages 
+W FG Hoty MM + V GEpTEH MD 


3(J+3 
~¢™ Gt ) wns) |= 4 Jal--v3 v (‘ME 
| + (‘IF 2) + ¥( - — V3 ¥(‘T*))]. 
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3 (J+3/2) (J -5/2) 
|- Td—3)  ¥¢ He) 


'3(J—3/2)(J+3/2) 
+ / — ry) ¥ (“TF 2) 


= al-v (Me) ld ¥(‘Hi2)—/3 v4) + vO). 


In case (a) and in case (6) the same formulae will hold for doublet states 
denoted by * except that on the right-hand side everywhere the eigenfunctions 
denoted by * are to be written. For the perturbations AA = 0 it seems adequate 
in Hund’s case (6) to use the simpler approximate formulae of the right-hand 
side. 

Using the above wave-functions as well as the operator equation (1) of 
Part I some perturbation matrix elements can be computed. 


2. PERTURBATIONS AA = 0 
For AA = 0, i.e. in case of perturbations between ~ states or II states, 
only the first term of operator I(1) supplies non-vanishing matrix elements. 


(a) ?=—*Y Perturbation 

Between 2 states of the same symmetry but different multiplicity there 
is no perturbation (in other words at most only second-order effects can 
occur). Thus for each component in the first approximation we have 


(7) H(? d+ 43+) = 0. 


For & states of different symmetry a different matrix element applies according 
to whether the wave-functions (1) or (2) are used for the ?2 state. The follow- 
ing expressions are obtained for the neneOn matrix elements: 


(8) V3 HCZ5 4 “Zisa) = —HC ZF ‘ZF4) 

™ H(*S zy “4 "SiH) = —V/3 H(? 75 = 32) = +hi/3/2 T; 
where in case of wave-functions (1), 
(9a) 7(2B+ 4EF) = \(ap>+a,), 


whilst in case of wave-functions (2), 
(96) r(2E*+ 4EF) = (1/+/3) A(a,—a,) 


(see Fig. 1). The a, and a, constants have similar meaning as in Part I and 
they will be dealt with later in the discussion of the selection rules. Since 
both states belong to Hund’s case (6), the matrix elements connecting levels 
that violate the selection rule AK = +1 are zero. 


(b) 2=—?3' Perturbation 

Between 22 states of the same symmetry the electrostatic interaction and 
certain terms neglected at the separation of the wave equation supply a 
matrix element independent of the rotational quantum number J according 
to the selection rule AK = 0. 
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Teh Ith I-% I-% 





Fic. 1. 


Between °> states of different symmetry the operator I(1) gives matrix 
elements of first order of magnitude; these are 


(10) H( zz, 255) = ACSF, “E5y) = FG/2V8)r 


J-} 


(see Fig. 2), where, according to whether we use the wave-functions (1) and 
(2) fora (A, \, ¢) configuration, the value of 7 will be 


sa r(2z+ 27) = 0, r(2E+ 22**) = r(a,+a,), 
- r(2Bet 2E*F) = (1/4/3)M(ap—ay) 
E ay’ 
a 
Jt’ 
J-% 


Ith JI-h 





Fic. 2. 


Hence the *2* or ?S~ states corresponding to the wave-functions (1) can 
perturb only the °2*~ or ?2** states corresponding to the wave-function (2), 
whilst ?2**+ or ?2*~ states can interact with each other as well as with ?2~ or 
°Y* states. The selection rule AK = +1 holds throughout. The selection rule 
relative to the states belonging to electron configurations of different X's 
will be dealt with in the following section. 

(c) 42—‘d’ Perturbation 

Between ‘> states of the same symmetry the perturbation terms previously 
mentioned in case of ?Y states give matrix elements independent of the 
rotational quantum number J according to the selection rules AK = 
Between components with AK = +2 only second-order effects occur. 

For different symmetries I(1) gives the following matrix elements: 


(12) 2H (‘ZFis/2 *Z544) = VW3ACIH, ‘2F_,) = 20 (SH, ‘Shse) = tir 


5 








KOVACS: SPIN-ORBIT INTERACTION. II 335 


(see Fig. 3), where 


(13) r(4Z+, 4-) = (1/V3)A(a,>—a,). 
" 
E eae 
Je% 
I-% 
." 1% 


eh Ih eh F-% 





Fic: 3. 


In Fig. 3 perturbations of the same magnitude are denoted by the same sign. 
Within the selection rule AA = +1 the operator I(1) allows for an inter- 
action only between ‘+ and ‘=~ states belonging to the (A, A, «) configuration 
of the same X. 


(d) *1I—‘I Perturbation 

Here various cases must be distinguished, on the one hand owing to the 
different *II and ‘IT wave-functions belonging to the (A’, \’, 7) configuration, 
and on the other hand owing to the fact that *IT and ‘II states can occur in 
either Hund’s case (a) or (6). Therefore the perturbations are first investigated 
in general according to the various Hund’s limiting cases, and later expressions 
of the matrix elements are given for special cases. 

(1) *II(a)—*II (a) perturbation.—In this case for both the *II(a) state and 
the ‘II(a) state the wave-functions (5) must be used. Thus we obtain: 


(14) HUF. ME») = —W CMF WH) = (1/6) + 


(see Fig. +). For these perturbations the selection rule AQ = 0 is characteristic 
since both states contributing to the perturbations belong to Hund’s case 


(a). 





Fic. 4. 
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(2) *I1(6)—*II (a) perturbation.—In this case using the simplified form of 
the wave-function (6c) for ?II(b) states and the wave-functions (5) for the 
4II(a) states we obtain (see Fig. 5): 


(15) H(nt,, Nt.) = H(ne,, ‘MF = ACUF, ‘WF2) 
= —H (mt, ‘1*) = (1/2V3)r. 


2 
Tt) “ 












Fic. 5. 


(3) *II(a)—‘II() perturbation.—Here we take the wave-functions (5d) for 
the *II(a) states and the simplified form of (6d) for the 4II(6) states (see Fig. 
6): | 


(16) HCI. “Ws) = /3 HM. ‘Mt,) = —/3 HC. ‘ME,) 
= —H(M¥. ‘Mt 3) = HME ‘MWe,s2) = —/3 Hie ‘nF,,) 
= —V3 A(t *Nt,) = HCP ‘Mt 3/2) = —3+. 


Ieh 
Joh 


I-k 





Fic. 6. 


(4) *11(6)—*II(6) perturbation.—Here we use the wave-functions (6c) for 
the 7II(4) states whilst the simplified form of the wave-functions (6d) is used 
for the *II(6) state and we obtain (see Fig. 7): 


(17) HCP, ‘Ws2) = —V3 HCnF, ME) = V3 WCF, ‘TF,,) 
= —H (Mt, ‘MW s2) = —(1/2V2)r. 


Since in this case both states belong to Hund's case (6) the selection rule 
AK = +1 will hold. 
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It may be noted if the ‘II state belongs to Hund’s case (a) only the two 
middle components will be perturbed, while if the ‘II state belongs to Hund’s 





Fic. 7. 


case (b) the two middle and the two extreme components will suffer different 
perturbations. Consequently if the *II state does not intersect the ‘II state, 
the two middle components will be continuously displaced relative to the 
two extreme components. On account of this displacement the ‘II term for- 
mula will not be exactly valid. This is the case for the ‘II state of the O2* 
molecule (Bud6é and Kovacs 1954). 

(5) The value of the r constant in different cases.—-As we have seen for the 
(A’, \’, r) configuration, four *IT and two ‘II states and for \’ = 1, six ?II and 
three ‘II states arise. Not all of these *IT and ‘II states interact with every 
other one. One can see that 


(18) r?N®, 41) = 2, 47) = 0 
unless 1 = k. This means that in several cases even between *II and ‘II states 
belonging to the same electron configuration no perturbations will exist. 
When i = k, the value of 7 can be different for the various interactions: 
721, 41 )) a —7(?®, 4] (2)) = d’(ap+a,), 
(19a) r?M*™, 417) = (1/4/3)[2a,—d’(a,—a,) ], 
r?2M*®, 4) = (1/+/3)[2a, +2’ (a,—a,)], 
whilst for \’ = 1, the following is added to these: 
(198) 72, 41) = a,—a,, r(21*, 411) = —(1/+/3)[2a,+a,+a,], 
where the meaning of a, is analogous to the values of a, and a, and its detailed 
form will be given in the next section. 


(e) *11—*II’ and *11—-*Il’ Perturbations 

For ?II—Il’ and ‘II—‘II’ perturbations if both states belong to Hund’s 
case (a), the electrostatic interaction and the terms neglected at the separation 
of the wave equation as well as the operator I(1) give matrix elements in- 
dependent of the rotational quantum number according to the selection rule 
AQ = 0. If both states belong to Hund’s case (6) then in the matrix elements 
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corresponding to the selection rule AK = 0, terms occur that are independent 
of J and have values determined by the electrostatic interaction and the 
certain terms neglected at the separation of the wave equation, whilst the 
perturbations between the components AK = +1 are given by the elements 
of operator I (1). 


3. THE VALUES OF THE CONSTANTS a,, az, AND a, AND THE 
SELECTION RULES FOR PERTURBATIONS AA = 0 
In addition to the selection rules already mentioned, by a more detailed 
consideration of the constants further selection rules can be obtained. The 


constants a are given by 


* * * 
(21 ) a; = f KR, sp) f Trier. lr. 72 Tr' tr. eed Pri. lip APn’ ip: Vip dr; 
* 
x f nig. Ug Qn’ kas Ugg AAT Rn. (p)dp, 


while for a, and a, the expressions I(14) will hold if the integral 


rx 
J Tnirrlrdr Tr'irsUtpeMr dr, 


is added. These expressions, similar to the previous ones, are diagonal in \ 
and according to the assumptions of Section 1.5 are diagonal in / as well; 
hence AX = Ad,, = 0, Alip = Alea = Ali, = 0. This means that in this case also 
an interaction can occur only between 2 states belonging to the same (X, A, o) 
configuration and between z states belonging to the same (X’, \’, 7) configura- 
tion. Thus a large number of the possible perturbations is excluded. 
According to the previous assumptions the expressions for a,, ad), a, can be 


written in an even simpler form: 
Ay = a(Niy, 2’ tp) f R*, R,. dp, 
(22) Be =A (Mee, ® tq) f R* oar, lds 
a, = a(n;,, 2 tr) f rR. 5 Rizr. dp, 


where each constant is diagonal in the other two principal quantum numbers 
not written out. 

More explicitly, the situation with regard to the principal quantum num- 
bers is the following: If two states are examined in which the principal quantum 
numbers of all the three electrons or only of two electrons differ, then each 
of the constants a,, a,, and a, will be zero and between such states operator 
[(1) does not give rise to perturbations. If two states differ only in the principal 
quantum numbers of one electron, two cases are possible: (1) An,, # 0 and 
An;, = An,, = 0 and thus a, ¥ 0, a, = a, = 0. In this case operator I(1) 
supplies no perturbing matrix elements between Y states, because in these 
only the constants a, and a, occur. (2) An,, = 0 and An;, = 0, Any, # 0 or 
An;, #0, An,, = 0, thus a, = 0 and moreover a, = 0, a, # 0 or a, ¥ 0, 
a, = 0. In such cases every perturbation dealt with until now is possible, 
with the restriction that in the matrix elements one of the constants is always 
zero. Finally if the principal quantum numbers of all the three electrons are 
equal in the two states, i.e. if the states belong to the same electron con- 
figuration, we have An,, = Ang = An,, = 0, and therefore a, # 0, a, # 0, 





KOVACS: SPIN-ORBIT INTERACTION. II 339 


a, ~ 0, hence the relations given in the previous sections are fully valid. Only 
in this case can a relationship be established between the matrix elements and 
the multiplet splitting constants 42, and Ag. For a (¢, ¢, 0) configuration 
only ?2+, *2*+, 42+ states are possible. According to operator I(1) these = 
state exists. 


states cannot interact, since in such a configuration no 2 
If in a (A, A, 7) configuration with \ # 0 two electrons are equivalent, only 
*y*-, 42> states exist. Among these, according to (1), (2), and (3), states 
of different symmetry can interact and in this case a, = ay. A perturbation 
can also occur between the above states and states belonging to another 
(A, A, ¢) configuration which differ from the previous ones only in that An;, # 0 
or An,, # 0, whereas if An,, ¥ 0, then a, = a, = 0 and each matrix element 
would be zero. 
If X = 0 and the two a electrons are equivalent, then only one ?=+ state 
arises and according to I(1) this state cannot interact with other > states. 
If in a (A’, A’, z) configuration the two electrons are equivalent, only one 
“TI, one ?1I*, and one ‘II state arise assuming that A’ > 1. For \’ = 1 a further 
"II state is added to these, whilst for \’ = 0, only one ?II state exists. 
Between states of the same configuration and \’ > 1 the perturbations 
dealt with previously are possible with the variation that everywhere a, = ay. 
In the case of perturbations between states of different configurations if the 
principal quantum number of the II electron varies, i.e. if An,;, # 0, and there- 
fore ad, = a, = 0, an intercombination perturbation can occur only between 
*11* (but not 71) and I states. If only the principal quantum number of one 
of the X’ electrons differs, i.e. An;, # 0, or Any, ¥ 0, and therefore a, = a, = 0 
or a, = a, = 0, each of the perturbations which have been dealt with can 


24 
= 


’ ’ 


occur but two constants are always zero. Since for \’ = 0 there is only one 
“II state, this state cannot interact with any ‘II state. 

Summarizing what has been said before for AA = 0 in Hund’s case (a) 
operator I(1) gives matrix elements: 


(23) H|(nip, Biss Nip; CO ips Nkgs bias Aka CKaqs Nir, Dae ) TER Tir; n, v, A, >) 


, , , ee , , oe , , Pee eel os 
xX (n soy ips r ips O ips N kas Vin, ke oO kqs n pel E aac rN iss F te, By DT, A > = )] x 0 


if among Anjy, Any 7, An;, at least two are zero, Aly = Ak, = Ali, = 0, Ary 
Adxrqg = Adi, = 0, i.e. AA = Oand Aoi = Aor, = Ao;, = 0, 1.e. AZ = O. 
In Hund's case (6) instead of AY = 0, the relation AK = +1 holds. (For 
states of the same A, i.e. AK = 0, only the electrostatic interaction and 
certain terms neglected at the separation of the wave equation give matrix 
elements, whilst operator I(1) does not give any.) 


1. AA = +1 PERTURBATIONS 
In these cases the bracketed part of operator I(1) supplies the perturbing 
matrix elements. Here the following cases can occur: 
(a) ?2—4II Perturbation 
(1) ?23+—4II(a).—If the ‘II state belongs to Hund’s case (a), 


(24) V3 HEF, 14) = WC2F 4 “Ty) / | 
= \/3H (st, ‘y) = —H (st, ‘M,) = -77 2/3 
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(see Fig. 8), where the 2 exponent is 0 or 1 according to whether we are 
state. The value of the constant £ appearing here 


dealing with a 22+ ora? 
as well as later on will be dealt with below: 


y 
~_ 
. 
_ 


Be 
ey 








Ta) 





: . hh 
A I-22 
iF 





Fic. 8. 


(2) *2+—4II(b).—If the ‘II state belongs to Hund’s case (0), 


(25) 2H (ZF, “Ut,s)2) = -V/3 H(2F 44 “T4) 
= 24/3 H(*x#,, ‘W#,) = 2/3 H(2t, ‘W) 
V3 HC zt, ‘MF_) = 2HC?2¢, ‘MF_a/2) = 17 (2/V3)E 


(see Fig. 9). Here as in the following figures the perturbations of the same 
order of magnitude are denoted everywhere by the same sign. Since both 
states are in case (6) the selection rule AK = 0, +1 will hold. 


De 





(b) 42—*II Perturbation 
Here again two cases can be distinguished: 
(1) 4=4—*I1(a).— 

(26) ~/BH (‘2d Md) = \/BA(Th22 4) = HC BE, M2) 

= —3H(‘2t,, m4) = H(‘2E, Wi) = 3H(‘22., MN) 
= V/3H(‘23,3/2 Min) = —\/38H (24,372 1) = — (07/2) 


(see Fig. 10). 
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IH Joh bh Ty 





Je5h Io% J-% I-53 


io 





(2) 4>+—2]11(b).— 
(27) 2H(‘S$.s,2 ME,) = 2\/3H(‘2,, ME,) = V/3H(‘2H,, ME) 
= \/3H(‘st, M4) = 2\/3H(‘2E, My) = 2H (“Sd ,9/2 M4) 


= —1°(2/V3)E 


(see Fig. 11). 





ely InN I-% Th 
Wh eh I-% I-% 





(c) *2—*II Perturbation 


Between Y and II states of the same multiplicity the operator I(1) gives 
matrix elements independent of J, whilst the terms neglected at the separation 


of the wave equation give elements depending on J. Operator I(1) can supply 
terms non-diagonal in A and ¥ unlike operator A (LS), as can be seen in detail 


in the next section. 
(1) ?3+—2II(a).— 
HC2#,4 WF) = HCZFy M2) = —P2n»/ (I-43) (J+3/2) , 
(28) = HC 3544 1) = 2° 1&+2n(J+3/2)], 
HS, M$) = i*[&-2n(J—1/2)] 
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(see Fig. 12). The value of the constant 7 will be dealt with in the following 
section. 


ar 





Fic. 12. 


(2) ?34—211(6).— 
rc J+3/2 
av [e+4(J+4)] 


= op let an FI 2)(J+3/2)] 


H(*3 Li “"TIF,5) 


= re: 
(29) H(23., TG) = iy: oe ‘at 
i ‘ J4+-3/2 i” 
,» Ty) = -— = ——~f 
H(?> = IT 744) V2 yee 9 § Vo 
9 oi ; ing ie ee 
H (=F, Ty4) = V2 VV 743 [§—4n(J+3)] 


wicad ‘ [é—4nr/ (J—1/2)(J+1/2)] 


see Fig. 13). For perturbations with AK = 0 the matrix elements depend on 
J whilst for perturbations with AK = +1 they are independent of J and have 


approximately the same order of magnitude 


he 





Meh I-% 





Fic. 13. 


(3) The A-type splitting of the *II state in the limiting cases.—*Il states can 
be perturbed by ? and ‘4 states. As can be seen from eqs. (26) and (27) the 
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4X states, owing to the equality of the corresponding matrix elements, cannot 
cause the splitting of the *II rotational levels of + and — symmetry. Thus the 
A-type doubling in *II states can be caused only by the more distant 22 states. 
Applying the general formula described in I(24), formulae which can be used 
in the two ideal limiting cases and even in the neighborhood of them can 
be given as follows: 


In Hund’s case (a) 
hAv(7II3,2) = 0, 
(30) ; 
hAv(?My) = 2(C,+C2)(J+34), 


whilst in the neighborhood of case (6) 
hAv(?ly43) = Ci(J+3/2)+ Co(J+ 3) (J+3/2), 


(31) 
hAv(*My-3) = Ci(J— 3) —C2(J— 3) (J +3), 


where C; and C, are identical with the expressions given in I(26a) with the 
difference that doublet terms are to be written everywhere instead of triplet. 


(d) *>—*II Perturbation 

For this perturbation special attention has to be paid to the elements of 
operator I(1) that are non-diagonal in A and ¥. For the configurations under 
consideration, operator I(1) supplies only a third of the contribution supplied 
by operator 4 (L.S). Here smaller values are obtained owing to the fact that 
‘Y as well as ‘II states are in intercombination interaction with other states 
and the remaining two thirds of the interaction is consumed there. (For details 
see the next section.) 

(1) ‘s*—‘*TI(a).- 


zs ,f J+5/D0-3/0~3 
H(‘3%,3). ‘Mze) = i7y V (J+5/2)(J—3/2)(J—>3) 


J+1 / 
= inv (J+5/2)(J—3/2) 
ee 4 + De 5 
H (2513/2 Ts;2) = ey ht ns 2 [e4+-69(T+5/2)] 
= si [§+6n(J+5/2)] 
tint tt, _ © ,/I+3/2 P 
HT ( — J+3/2 T+; ) = a 4 J+1 [§+3n(J+5/2 2) } 
= — le+3n(J+5/2)] 
/ T+3/2 


H(‘34,3)2 lt) = 5 a = [€+2n(J+5/2)] = 5 [E+2n(J+5/2)] 


J+1 
/3(J+5/2)(J—3/2)(J—3 


H(S$4 M2) = Po VY -— 


= inr/3(J+5/2)(J—3/2) 








344 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


oe 


H(‘23,4 ‘Td ) = ae Nieponst9/2)) = © [E+2n(J+9/2)] 


H(*2544 ‘4 = F 4/ (T+8/2 §—39(J—3/2)] = [$-3n(J—3/2)] 


J 
; (Taal 2 i” 
H (2b TE) = 554 =[E+6n(J+})] = - a [E+6n(J+3)] 
ci Me ae )(J+3/2)(J=3/2) 
(32) H(CZ74 Ts;2) = tn / I+) 


ai *9VJ/3(J+5/2 2 (7-3/2 
ae eS I4ae . * 
H(*3F-4 “TW2) = = [E-2n(J—7/2)] = a n(J—7/2)] 
gf [E+3n(J+5/2)] = y E+3n(J+5/2)] 


4 r . J Bk cots 
HC33 ‘T,) = “fa ile-m —6n(J+3)] = ~3v8 [§—6n(J+4)] 


H (‘3 3/2 Wey) = in aaa nee 


H(*3#, Wy) = - 


= Par (J+5/2)(J-3/2 


H(‘s3 = 3/2 “Tiys) = 24/3 V 7 = [é- 6n(J—3/2)] 
i” 
= 5073 [E-Gn( J -3/2)] 
| ® ,/F=3 
H (‘25a ‘W}) = ——Jy Af 7? E-3n(J-3/2)] 


fe 
= — 9 e-3n(J-3/2)] 


itn 


[E-2n(J-3/2)] = 5 [e-2n(J-3/2)] 


ee 
H(*2 Zy-3/2 = 1) = ar ag 


(see Fig. 14). According to these operations each perturbation is possible 
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because the *2 state belongs to case (6) and the ‘II state to case (a). Every 
matrix element depends on J and those which are approximately the same 
are connected in the figure by a thin line. 


(2) ‘>+—‘TI(b).— 


4 — i* /JI+5 ke 

2742/2 32) =—- = J 
HT ( J+3/2 IT 7+-3/2) V/2 VV J4+3 [é +4n(. +3/2 )] 

— 5 etn J+5/2)(J+3/2)] 

(Awe 4 F = A as J(I+3)_ = _ PP ‘ 
ACh My) = 5 V 3741)U+3/0) * A 
H(‘st, 3/2 ‘at 1) = 0 
H(‘s4, 3/2 “ 2) = 0 

a ee ae ae I(I+5/ 2) ie a ‘ 

H('2F44 M7432) = /2 V3 (J+1) (J+3/2)*~ ~ye* 


f J+9/2 é ~~. | 
ni et Any/ (J J+3/2) 
v [oplttes{+uvon (J+ 


scat lite inn/(I+3, —_ 
V2 


H(‘s 2. i "TIF,4) 


1 
H('zby ‘nd ) = #2 4/6 ae e = 


(33) H(‘33,, "TIF_-3/2) = 0 
H(z. } Ta = 0 


‘ 2 V2 (J=43)(J4+3/2",_ 2 v2 
HS34 M4) = — 3 2 / J50G ay ee es 


sai 2 J—7/2 eer 
H(3$ ‘TG 4) = + [dle pan VOFDT-D) 


case g ; Cri’ = 
~ a he tn /(J+3) (J—3) 
H(3# ME sa) = 5 4/ GEDI=8P) , 

a Re Viv 3J(J-31) * VW 
H (‘SE 3/2 “M,3;2) = 0 


H(‘3 2 3/2 “THE, ;) = 0 


ic J+1)(J+3) i 
B(*ZF- 2 “TE y= - a ‘- Hn ee i= es 
, £ 3/2 
H(‘23-a2 ‘W-a2) = 5 le 4a) 
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(see Fig. 15). According to the selection rules valid for case (6) the matrix 
elements AK = +2, +3, are zero, the matrix elements AK = +1 are inde- 
pendent of J, and four of the six perturbations are approximately equal while 
the other two have a different value but are approximately equal to each other 
(in Fig. 15 these are connected by thin lines). The matrix elements AA = 0 
depend on J and differ from one another. 


IeVe Ir J- I-%y 





Teg Hh J-% J-% 





Fic. 15. 


(3) The A-type splitting of the *II states in the limiting cases..-According to 
what has been said in points (a) and (d) perturbations can be caused in a 
‘II state by 2+ and 4Z+ states. (In our model *S* states do not occur.) The 
>~ and 4+ states lying on one side of the ‘II state and the ?2*+ and 4 7~ 
states lying on the other side will give rise to a shift in the same direction. 
Taking into consideration what has been said above according to I(24) as 
well as (24) and (32) for Hund’s case (a) we obtain: 


hAv(‘IIs/2) = 0, 
hAv(‘I3,2) = 0, 

(34) hAv(*I,) = ((4/3)C1 + 4C2)(J+34), 
hAv(‘Il_,) = 0, 


and in the neighborhood of Hund’s case (6) according to (25) and (33): 


4 1. J+5/2 5 = = 
hAv( IT 743/2) = alo Fit + Cy\(J+5 2) + C.(J+3 2)(J+5 ah 

‘ lL. J-82 <p ; ‘ . sas 
hAvC y+) = go —, 3Ci(J+9/2) — C2(J+3)(J+3/2), 


(35) lL. J+5/2 - i ae 
hAv(‘My-) = Co 3 — 4C,\(J—7/2) + Cx(J-3)(J+4), 


. J—3/2 
hAv(*Iy—3,2) wa 30 == : =i + Ci(J—3/2) = Co(J—3/2)(J—3), 
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where Cy, Ci, and C2 correspond to the value given in 1 (26a) with the difference 
that the triplet and singlet notations should be replaced by the quartet and 
doublet notations. 

As an example of equation (35) it may be mentioned for instance that for 
the A-type splitting of the II state of the O.*+ molecule the difference between 
the experimental and the theoretical values at J = 40 (extrapolated) is small. 
[| Av(4IT743/2) = 0.84 cm.~! (theor.) and 0.56 cm.~! (exp.), Av(*Iy43) = —0.33 
cm. (theor.) and —0.30 cm. (exp.), Av(4IIy-3) = 0.03 cm.~ (theor.) and 
0.02 cm. (exp.), Av(*Hy—3;2) = 0.04 cm. (theor.) and 0.00 cm! (exp.) 
(Budé and Kovacs 1944).] 


5. SELECTION RULES AND THE VALUES OF THE CONSTANTS 
—& AND n FOR AA = +1 PERTURBATIONS 
The values of the constants € and 7 occurring in AA = +1 perturbations 
are the following: 


, 1 ; . i ’ 
(36a ) E[(AAa) (AAT) ] = 2/2 | ie... tp) | O¥ntec tar 0B fit grt n't. Ver 1 4T1 


x j Praca lip.d Pn’ ip. Uip.d dr, j a lig. A In’ kg. Ukgd dreRy’, » (p)dp, 
: 1 ' ‘ 
(365) n[(Qdo) AAm)] = 2/2 f Re, s(o) | nie tir OFF 8)’ tr tr t1ET 1 


x j Pi, Iip.A Pn’ ip. Wey dr dr; f oie. leg. ¥ Yn’ kq. lig.d d7,B (p) Ry’. v’ (p)dp, 


where B = h?/82°yup? and yu is the reduced mass of the molecule. 

For matrix elements containing £, owing to the integration carried out pre- 
viously over the spin-coordinates, two of the spin quantum number pairs 
(Gin, O ip)y (Oxq, O'eg)) (Grr, 777) have to be always equal to one another, and 
where the two members of a pair are different from one another, the difference 
can only be +1. For matrix elements in which the constant 7 occurs, all three 
pairs of spin quantum numbers have to be equal. Asa result of the integrations 
over ¢ (as can also be seen from the expressions of the constants) one finds 
that AX = Oand Ad,, = +1. For \ = 0 the expressions (36a) and (366) have 
to be multiplied by 2. 

The preceding selection rules expressed by the resultant quantum numbers 
of the molecule are the well-known selection rules: for £ AA = +1 and AL 
= #1, ie. AQ = 0, and for 7 AA = +1, AY = 0, i.e. AQ = +1. 

For the perturbations dealt with in Section 4 a more detailed examination 
of the integration over yg and the spin coordinates shows that not every pertur- 
bation from among those which satisfy the above general selection rules is 
possible. Thus in the case of perturbations between terms of different multi- 
plicity the 2? state described by the eigenfunction (1) cannot interact with 
any “II state, hence /7(°X, ‘Il) = 0; only between ?* states having the eigen- 
function (2) and ‘II states is an interaction possible. Similarly between *II 
and 45 no interaction exists, i.e. A7(7M, 42) = 0, but there is interaction 
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between *II* and ‘2 states. In the case of states of the same multiplicity no 
perturbation occurs between ?Z and *II* states nor between ?=* and *II states, 
ie. Z(?2, *11*) = H(?>*, 211) = 0. This statement holds not only for matrix 
elements in which the ~ constant occurs but also for those in which the 7 
constant occurs. 

There is a difference between the (*2, *II) and (?2*, *II*) perturbations. 
For the matrix elements (*2, 7II) perturbations containing the constant £ 
the operator I(1) gives the same values as the operator .1(LS) whilst in the 
case of perturbation (*2*, ?II*) as well as of (42, ‘II) a difference appears. 
More explicitly: H(?=*, ?II*) = 3 H(?z, 711). In general it can be stated 
that between two states of the same multiplicity which owing to configurational 
reasons are not in AA = +1 intercombination interaction with other states 
(for instance *Z and *II) the operators I(1) and A (LS) give the same matrix 
elements, whilst between two states of the same multiplicity which can be in 
AA = +1 intercombination interaction with other states (for instance 
?S*, *11* or 4, 411, as well as from I the *2, *II states) smaller values are given 
by operator I(1) than by operator A (LS). The 71 and ‘1 states occurring 
for \ = 1 cannot interact with any of the ?Z, ?2*, 4D states. 

If the assumptions of Section 1.5 are introduced, additional selection rules 
can be obtained. For the non-vanishing elements we have An;, = An,, = 0, 
and Alj, = Al,, = Al;, = 0. In such cases the equations (36a) and (376) 
can be written in even simpler forms 


(37a) E[(AXo) (AXm)] = <6 la(f+e)1. R*,o(e) Rn’ (p)dp, 
(376) n{(Ado) (AAT)] = a [F+eh-f Rs.(0)B(o) Rev (0)de, 


where 
(38) (a(f+g)], = fr. ur.o(P a 8 )a(r1) (fitgetn're, Vir.t1(%y 8) sin 8,7; dr dd. 


Here the distance of the electron from the origin has been designated by 
F in order to distinguish it from the eigenfunction 7. 

If according to [(29) and I(30) in the r eigenfunction the two variables are 
separated, we find 


(39) [a(f+g)], = F a(n, n',) V1 (L+1), 


where the expression a(n,,, ”’;,) is analogous to those occurring in I (35). 

As we have seen before for the principal quantum numbers only An,, may 
differ from zero. If it is different from zero the constant a(2,,, n’;,) cannot be 
related to the multiplet constants. If however An,, = 0, then the — and 7 
constants will take the following form: 


ELORe)(MAA)] = F 575 a VIATED, 
(40) ./9 — 
n[(Ado) (AAw)] = =a Ben/I41), 
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where 
(41a) a, = a(ni;. Nir) f R*, »(p) Rn. o(p)dp, 
(41d) Bsn = f R*.o(p)B(p) Rn, o(p)dp. 


Summarizing what has been said before for AA = +1 in Hund’s case 
(a) operator I(1) gives 


FA (10 55s Lips No Fins Megs bags Nv Tags Mier Lees O, O73 2, V, A, Z) (Mens D5) Ay O55; 


up? ip? 


(42) 
Nias lags Xs Thgi Mrs Urys 1, 03,3, 0, A+1, ZF1)] + 0 


if for the principal quantum numbers An; = An,, = 0 (An;, may be zero 
but might also differ from it), and if Alj, = Ak, = Al;, = 0; AX = 0, AdA;, 
= +1 that is AA = +1; furthermore from among Ao;j,, Aox,, Ao;, two are 
always zero, whilst the third is equal to +1, that is AY = 1 and thus AQ = 0. 

In Hund’s case (6) instead of the selection rule AY = +1 the selection 
rule AK = 0, +1 holds. 

In Hund’s case (a) for matrix elements depending on J and neglected at 
the separation of the wave equation the same selection rules are valid except 
that now Agip = Aox, = Ao,, = 0, hence AY = O and thus AQ = +1, whilst 
for Hund’s case (8) selection rule AK = +1 takes the place of AY = 0. 

In Table I below the states between which perturbations can occur are 


indicated. 

















TABLE I 
Electron 
config. Molecular states 
lot 2+ = 

¢ ¢ 4 = : = : >a 

27, (1) 2,4 (2) ony * (1) oy * (2) 477 ()) 4,,(2) 
oot i =i; , ae = B e ih =e 

Jet 2e- Jan + 2a *— dent 4a 

"TO a a = . a p 

2y_ (1) 2y_7(2) 247 (3) 2 (1) 2p gp (2) 2p *(3) 4 yy (1) 4y_(2) 477 (3) 
nur i; Ge 3 Oy Se ee a ee ee gy Ee in 





The connecting lines in Table I connect states which can perturb one 
another except that perturbations between states in the same row are not 
indicated in order not to overcrowd the diagram. The (2 and II) states 
occurring here—with the exception of those indicated—cannot interact with 
(II and 2) states of another configuration. Thus an explanation is found for 
cases in which perturbations otherwise possible do not occur. For values 
? = 1,2 the matrix elements 7(2, II) differ only in sign from each other. 

If two electrons are equivalent, then Table I must be replaced by Table 
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II. The first two rows illustrate the case of one valence-electron, because two 
equivalent o electrons form a closed shell in the molecule. 











TABLE II 
Electron 
config. Molecular terms 
2 2a 
og = 
2 2 1 2 2 
o TF o.- = nn” 
ro a > = > = Scare 
2 2. 1) 2 2 2 3 2H 277 * (2. 4 1 4 2 
nx Rate? ee? Ee ee. ee eS wate 





6. RELATIONS BETWEEN THE MATRIX ELEMENTS AND THE 
CONSTANTS OF MOLECULES 
If the diagonal elements of the operator I(1) and .1(LS) are calculated 
and equated, relations are obtained between the experimentally measurable 
multiplet splitting constants of 7II and ‘II states and the constants a. More- 
over with the knowledge of these the regular or inverted character of the 
II and ‘II states can be established in certain cases. (The places where this 
can be done unambiguously are denoted in Tables I and II.) 
Ay =a, Ame) = 3[—a,4+2A(a,—a,)], Asn) = 3[a,+A(G,—a,)]; 
Ay) = 4, Ame) = 3[—a,—2A(a4,—a,)], Atm = 3[a,—A(G,—a,)]; 
Ary) = —a,, Ames) = 3[4,+2(a4,+4,)],  -14ns) = 3[—a,+ (4,+4,)] 


where the values of d,, @,, and a, correspond to the a constants on the right- 
hand side of (22) if there the principal quantum numbers distinguished by a 
prime are equal, similar to those without a prime. If two states for which 
Ani, = Any, = An,, = 0, i.e. states belonging to the same electron con- 
figuration, perturb each other, then the difference between the constants of 
the matrix elements and of (43) can arise only from the fact that in (43) the 
integral of the vibrational eigenfunctions has been assumed to give exactly 
unity, whilst the actual value of the integrals occurring in the matrix elements 
depends on the relative position of the potential curves of the two states. 

In the case of two equivalent electrons d, is equal to d, and thus they do 
not occur in (43); in such cases for each II state it doubtless can be decided 
whether it is a regular or an inverted state. 
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BAND STRUCTURE OF RHOMBOHEDRAL GRAPHITE! 


R. R. HAERING? 


ABSTRACT 


The band structure of rhombohedral graphite has been investigated using 
the nearest-neighbor tight-binding approximation. The resulting behavior of 
the z-bands near the Fermi surface is more complex than in the case of the 
Bernal stacking. The two z-bands still touch, but the touching points no longer 
lie on the edges of a hexagonal prism in k-space. Instead, they lie on cylinders 
whose axes are the edges of the hexagonal prism. The radii of these cylinders 
are proportional to yi, the nearest ‘out-of-plane’ exchange integral. The 
de Haas — Van Alphen effect in the rhombohedral structure may be expected to 
vield useful information about the magnitude of ¥:. 


1, INTRODUCTION 

The band structure of graphite has received a great deal of attention in 
the literature (Wallace 1947; Coulson and Taylor 1952; Lomer 1955; Slonc- 
zewski 1955; Slonczewski and Weiss 1955, 1957; Johnston 1955, 1956; Corbato 
1956). Most of these investigations are extensions of the work of Wallace 
(1947), who used a nearest-neighbor tight-binding approximation based on 
the 2p, orbitals. All of the above calculations assume that graphite crystallizes 
ia the Bernal structure (Bernal 1924). In this structure, the carbon atoms 
within a single layer of graphite form a plane hexagonal lattice (assumed to 
be the xy plane), in which the atoms are separated by 1.42 A. The layers are 
then placed on top of each other in a sequence which we may denote by 
1—-2-1-2...., indicating that alternate layers have the same projections on 
the xy plane. The interlayer separation is 3.35 A and the unit cell contains 
four carbon atoms. 

The X-ray and electron diffraction spectra of graphite show a number of 
extra lines which cannot be attributed to the Bernal lattice. These extra 
lines were first observed in the electron diffraction spectrum by Finch and 
Wilman (1936), who explained them in terms of diffraction from very thin 
crystallites (thickness ~4 unit cells). Taylor and Laidler (1940) and Lipson 
and Stokes (1942) also observed such reflections by X-ray diffraction. The 
sharpness of the X-ray lines suggested that one was dealing with much larger 
crystallites, so that the explanation of Finch and Wilman (1936) had to be 
rejected. Lipson and Stokes (1942) showed that the extra reflections were 
consistent with a rhombohedral structure of graphite, which differs from the 
Bernal structure in that the stacking sequence is 1—2-3-1-2-38.... Line 
intensities suggested that graphite was a mixture of the two structures in 
the following ratio: 80% Bernal, 14°% rhombohedral, 6% disordered. Lukesh 
(1950, 1951) found furthermore that the bond angles were not quite equal 
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and that impurities introduce a superstructure with a 15-ring periodicity 
which could also yield extra lines. However, his conclusions were tentative, 
and these points have not yet been settled conclusively. 

All of the effects mentioned here may be important in a complete discussion 
of the band structure of graphite. Here, however, we shall be concerned only 
with the effect of the rhombohedral modification. 

Hennig (1957) has pointed out that natural graphite is a mosaic combination 
of the 1-2-1-2 and the 1-2-3-1-2-8 structures, the latter comprising perhaps 
5-10% of the total crystal. On heat treatment, the rhombohedral structure 
disappears almost completely. This finding is in accord with what one expects 
theoretically. The binary mixture of the two graphite phases may be treated 
in terms of a one-dimensional Ising model (Newell and Montrell 1953). An 
exact solution is possible, which yields the result that one of the two phases 
disappears almost completely as the temperature is raised. If the crystal is 
then quenched, a sample of nearly pure Bernal-type graphite is obtained. 
Hennig (1957) pointed out that it would be possible to make samples con- 
taining an appreciable fraction of the rhombohedral structure, so that the 
energy bands of this structure could be investigated. Since these bands depend 
on the same parameters as do those of the Bernal structure, one may expect 
to learn something about the numerical values of these parameters, particu- 
larly about the strength of the interplanar interaction. 

In the next sections we will outline a nearest-neighbor tight-binding calcu- 
lation based on the 2/, orbitals, in a manner completely analogous to Wallace's 
calculation (1947) for the Bernal structure. More detailed investigations of 
the Bernal structure (Coulson and Taylor 1952; Lomer 1955; Slonczewski 
1955; Slonezewski and Weiss 1955, 1957; Johnston 1955, 1956; Corbato 1956) 
have shown that the results of Wallace (1947) for two-dimensional graphite 
are unaffected by a more general treatment. However, for the three-dimen- 
sional structure, a more general treatment admits the possibility of ‘‘vertical”’ 
overlap of the m-bands, while in the Wallace treatment these bands touch, 
but do not overlap. The fault lies in the nearest-neighbor approximation. In 
this approximation, the vertical edges of the Brillouin zone are lines of con- 
stant energy, a condition not required by the crystal symmetry. A variation 
of energy along these edges results in the vertical overlap mentioned above. 
By analogy, we may anticipate the manner in which the band structure of 
rhombohedral graphite would change if the nearest-neighbor approximation 
were abandoned. 

2. RHOMBOHEDRAL GRAPHITE STRUCTURE AND BRILLOUIN ZONE 


The crystal structure of rhombohedral graphite is shown in Fig. 1. The 
interlayer spacing is c = 3.37 A, and the spacing between similar atoms in 
the same plane is a = 2.46 A. The unit cell for this structure contains only 
two atoms. The cell we choose is shown in Fig. 1. Relative to the cartesian 
coordinate system shown in Fig. 1, the primitive lattice vectors may be 


chosen to be: 
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oa CS sq). 

22> (22, Aa, 
a 

aed (ine). 


The volume of the unit cell is thus 


Ve = ler: (e2X3)| = $a2eV/3. 





Fic. 1. Structure of rhombohedral graphite. 


From the definition of the primitive vectors for the reciprocal lattice, viz. 


b,-,; = 276;;, 


1 S 
= 2 — - _--- 
b 2(. Vo a eer” 
1 1 & 
eee ee meen 
b: (15 tae: sg 


b; = 24( 0, 0, . 
c 


The volume of the first Brillouin zone is thus 


Vp.z. = |bi- (beXb;)| = 823/V,. 


one easily finds: 


The surface of this zone is formed by the perpendicular bisecting planes ot 
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the following set of vectors: +b;, +bi, +b., +(bi+b.+b;), +(bi+b:,), 
+(b,+b;), +(b.+b;). The resulting zone is shown in Fig. 2. The coordinates 
of the 24 corners of this zone are given in Table |. These corners are labelled 
1 to + according to their k, value. The six corners corresponding to a fixed 


ky 





Fic. 2. Brillouin zone of rhombohedral graphite. 


k. value are counted 1 to 6 starting at the positive &,-axis, and numbering 
counterclockwise. Thus, each corner has two labels; the first giving the k. 
value to which it belongs, and the second giving its position relative to the 


k,-axis. Some of the corners have been labelled in this manner in Fig. 2. 


3. BAND STRUCTURE OF RHOMBOHEDRAL GRAPHITE 

The object of a band-structure calculation is to give at least a qualitative 
description of the energy bands for the conduction electrons. By analogy 
with the Bernal structure, we may expect that the bands of interest in the 
rhombohedral structure are those formed by the 2, orbitals. We shall restrict 
our discussion to these orbitals, so that our problem will involve one electron 
per carbon atom and two atoms per unit cell. In particular, we shall be inter- 
ested in knowing the shape of the energy surfaces near the highest energy in 
the first Brillouin zone, since this region corresponds to the Fermi surface. This 
may be seen as follows: The density of states in k-space is uniform and is 
given by 


N(R) dk = (2V/8r*)a®k. 


Here V = VV, is the volume of the crystal and V is the number of unit 
cells. The number of states in the first Brillouin zone is just (2VV./82°) X Vg. z. 
= 2. Thus, there are 2.V states in this zone, or one state per atom. It follows 
that at 0° K. the Brillouin zone will be just filled, so that the Fermi energy 
is the highest energy in this zone. 

Following the usual tight-binding treatment, we choose an approximate 
Bloch function of the form 


(1) v = bitrds, 
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where 

(2a) o: = Do exp(ik-r4)x(r¥—-ra), 
A 

(2b) o: = >. exp(ik-rg)x(r—Frz). 
B 


In equations (2a) and (20) the sums are over all 4 and B atoms, and x(r—r,4) 
and x(f—fg) are 2p, orbitals centered on A and B atoms respectively. We 
now substitute (1) into the Schrédinger equation 

(3) Hy = Ey. 

Multiplying this equation in turn by ¢:* and ¢:* and integrating over all 
space, we find (if we neglect overlap of 2p, orbitals centered on different 
atoms): 

Hy,+ALM 2 = NE, 

Hy+d\He = NE. 


(4) 


In equation (4) we put: 


(5a) Ay, = Ho 


ll 


¥ explik: (ta—a)) fx re Hx e— rade, 
is e 


(56) Hy =H = > explik: (re—ra)] { x (r—r,4)Hx(r—rg)dr. 
A.B e 
It should be pointed out that Wallace (1947) also used the relation Hy; = Ho». 
Carter and Krumhansl (1953) have pointed out that, in the case of the Bernal 
structure, these energies are actually different by an amount of the order 
of y1, the nearest out-of-plane exchange integral. However, this observation 
does not apply to the rhombohedral structure, for which the relation 1,,;= He» 
is exact. 
The energy resulting from equations (4), by elimination of X, is 





(6) E=~ Hus| 5 Hy | . 
We now perform the double sums involved in (5a) and (56) by keeping only 
nearest and second-nearest neighbors in and out of the graphite planes. Con- 
sider, for instance, equation (5a). We fix the value of r4, for the moment and 
perform the sum over A keeping only terms for which (r4—Tf4-) has one of the 
following values: 0; +e), +e, +(41—2); es, +(43—t2), &(e3—4:). The 
sum over A’ then multiplies our result by NV. 

The various exchange integrals which occur in the above manner are listed 
in Table II. The notation corresponds exactly to that used by Wallace (1947). 
The various minus signs in Table II are convenient, in order that all y’s be 
positive. Keeping only the terms listed in Table II, one finds: 


(7) Hy = Eo— wy? cos(ak,)+4 cos( Wiehe) cos( 24+) 


‘\ ‘ 9 


+70 2 cos( tHe ck.) +4 cox( #2) cos( tte - cas) t , 











358 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


TABLE II 


NEAREST AND NEXT NEAREST NEIGHBOR EXCHANGE 
INTEGRALS 























(f4—Ta’) Sx*(r—ra)Hx(r —1r4)dr 
U Ey 
+71, £T2, +(t1 —2) —yo' 
+3, +(t3—%1), +(*3—2) vn’ 
(Tx—Ta) Sx*(r —ra)Hx(r —rp)dr 
4(%2—2%)) \ 
3(t1+%2) ( sli 
(41 —2t2) 
3(Ti +2) —Ts y1 


(t1—2t2) +75 | 
i(%2— 281) +5 
— 4(t1: +2) +43 


— (ste) _ ( ict) (2) 
(8) oo — Yo EXP V3 +2 exp 2/3 cos\ “5/5 


J ak, iak, k,\\ ’ 
+y, exp — its) 4.2 exo( ae ) cos( , +71 exp(—1ck,). 


The inclusion of the y;’ term presents no difficulty, but greatly complicates 
the analytical form of the /(k) relation. From now on, we shall ignore this 
term completely. It may be shown that the inclusion of this term would not 
alter our basic conclusions about the energy surfaces near the Fermi level. 
Equation (6) then becomes 


3ak 2, \ | 
(9) E(k) = Ko— 02 cos(ak,)+4 ene Sate) cos( tH) 


7 e 4 \ 

* la I+4 cos'( #) 4 | coal #) cos( ists 
‘ : ; 4 17! 
i raycos( tte tok.) +2 cos( te) cos( ack ) | ; 


The minus sign applies to the inside of the first Brillouin zone (valence band). 
The energy at various points in the Brillouin zone may now be computed 
from equation (9). Consider, for example, the energy at the center of the 
zone, O. Putting k; = k, = k, = 0, one easily finds: 


E(O) = Eo—6y0' +[38y0—1]. 
In a similar manner, the energy at each of the 24 zone corners is 
E(C) = Eo t2.865y0'+[0.3375y0+y:]. 


We may also evaluate the energy at the center of each of the 14 faces of the 
Brillouin zone. Denoting these points by fz, Pex, Fry depending on whether 
the face in question is a rectangle, a regular hexagon, or an irregular hexagon 
respectively, we find for the energies at these points: 
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E(Fr) = Eot2yo' [vot] (six faces), 
E(Fra) = Eo—6y0' [870+ 71] (two faces), 
E(Frxz) = Eot2y0'+[vo—71] (six faces). 











The significant feature of these results is that the two z-bands are separated 
by an energy of the order of yo. None of these points in k-space are of great 
interest, since they do not correspond to the highest valence band energy. 
This maximum energy is principally determined by the magnitude of jp, 
which should be as small as possible. In order to discuss the E(&) relation near 
the maximum valence band energy, it is convenient to introduce a volume in 
k-space which is bounded by a hexagonal cylinder whose height and volume 
are equal to the height and volume of the Brillouin zone shown in Fig. 2. The 
edges of this cylinder are given by: 


_ (22 2s ) . ( 4m ) 

Si (2s ’ 3a oe al 3a rhe}, 
_(=2 2m ) a (3 2 ) 

ve (33 ’ 3a Re)» Nae Se a 

ote ) a (2s. =2n ) 

“= (0, 3a ke) = avV/3" 3a she} - 


From Table I it is easily seen that the actual Brillouin zone would go over 
into this hexagonal cylinder if all a?/c? terms were neglected (i.e. in the limit 
c— o). Let us now evaluate the E(k) relation at the edges e of this hexagonal 
cylinder. It is readily shown that the energy along these edges is independent 
of k. and is given by 

E(e) = Eot3y0' +71. 


Since yo > Yo’, y1, this is by far the highest V-band (and lowest C-band) 
energy found so far. However, there is still an energy gap between these 
bands equal to 2y;. The question now arises whether there are any points 
in the Brillouin zone where the gap is zero (i.e. Hi. = 0). An exact solution 
of the relation H7,;. = 0 shows that there are such points, and that they lie 
very near the edges e of our hexagonal cylinder. We may therefore find these 
touching points by expanding the E(k) relation about the edges «. We put: 


k, = k,*+K,, k, = k,*+K,, k, = K,. 
Furthermore, we shall use the notation of McClure (1957) and put: 
y/3a[K2+K,?] = 3\/3a|K] = o, 
K,/|K| = sina, K,/|K| = cosa, K,c = £, ¥i/Yo = B. 
Energy will be measured from the point Eo+3y0’, so that we put Eo +370’ = 0. 


By expanding equation (9) in the indicated manner for each of the edges 
e,(4 = 1,2,...6), one finds: 


(10a) Ex, (0, @, £) = +yolo?+8?+28ecos(E—-a+¢,)]} (i = 1,3,5), 
(106) Ex, (0, a, £) = +y0lo2+82+2Becos(E+a+¢,)]} (i = 2, 4, 6). 
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In these relations, ¢; is a phase angle which is different at each edge: 1 = 7/6, 
63 = 52/6, 65 = 32/2; ¢2 = —7/2,¢4 = —71/6, be = 7/6. It is easily verified 
that only that portion of our hexagonal cylinder for which 


1 a 
—4n oi < +49 


lies within the first Brillouin zone. 

From equations (10a) and (100), we see that the bands will actually touch 
near ¢ = OQ, but only in special directions (fixed a). The condition for touching 
bands is 
(11) o°+6?+26c0cosé = 0, 
where @ is any one of the angles £+a+¢;. Solving for o, we find: 


(12) o = —Bet"?, 


For a physical solution, we require that: 


o = real, oi 0, 
which implies that 
(13) 6=7, G-=p. 
For any € in the range —3m < — < 3a the bands touch at six points in 


the Brillouin zone. These touching points lie on six cylinders whose axes are 
the edges «; and whose radii are ¢ = 8. The loci of these touching points 


are shown in Fig. 3. 





Fic. 3. Lines along which the two z-bands touch. The lines lie on cylinders ¢ = 8, whose 
axes are the edges e;. (The drawing is not to scale.) 


It may also be verified that the energy contours in a plane — = constant 
are circular about the touching point. If ¢ is o measured from the touching 


point instead of from an edge e€;, one finds: 


(14) Exccont(G) = 4 ye = £4\/3ay0\K —Kol. 


g 


In equation (14), Ao denotes the A corresponding to the touching point. 
Typical constant-energy contours for such a plane of constant & are shown 
in Fig. 4. The broken circles in Fig. 4 represent the cylinders o = 8. As & 


changes, the energy contours move along the broken lines (cf. Fig. 3). 
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Fic. 4. Typical energy contours in a plane € = constant. Broken lines denote the cylinders 
o = B. (The drawing is not to scale.) 


4. DISCUSSION 


The band structure of rhombohedral graphite which emerges from this 
simple calculation is surprisingly complex. The two z-bands still touch, as 
they do in two-dimensional graphite or in Bernal-type graphite in this approxi- 
mation (Wallace 1947). However, the touching does not occur at points of 
special symmetry of the zone. In fact, one sees from equation (13) that it 
is the value of the exchange integral y; and not any condition of symmetry 
which determines the location of the touching point. For a non-zero value 
of yi, the energy contours for energies near zero are highly warped cylinders, 
whose axes are the loci of the touching points (cf. Fig. 3). The warping becomes 
more and more severe as y; is chosen larger and larger. For y; = 0 the 
warping disappears and the energy contours are regular cylinders whose axes 
are the edges ¢€;, as in the case of two-dimensional graphite. The energy 
surfaces for E > y: also display only a slight warping. 

It seems that the complexity of this band structure may furnish a useful 
criterion for the magnitude of y:. Consider first a sample of rhombohedral 
graphite for which Ey < y;. Let us see what we may expect to be the result 
of a de Haas — Van Alphen effect measurement in this case. The periods of 
the susceptibility fluctuations are a measure of the areas of cross-section of 
the constant energy surface E = Ey and a plane perpendicular (in k-space) 
to the direction of the applied magnetic field. If the magnetic field is applied 
along the c-axis in graphite, one period will be observed in the de Haas — Van 
Alphen effect, since the six cross-sectional areas in question are all equal and 
correspond to circular disks (cf. equation (14)). However, if the magnetic 
field is tilted away from the c-axis, six different periods will result, because 
of the warping of the energy surfaces. 

Consider next a sample for which Ey > y:. In this case, the warping of 
the surfaces is very much less severe, so that essentially one period (i.e. six 
very nearly equal periods) will result in the de Haas — Van Alphen effect, for 
any orientation of the magnetic field (the period will of course be a function 
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of the direction of the field, but for any direction, only one period is present). 
If we may assume that the relative magnitudes of Ey and y; are the same 
for rhombohedral graphite as for Bernal-type graphite, then the case of 
Ey &y.1 corresponds to the McClure model of graphite (McClure 1957), 
while the case of Ey > y1 corresponds to the model of Haering and Wallace 
(1957). Actually, the model of McClure (1957) involves holes as well as 
electrons, owing to the ‘‘vertical’’ overlap mentioned in Section 1. If such 
overlap were accounted for in the present calculation, we would expect that 
the loci of the touching points (cf. Fig. 3) would no longer be lines of constant 
energy. This alteration does not markedly affect the above argument, which 
is based solely on the warping of the surface. The presence of holes as well 
as electrons would however increase the number of distinct oscillations to be 
expected (but equally for all orientations of the field). 

We conclude that a measurement of the de Haas — Van Alphen effect of 
rhombohedral graphite as a function of orientation of the magnetic field 
could yield in a convincing manner the relative magnitudes of Hy and 1. 
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MECHANICAL STRENGTH OF ICE FROZEN FROM 
AN IMPURE MELT! 


E. R. PouNDER 


ABSTRACT 


Ice frozen from a melt containing a small amount of a light alcohol, ketone, 
or ether is much weaker than normal ice, for a limited freezing exposure. The 
weakening effect is enhanced if the melt also contains traces of one or more of a 
variety of long-chain, organic substances, one of the most effective being sodium 
carboxymethyl cellulose. Inorganic salts, on the other hand, appear to be largely 
ineffective and may even cancel out the weakening action of the organic solvent. 
It is suggested that the organic additives are concentrated in the surface of the 
melt and remain as liquid layers between the ice crystals, thereby causing a 
large reduction in the shear strength of the ice sheet. Adsorption at the surface 
is discussed as a mechanism for concentrating organic substances in the surface 
layer and rejecting inorganic salts into the interior of the solution. 


INTRODUCTION 


When a solid is grown from a melt containing traces of impurities, it is 
well known that the crystal structure of the solid and its mechanical, electrical, 
and other physical properties may depend strongly on the nature and con- 
centration of the impurities in the melt. In the case of ice, Truby (1955) has 
shown that one fluorine ion per 5X10' water molecules alters the micro- 
structure of the ice drastically. Smith-Johannsen (1946) measured the adhesion 
of ice frozen on to a waxed aluminum plate from 10~* molar solutions of 
various salts. He found values for the adhesion ranging from 4250 for distilled 
water through 1300 for magnesium chloride solution down to less than 100 
for solutions of strontium chloride and other salts, values of adhesion being 
quoted in grams per square centimeter. 

Moesveld (1934, 1937) showed that if minute quantities of certain organic 
materials having long carbon chains were added to the surface of still water 
as it was starting to freeze, the resulting ice consisted of largely non-coherent 
crystals so that the shear strength of the ice sheet was negligible. He reported 
that the most effective chemical additives were sodium and ammonium salts 
of stearic and alginic acids and mixtures of these salts. 

Inorganic salts, particularly chlorides of calcium and sodium, have been 
used for many years to combat icing problems but since the useful property 
of these salts is the depression of the freezing point the quantities which 
must be employed are large. In addition, the corrosive properties of salt 
solutions (which need not be stressed to anvone driving a car in Montreal) 
virtually rule out their use in connection with aircraft. The use of organic 
additives offered possibilities of the need for smaller quantities of anti-icing 
material and freedom from corrosion problems, and this study was undertaken 
to confirm and if possible extend the investigation of Moesveld. 

‘Manuscript received October 24, 1957. 

Contribution from the Ice Research Project, Department of Physics, McGill University, 
Montreal, Quebec. 
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METHOD 

Numerous freezing tests were carried out in a commercial deep-freeze unit 
which could be set to temperatures in the range —5° to —30° C. The tem- 
perature in the box was measured with a resistance thermometer and it was 
found that the actual temperature varied within a range of about two degrees 
on either side of its mean value, as the refrigeration equipment cycled under 
the control of its thermostat. Freezing was carried out in glass or plastic 
tanks, insulated on the bottom and sides with 1 in. of a foamed plastic with 
a thermal conductivity A = 0.25 B.t.u./ft.2/hr./° F./in. For the freezing 
times involved, this insulation ensured that heat transfer was virtually limited 
to the top surface of the water or ice, as in the freezing of natural bodies 
of water. 

The normal procedure was as follows. The substance being tested was 
dissolved in water or in water and alcohol. The freezing tank was partially 
filled with distilled water and placed in the freezer to cool, the water being 
stirred frequently until it was at a uniform temperature below 4° C. but with 
no ice present. A quantity of the solution was then sprayed uniformly over 
the surface of the water from a plastic squeeze-bottle. Several polyethylene 
rings of 10cm. internal diameter were placed in the tank, resting on the 
bottom and projecting slightly above the surface. The tank was then allowed 
to freeze for the desired exposure, usually from 40 to 100 degree-hours. The 
exposure was calculated as the product of ambient temperature in the freezer 
(degrees Celsius below zero) multiplied by the time the tank was left in the 
freezer after the solution was added. 

When the tank was removed from the freezer, the disks of ice contained in 
the polyethylene rings could be removed easily from the tank and tested. The 
thickness of the disks ranged from about 5mm. to 1.5¢cm., depending on 
both the exposure and the nature of the additives used. Each disk was placed 
centrally on a padded ring stand with an internal diameter of 5 cm. A padded 
brass pin of area 10 cm.? was placed on the center of the disk and loaded with 
a lever arrangement until the ice broke. The loading pressure for fracture 
of the ice (Pr) was used to compare samples. It is a satisfactory relative 
measure of strength, but is difficult to relate to the ultimate tensile and shear 
strengths. For any series of tests the degree-hours of freezing must, of course, 
be held constant. Despite the relatively large area of the loading pin relative 
to the ice sample, fracture patterns observed were typical of those found 
when a “‘point”’ load is applied to an infinite sheet of ice. Fig. 1 is a photograph 
of the assembled pieces of a disk which had been fractured. It will be noted 
that most of the breaks occur along radial lines. With some of the weaker 
ices, as the load increased the first crack was a rough circle centered about 
the point of application of the load. A small increase in load then resulted in 


complete failure, again along radial lines. 


RESULTS 
Initially, freezing tests were carried out using solutions of various chemicals 
in water only. The tests were of two types: volume tests in which the additive 
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Fic, 1. Fracture pattern of a disk of ice loaded centrally. 
Fic. 2. Lower surface of ice sheet containing alcohol and soluble starch. 

Fic. 3. Lower surface of ice sheet containing alcohol, polyvinylmethyl! ether, and sodium 
carboxy methy! cellulose. 
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was dispersed uniformly throughout the volume of the water before freezing, 
and ‘‘surface’’ tests in which a technique, described under Method above, was 
used to ensure a Maximum concentration of the additive in the surface layer 
of the water. For a given quantity of additive, the surface tests were from 
two to five times as effective in reducing the strength of the ice, showing 
clearly that we are dealing with a surface phenomenon. Results were dis- 
appointing with water solutions, the resulting ice being at best only 30° 
weaker than pure ice. However, it was found that traces of alcohol (of the 
order of 5 milliliters of ethyl alcohol applied to a surface of 800 cm.?), either 
by itself or in combination with other additives, gave extremely weak ice. In 
some cases the strength was less than 10% of that of normal ice. 

Table I lists the results of a series of tests. One gram of each additive was 
dissolved in 75 ml. of water and 25 ml. of industrial ethyl alcohol. Twenty 
milliliters of the resulting solution was sprayed on a tank of surface area 
800 cm.? and containing 2.5 liters of distilled water. The concentration of 
additives was thus 2.5 g. of solid and 50 g. of ethyl alcohol per square meter 
of surface. In all cases the freezing exposure was 60 degree-hours at an ambient 


TABLE I 
FRACTURE LOADING PRESSURES OF IMPURE ICE 








Additive in 25% alcohol solution Pr Average % deviation 
Pure water (no alcohol) 1.000 9 
20 ml. of 25% C2:H;OH 0.2638 18 
20 mi. LiCl 0.139 15 
20 ml. methyl cellulose 0.180 8 
10 ml. sodium stearate + 10 ml. PVM 0.185 0 
10 ml. PVM + 10 ml. D-235 0.185 16 
20 ml. T-253 0.202 Is 
20 ml. PVM 0.211 17 
20 ml. ethyl cellulose 0.221 12 
10 ml. sodium stearate + 10 ml. D-235 0.224 6 
10 ml. PVM + 10 mi. L-245 0.255 20 
20 ml. soluble starch 0.262 5 
10 ml. sodium stearate + 10 ml. agar-agar 0.275 6 
10 ml. soluble starch + 10 ml. D-235 0.288 3 
20 ml. sodium stearate 0.302 29 
20 ml. agar-agar 0.325 29 

345 28 


10 ml. sodium stearate + 10 ml. sodium alginate 0: 


temperature of —25°C. This resulted in ice sheets of thickness 7 mm. for 
pure water and about 8 to 10mm. for various additives. The first column 
of the table lists the additives used. The abbreviations D-235, L-245, and 
T-253 refer to grades of Carboxel, the trade name used by Chemical Develop- 
ments of Canada, Ltd., for sodium carboxymethyl cellulose. PVM is an 
abbreviation for polyvinyl methyl ether. In the second column, the average 
value of Py, the loading pressure for fracture, is listed after being normalized 
to unity for pure ice. Each figure is the average of from 2 to 10 measurements 
and the third column gives the average percentage deviation from this average 
value, as a measure of the reproducibility of the tests. 
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The appearance of the impure ice differs greatly from that of natural ice. 
Pure ice is clear and glassy with perfectly smooth top and bottom surfaces. 
The impure ice is opaque and the top surface slightly roughened, but the 
most striking difference is in the bottom surface, which is irregular and often 
covered with stalactites and ridges of clear ice, ranging in horizontal length 
up to an inch or more. The ridges show no preferred orientation in either 
vertical or horizontal planes. They sometimes form irregular polygons bounding 
empty spaces in the ice which extend almost to the top surface. Figs. 2 and 3 
show two examples. Fig. 2 is a photograph of the lower surface of ice grown 
from a melt containing alcohol and soluble starch, whereas the additives in 
Fig. 3 were alcohol, PVM, and L-245. A polariscopic study of crystal sizes 
and orientations will be reported in a later paper. 

All the additives listed in Table I are organic compounds with the exception 
of lithium chloride. The peculiar effect of the lithium chloride was discovered 
in a series of freezing tests carried out with a variety of inorganic salts as 
additives, both in aqueous and ethyl alcohol solutions. These tests were 
undertaken to see if there was any correlation between the reduction of 
adhesion reported by Smith-Johannsen (1946) and cohesive strength as 
measured here. No correlation was found. The concentration of the salts 
used (various chlorides, nitrates, and acetates) was 6.2 g. per square meter 
of surface, some five times greater than in the tests summarized in Table I. In 
aqueous solution, lithium chloride reduced the strength of the ice by some 
20%, whereas the other salts had little effect, if anything usually increasing 
the strength slightly. Similar results were obtained with alcohol solutions. 
Lithium chloride reduced the strength sharply to about half of the strength 
obtained with alcohol as the only additive, whereas the other salts had little 
effect or increased the strength. 

Other variables investigated were concentration of additives, variation of 
exposure for a fixed ambient temperature, variation of ambient temperature 
for a fixed exposure, and the effect of substituting various organic solvents 
for the alcohol. 

To investigate the effect of concentration, equal weights of PVM and D-235 
were dissolved in a mixture of 75% water and 25% ethyl alcohol (by volume) 
in the proportion of 1 g. of solids to 100 ml. of solvent. Varying amounts of 
this solution were added to distilled water and frozen for 60 degree-hours in 
each test. Fig. 4 shows the variation of fracture loading pressure with con- 
centration, in milliliters of added solution. Curve A gives the results obtained 
when the additives were concentrated in the surface and Curve B when they 
were diluted uniformly throughout the volume of the freezing tank. The 
experimental spread of values is indicated by the lengths of the vertical 
lines. 

Using the same stock solution described in the last paragraph, a series of 
surface tests were made for different exposures. Five milliliters of the solution 
was added to the water and the tank was allowed to freeze at an ambient 
temperature of —25°C. Results are shown in Fig. 5. Curve A is a plot of 
relative fracture loading pressure vs. exposure, normalized so that unity gives 
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the strength of a sheet of ice formed from pure water in 60 degree-hours, and 
curve B shows the way in which ice thickness varies with time. Note that 
logarithmic scales have been used in this figure. Experimental spreads have 
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Fic. 4. Loading pressure for fracture as a function of concentration of additives. Curve A, 
surface concentration. Curve B, volume dilution. 
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Fic. 5. Effect of freezing exposure. Curve A, fracture loading pressure. Curve B, thickness 
of ice sheet. 


not been shown but are comparable with those in Fig. 4. Thickness could 
not be measured very accurately because of the unevenness of the lower 
surface of the ice sheets. If we let E stand for the exposure in degree-hours 
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and ¢ for ice thickness, the data of Fig. 5 can be represented by the approximate 
equations 
Py, = aE", t = bE®-®, 


where a and 6 are constants. The simplest theoretical expression for the rate 
of ice thickening is t « E®*, An empirical equation frequently used is 


(P+3t) « E. 


The index of 0.8 obtained here shows that this highly artificial ice thickens 
considerably more rapidly than ice frozen from pure water. 

A series of tests was made in which both pure water and water containing 
additives were frozen for a fixed exposure but at various ambient tempera- 
tures. The results indicated a very slight reduction of strength as the ambient 
temperature was increased. The reduction was only slightly greater than the 
uncertainty of the experimental measurements so the only conclusion to be 
drawn is the fairly obvious one that exposure is the major factor in deter- 
mining the strength of the ice sheet grown from a given melt. 

Finally, a series of solutions was prepared in which | g. of D-435 (a grade 
of sodium carboxymethyl] cellulose very similar to D-235) was dissolved in 
75 ml. of distilled water and mixed with 25 ml. of various organic solvents. 
‘‘Surface”’ freezing tests were then made with 25 ml. of the resulting solution 
added to 2.5 liters of distilled water and frozen for 45 degree-hours in a cold 
room at —15° C. Results are shown in Table II. 


TABLE II 
FRACTURE LOADING PRESSURE AND PROPERTIES OF LIQUID ADDITIVES 








Relative Melting point Boiling point Specific 
Additive Pr mn <. m ~ C. gravity 
Pure ice 1.000 
D-435 alone 0.750 
Methy! alcohol 0.062 YS 65 0.792 
Ethyl] alcohol 0.095 114 7s 0.789 
Isopropyl alcohol 0.065 90) 82 0.785 
n-Propy! alcohol 0.143 127 Q7 0. 804 
Isobutyl alcohol 0.162 108 108 0 802 
Ethyl ether 0.168 —116 35 0.708 
Isopropyl ether 0.475 60 6o 0.725 
n-Butyl ether 0.647 OS 142 0.773 
\cetone 0.114 OS 56 0.791 
Methyl ethyl ketone 0.060 87 80 0.805 
n-Heptane 1.01 a] 98 0.684 
Ethylene glycol 0.378 16 10S PES 
Methylene chloride 0.374 Q7 10 1.336 


DISCUSSION 


\ny foreign matter dissolved in water reduces the freezing point of the 
solution below that of pure water. The reduction is proportional to the number 
ol ‘particles’ (molecules or ions) present per unit volume. For equal masses 
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of organic and inorganic materials, the latter usually produce much greater 
lowering of the freezing point because of, firstly, their normally much lower 
molecular weights and, secondly, their high degree of dissociation in aqueous 
solution. For the concentrations of additives used in this experiment the 
over-all lowering of the freezing point is, at most, a small fraction of a degree 
but, on a microscopic scale, if a mechanism exists for concentrating the 
additives then it is possible for pockets or layers of solution to exist whose 
freezing point is lower than the temperature of the surrounding ice. 

It isa known property of ice that a growing crystal rejects impurities. Hence 
when two crystals grow towards each other the layer of liquid between them 
becomes a concentration layer for impurities. It is suggested that the reduction 
in strength of ice sheets observed in these experiments results from the pre- 
sence of liquid layers separating crystals and permitting them to slide fairly 
readily with respect to each other. These intercrystalline layers must usually 
be of microscopic thickness, although occasionally regions develop in the ice 
sheet, below the thin surface crust, in which evidently the melt contains too 
high a concentration of additives to freeze. These pockets may be several 
millimeters across and can be seen in Fig. 3. 

If this suggestion is valid we must look for an explanation of the ineffective- 
ness of inorganic salts relative to organic materials in weakening ice. A 
mechanism for concentrating organic impurities in the surface and inorganic 
impurities in the body of the solution is known to exist. Adam (1941) dis- 
cusses adsorption on liquid surfaces and points out that “organic substances 
with more than two carbon atoms for each soluble group in the molecule 
usually show marked positive adsorption’, i.e. tend to concentrate in the 
surface layer. Inorganic salts usually show negative adsorption. Adsorption 
phenomena are explained in terms of strength of intermolecular forces. Aqueous 
solutions of inorganic salts have a surface tension greater than that of water 
so the forces of attraction on the salt molecules or ions must be greater than 
those on the water molecules. Hence intermolecular forces tend to draw the 
salt into the interior. An illustration is the freezing of salt water. When normal 
sea water of salinity 35 parts per thousand is frozen the resulting ice has a 
salinity of only about 4 parts per thousand. 

Organic compounds usually lower the surface tension and hence are con- 
centrated in the surface. The technique of applying the additives in the 
“surface” tests tended to enhance this concentration. The water was cooled 
to its maximum density before adding the solution, which in most cases had 
a density less than that of the water. This is probably part of the explanation 
lor the effectiveness of alcohol as an additive. In addition, because of its low 
freezing point the alcohol must aid greatly in the formation of the inter- 
crystalline liquid layers. The tests summarized in Table Il were undertaken 
to try to sharpen the understanding of the role of the alcohol. Tests were 
made using two or more members of the homologous series of alcohols, ethers, 
and ketones plus a few assorted organic solvents. A tew conclusions can be 
drawn. The liquid organic additive should be soluble in water-—-heptane, 
which is virtually insoluble in water, had no weakening effect at all and the 
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ethers, which have limited solubility, were the least effective class of the 
solvents tried. The alcohols and ketones appear to be of the same order of 
effectiveness in weakening ice. Within a homologous series, there is a definite 
decrease in strength reduction of ice as one goes to higher members of the 
series. The reversal of this generalization for the two ketones is probably 
caused by the high vapor pressure of acetone, which may well evaporate in 
large part before the ice sheet forms to trap it. Isopropyl alcohol appears to 
behave anomalously, and the large difference obtained with the iso and 
normal forms of this alcohol is interesting. In this connection it is worth 
mentioning that the number of samples of ice tested in preparing Table II 
was small (in most cases only two for each solvent) and that the general 
experience in the experiments reported in this paper is one of poor repro- 
ducibility. Samples prepared, frozen, and tested under apparently identical 
conditions sometimes show deviations of up to 50% in the results. Individual 
results in Table I] might thus be modified if more extensive testing is done, 
but it is felt that the trends indicated are reliable. 


These experiments were supported by the Defence Research Board of 
Canada through D.D.P. contracts C.69-10-157 and C.69-700102. The 
assistance of P. Protopapas and F. G. J. Perey in performing some of the 
measurements is gratefully acknowledged. Mr. H. O. Stahl, formerly of 
Chemical Developments of Canada, Ltd., kindly suggested and supplied some 
of the chemicals used. 
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NOTE ON THE FORCED AND DAMPED OSCILLATOR 
IN QUANTUM MECHANICS! 


Epwarp H. KERNER 


ABSTRACT 


The wave equations for the forced, damped, and forced and damped oscillators 
are solved in closed form for an arbitrary forcing function; the solutions produced 
being in one-to-one correspondence with the stationary states of the unforced, 
undamped Hamiltonian Ho. The quantal motion is closely connected with the 
classical: for the forced oscillator the probability density is that of Ho but moves 
as a whole with the classical motion; for the forced and damped oscillator this 
motion is accompanied by a contraction progressing eventually into a delta 
function at the classical position. Transition probabilities between states of Ho 
are computed in the case of forced motion and depend solely on the classically 
acquired energy of the oscillator at any time. The transition probability vanishes 
strictly only when this energy has a value falling at the roots of a Laguerre poly- 
nomial associated with the transition. The classical dipole radiation emitted by 
a disturbed oscillator is, when the damping force is identified with the force 
of radiation damping, that of the classical oscillator: a shifted and broadened 
line. 


INTRODUCTION 

The customary development of perturbation theory for time-varying per- 
turbations ordinarily leaves unanswered the question of the long-time behavior 
of the perturbed system. Even supposing a perturbation calculation to be 
carried out to any desired order of accuracy, there remains the serious problem 
of the convergence of the calculation, about which little seems to be under- 
stood. We should like here to sketch the rigorous solution of the quantum 
motion of an oscillator exposed to an arbitrary time-varying external force. 
Special though the oscillator is, one may hope from a full knowledge of a 
simple problem to learn something of the nature of the answers to the more 
general questions posed above, or perhaps to replace these questions by more 
relevant ones. It will appear below, for example, that the effect of a perturbing 
force is intimately related to the corresponding classical motion, and it is not 
at all clear that discussing the perturbation in terms of transition probabilities 
between unperturbed states is in general very illuminating. 

Additionally it will be pointed out how the motion of the damped and 
simultaneously forced and damped oscillator may be solved. We shall have 
altogether a little isomorphism of the simple classical motions that are so 
instructive and their equally simple quantal counterparts. Our object: will 
be principally mathematical for the present, namely to demonstrate techniques 
for managing the wave equation; in another report it is planned to work out 


illustrative examples and physical consequences in greater detail. 


‘Manuscript received October 17, 1957. 
Contribution from the Physics Department, University of Buffalo, Buttalo, New York 
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FORCED OSCILLATOR 


For the Newtonian motion of an oscillator subject to an external force 
F(t) we have 


(1) mMéotkxo = Fit), 
and for the quantum motion, 

i? ay ‘ , oy 
2) ——— ——_; 1 xv —xF = - 
(2) —- aye + Laks vF(t)]y = th rm 


We may construct solutions to (2) that are in one-to-one correspondence 
with the unforced stationary solutions (F = 0) as follows. Write y = 
xexp[xg(t)] and then x = ¢(x—u(t), t), where g and u are to be determined. 
The Schrédinger equation transforms* into 


_h a (: ae ) a¢ 
2m age ‘oa m® 0g 


h? 


+4 pee +e (eu — Pie) -+( eu Fu ite; “)t¢ age, 
zm ot 


where & = x—u(t). Choosing now g and uw so that the coefficients of 06/0 
and &@ vanish, 


—mu = ihg, 
(3) ku— F—ihg = 0 or mi+ku = F, 
one is left with 
(4) —e jet Uke +3] 4 = ih =. 
6(t) = teu? — Fu—ihug—2-¢? 
2m 


= smu — ku’. 


Thence, € and ¢ variables being separable, 
7 s if . ( + S22 +) 
bd =, exp 7 [6(¢) +E, |dt ¢ exp(—gaé )IT, (aé), 


at = mk/h’, E, = (n+4)hwo (wo = [k/m]}), 


° . i? _ 
IT, = nth Hermite polynomial, N,2 = a/w?2"n! 


In words, a class of solutions of equation (2) is that formed from the stationary 
states of the unforced problem, with x replaced by & and with a phase 
exp(—1 h)({6dt+xg) ; by equation (3) & is nothing other than x—xo(t), and 
by (4) 6 is the classical Lagrangian for the unforced motion written as a 
function in time of classical forced position and velocity. 


*\n alternative route for transformation is via contact transformation of the classical 
Hamiltonian 














KERNER: QUANTUM MECHANICS 373 


Conventionally we may suppose an initial state ¥(x, 0) that coincides with 
one of the eigenfunctions v,,(x) of the unforced Hamiltonian. Then if the 
classical motion (1) be restricted to the initial conditions x9(0) = 0 = %o(0), 


V(x, t) is 
(5) Veet) = Nwexp( ; po(t)x exp —4 J (-+2,)d) 


X exp{ — 4a"[x—x0(t)]?} Hm (alx—xo(t)]), 


po denoting classical momentum mxXo(t). Plainly these eigenfunctions form a 
complete set at any time ¢-and a superposition of them can represent an 
arbitrarily specified wave function at any time. Though thy = H(t)W can 
possess no “‘true’’ stationary states, it is clear that in some suitable &(x, t)- 
space it may do so. 

The meaning of the é-stationary states is at once visible: #(f) is just x»(t): 
the probability density 


lal? = Nm exp| —o°[x—xo(t)]?} Hn? (alx—xo(t)]) 


dances a classical dance centered at the instantaneous classical position xo(¢), 
moving in toto classically. This generalizes the well-known “oscillating wave- 
packet” solution of the unforced problem. 

The computation of the transition probability 


Pin = 





%00 | 2 
\ Ym (x, t) tn® (x) ds | 3 


giving the probability that at time ¢ the oscillator is in the unperturbed state 
?, if initially it was in the state v», is readily performed using the generating 
function for the Hermite polynomials. The result is 


ML. ) v v- \ 12 
Pov = e a €y "hts *(€0) > 


where 
classical energy, 3mxo +3kx 0° 


€y(t) = 
quantum energy, hwo 


and » is the greater, uw the lesser of m, n; L*~* denotes the associated Laguerre 
polynomial. There follows the rigorous selection rule that the transition prob- 
ability vanishes only when the classically accumulated energy of the oscillator 
in units of hwo falls at a zero of the appropriate Laguerre polynomial. The 
discussion of several interesting illustrative examples will be taken up at a 


later time. 
DAMPING 


For the linearly damped motion of a particle in a field of force I(r) the 


Newtonian equation 
mf+yr = —vV(r) 
may be cast into Lagrangian form by means of 


L = exp(At)(4mr?— Vi(r)), \ = y/m. 
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A brief quantum-mechanical discussion of the damped oscillator was given 
first by Kanai (1948). Havas (1956, 1957) has considered quite generally the 
theory of multipliers, such as exp(At) above, that allow a Lagrangian formula- 
tion of a broad class of Newtonian problems not otherwise fitting into the 
Lagrange-Hamilton scheme of mechanics; and he has emphasized the im- 
possibility or else ambiguity (due to the multiplicity of possible integrating 
factors) of proceeding to quantization by the usual rules.* Although the 
meaning of doing so is not altogether clear at this point, we shall in the 
conventional way quantize the damped motion described by L and show 
that physically reasonable results are consequent. 

By the customary route, the Hamiltonian implied by Z is 


9 


(6) f= = exp(—At)+ V(r) exp(‘), Pp = mr exp(AL), 

and is neither the energy nor any constant of the motion, but is just the 
generator of motion. To place HZ in a more familiar light, make the contact 
transformation 


r = exp(—3A4)R, p = exp(3r2)P, 
giving the new Hamiltonian 


G= a + 3AR-P+exp(At) V(exp(— $d2)R), 


in which the whole burden of time dependence is in the potential-energy 
term. When V is homogeneous of the second degree, the V term is V(R), and 
G, being free of t, is a constant of the motion and is an energy though not 
the energy. For the one-dimensional oscillator, for example, 

- 
(7) Com +3XRP+4RR’. 


2m 


The further transformation in this case generated by the generating function 


F(R, 1) = — 4 10°+" Ro 
4mw w 
gives a new Hamiltonian 
) ian tie aman 
(8) P= 5 +4KX 


in the new momentum IT = (w/wo) P and conjugate coordinate Y = (A/2mawo) P 
+(wo/w)R, where wo is the undamped frequency (k/m)', w the damped fre- 
quency (K/m)? = (k/m—)?2 1)}, 

It is interesting to see here that the Gibbs statistical mechanics is applicable 
to ensembles of weakly coupled damped oscillators, just as for the undamped 
ones; for in the phase space (II;, X 4) Liouville’s theorem holds and I is con- 
served. Hence, for example, in a microcanonical ensemble mean values for 


*J am indebted to Dr. Havas for correspondence on this matter. 
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any quantity of interest may be found, in a parallel with the usual computa- 
tions; and in a canonical ensemble, with a “‘heat’’ bath of damped oscillators, 
there is a modulus of the distribution playing a role analogous to temperature 
and telling the preferred direction of flow of ! between ensembles with different 
moduli. There is, in short, a kind of thermodynamics of macroscopic systems 
damped in their microscopic coordinates. 

Turning now to quantization we have that, as a fundamental proposition 
stemming from the commutation rules, p may be represented as —ihv, and 
that, because only of the meaning of the Hamiltonian as generator of the 
motion, the Schrédinger wave function is controlled by thy = Hy, with H 
as in (6). That ¥ retains its meaning of positional probability amplitude is 
indicated by maintenance of probability conservation (since H is Hermitian) 
and of the validity of Ehrenfest’s theorem, to the effect that the quantal 
mean position is Newtonian. 

For the oscillator, either by direct transformation of the wave equation or 
by use of the transformed Hamiltonians (7) or (8), the states corresponding 
to the undamped ones come out to be 


Yn = exp (—! ct tn)r—(+ 4g") expe HT,,(8x exp(4M)), 
where the ¢«, are eigenvalues of T, ha(w+3), and Bt = mK /h? (K > 0). The 
factor exp({At) guarantees the time-independent normalization of y,. These 
are stationary states in R-space: a measurement of G for them is certain to 
give the value e,. So to speak the moving states are stationary with respect 
to the moving Hamiltonian G. The position—momentum uncertainty product 
is AvAp = h(wo/w)(n+3) while the mean energy (3mx*+$kx?) becomes 
hao(wo/w)(n+3) exp(—At). Altogether the quantum motion, like the classical, 
winds down (or up, for increasing ¢), |Y,|? shrinking eventually into a delta 
function. So long as quantum mechanics is linked to classical mechanics via 
Ehrenfest’s theorem such a result is inevitable; the ground state of the un- 
damped oscillator can not in the course of time be the favored one, as it 


might at first sight be expected to be. 


SIMULTANEOUS FORCING AND DAMPING 


The forced and damped oscillator, controlled classically by 
(9) mo+tyXotkxo = F(t), 


has the Hamiltonian 


») 


i r exp(—\t)+ (4kx°—xF) exp(dt) 


with Hamiltonian equations equivalent to (9). 
The Schrédinger equation becomes 


h° avy 1 2 . . oy 
= xp(—) + (1by*—xF) ex i ; 
dm OXPC WO ewe + (3kx°—xF) exp(rAt)y = th at 
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First put ¢ = x exp($M) and y = U(f, t); then with U = V exp[—(im/4h)¢], 
the transformed wave equation is 
I i ae 
—a— sa t [KE +4hA— F exp(Zry)cE]V = ih, 
ae agit l C+4ihy exp(3A4)¢]V ih 
with K as previously defined. From this point on, with the problem reduced 
to an equivalent forced and undamped oscillator problem in ¢-space, we pro- 
ceed as before, placing V = W exp[¢f(t)] and W = S(¢—wi(t), t) to give (» being 
¢—w(t)) 
has Te ae. 
He ait laka +A()]S = ih 
C 


2m dn- . 
provided f, w satisfy 
—mw = thf, 
(10) Kw-— F exp(3M) —ihf = 0 or mw+Kw = Fexp(4n), 
and A is 
A = Fith\+3mu?—} Kw". 
Finally, for K > 0, 


SS exp( 5 J(Até at) exp(— 38°n")H, (Bn). 


For K <0 (no oscillation at all classically in the damped but unforced 
motion) S is of continuum type. By (9) and (10) we identify 


w(t) = xo(t) exp(3At), 
A(t) = LotsmrAx 0X0 exp(At) +17hX, 


where Lo represents the classical Lagrangian for the damped but unforced 
motion as a function in time of the damped and forced position and velocity. 
Altogether, for the y, that are in correspondence with the v,, we get 


Y, = exp f Af Adt-tent-+imdx* exp(At) —x[pot3dAmxo exp) 
X exp{ —487[x—xo(t)]° exp(At) |, (8[x—xo(t)] exp(42)), 


where po is mxXoexp(At), the momentum conjugate to xo. The probability 
density again moves as a whole in the classical rhythm, simultaneously de- 
forming by a scale change exp(3AZ) in the x-scale. Asymptotically for ¢ > 7, 
lYn|? becomes a delta function centered at xo(t): for sufficiently long times 
the classical and quantal motions coincide. This does not, formally, contra- 
vene the uncertainty principle, for, as a calculation shows, the momentum 
becomes wholly indeterminate in this limit. On the other hand we recognize 
the gradual transition into the classical regime by the disappearance of h 
from the commutation rule 






x(mx) —(mx)x = th exp(—)t) 


connecting position and mechanical momentum mx. 
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Suppose the damped oscillator is hit by any transitory disturbing force F, 
F(t) ~ 0 (¢ > 7). The oscillator motion subsequent to time T decays classic- 
ally according to 

Xor = Xo(T) exp(—}A) cos wt. 


If the oscillator is charged with charge e, we may interpret —yx« as the 
approximate force of radiation damping, and may examine the Maxwell fields 
emanating from the quantal charge-current distribution (Schrédinger 1926). 
The charge density e|y,|? has a dipole moment ex»7; an observer in the radiation 
zone will see a dipole line centered at w with a width \; that is, he will see 
the classical result of a shifted and broadened line. Quite generally, in fact, 
the dipole field will be the same as that calculated classically. In the case 
that the oscillator is damped but not forced the quasi-stationary states y, give 
a nearly static charge distribution having no dipole or other odd moments. 
What radiation there is is confined mainly to the induction field; so the 
oscillator as it damps out creates variable fields predominantly in its immediate 
neighborhood, i.e. is surrounded by a ‘virtual’ field of varying energy 
content, one that can affect closely approaching particles but not distant 
ones. 
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AN INTERPRETATION OF THE LOW-LYING EXCITED STATES 
OF Mg” AND AI*! 


A. E. LITHERLAND, H. McManus, E. B. PAUwL,? 
D. A. BROMLEY, AND H. E. Gove 


ABSTRACT 


A description of the experimental results obtained for Mg?> and Al?> is given 
based upon the assumption that the excited nuclear states are rotational states. 
In Mg®® and Al?® members of four rotational bands can be identified. The band 


based on the ground state can be assigned AK = 5/2 whilst the three excited- 
state bands, two positive parity and one negative parity, can be assigned K = 
1/2. The dipole y-ray transitions between the K = 1/2 bands and the ground state 


K = 5/2 band are forbidden on the rotational model and it is a salient feature 
of the y-ray cascading in Mg** and Al® that the ground-state transitions are 
always amongst the weakest transitions from excited states assigned to K = 1/2 
bands. Besides giving a quantitative account of the y-ray decay of the excited 
states the rotational model also predicts the number and type of the bands 
observed. The experimental values of the Mg?(d, p)Mg®® deuteron stripping 
reduced widths can also be approximately predicted by the model. 


I. INTRODUCTION 

The low-lying excited states of the nuclei Mg*® and Al*> have recently been 
studied extensively and a large body of information gathered. This information 
comprises evidence from proton and neutron capture y-ray studies (Ager- 
Hanssen et al. 1956; Campion and Bartholomew 1957; Bromley et al. 1957; 
Craig 1956; Endt and Braams 1957; Green et al. 1955; Litherland et al. 1956; 
Varma and Jack 1956); deuteron stripping and proton scattering (Endt and 
Braams 1957); inelastic proton scattering by Gove et a/. (1956-57); and beta- 
decay measurements by Maeder and Stahelin (1955) and Gove ef a/. (1956-57). 

The experiments on Meg?*!(p, y)Al® described by Litherland et a/. (1956) 
first suggested that the excited states of Al** could be described in terms of 
overlapping bands of states having some of the characteristics of rotational 
states and a brief description of the states of Al*® as rotational states was 
given. It is the object of Section II of this paper to apply the rotational 
description to the states of Mg*> and Al** in more detail and of Section III 
to make a comparison of some of the properties of Mg*> and Al*® with those 
calculated from the eigenfunctions, for a nucleon in a spheroidal potential, 
determined by Nilsson (1955). 

The success of the description for Mg*> and Al’** has resulted in its extension 
to neighboring nuclei. Attempts to do this have been made recently by Paul 
(1956) for F!°, Sheline (1956-57) for Al?8, and Bromley et a/. (1957) for S?’°. 
Though in each case there is less evidence than in the case of Mg®® and Al’ 
the success obtained seems to justify further experiments to test the inter- 
pretations. 
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II. ROTATIONAL DESCRIPTION OF THE STATES IN Mg®® AND AI? 


1. The Sequences of States 
A summary of the experimental results on the y-ray cascading in Al** is 
shown in Fig. 1* and a separation of the low-lying states of Al®* into bands 


of rotational states is shown in Fig. 2. The evidence for this separation has 
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Fic. 1. A summary of the decay scheme showing all states in Al?* which were studied and 
their spins and parities where they could be established. The principal gamma transitions from 
each state are shown with the measured branching ratios. The energy of the excited state is 
given in Mev. 
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The decomposition of the low-lying states in Al?> into bands of rotational states 
with the parameters used to fit the bands. Equations 1, 2, and 3 given in the text are used. 
The values of B marked with an asterisk show that the same vibrational parameters were used 
as for the sequence of states in Mg**. 


*Where the energies shown in Fig. 1 differ from those on subsequent figures the energies in 
Fig. 1 are to be regarded as more accurate values. 
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been discussed by Litherland e¢ a/. (1956) and consists mainly of total angular 
momentum and de-excitation y-ray branching-ratio determinations. The 
states shown as dashed lines are predicted on the basis of the listed parameters 
and Equation 2. Evidence for the predicted 7/2+ state near 2.80 Mev. has 
recently been obtained by Bromley et al. (1957). Fig. 3 shows a summary of 
the present knowledge of the excited states of the mirror pair Mg? and 
Al®?®, Some of the details of this figure have been discussed recently by Gove 
and Litherland (1957). 
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Fic. 3... A comparison of the mirror pair of nuclei Mg?> and Al®. 


In describing the states of Al** as rotational states, the concepts developed 
by Bohr (1954) and Bohr and Mottelson (19538) are used. They considered 
the model of a nucleus of stable spheroidal shape capable of executing collective 
rotations about an axis perpendicular to the axis of symmetry.* It has been 
shown by Bohr (1954) that such rotations, for the case of an even-even 
nucleus, have an energy spectrum given by 


(1) E(J) = A[J(J+1)], 


where J is the total angular-momentum quantum number of the nucleus and 


*The nuclei in the region of A = 25 may be considered to be prolate spheroids since the 
electric quadrupole moments of Na® and Al?’ are positive. 
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A is constant for each band and is equal to h?/2/; J is the effective moment 
of inertia of the nucleus undergoing collective rotations. 

For some time Mg** has been known to exhibit the 0+, 2+, 4+ sequence 
of spins and also to have approximately the correct energy-state spacing given 
by Equation 1. The small discrepancy may be attributed to a vibration- 
rotation interaction term (cf. Fig. 2). Bohr and Mottelson (1955) also give 
the expected energy spectrum for an even-odd nucleus: 


(2) E(J) = A[J(J+1) +6x¢.1,2a(— 1)7t+!?(J+1/2)] 
+ BIJ (J+1)+6x,12a(—1)7+!/?(J+1/2)]?. 


ig. 4 shows a plot of Equation 2 as a function of a for various values of J 
with A = 1/2 and B = 0. Using Fig. 4+ the order and relative spacing of the 


STATE ORDER IN K=!4 SERIES 2 
AS FUNCTION OF PARAMETER @ 


E (u)=A y(utntal-1y?*2 (y+) 






LOWEST STATE 






Fic. 4. The order of the excited states for a rotational band with Q = 1/2 is displayed as a 


function of the parameter a. For a = 0 the order of the states is given by J(J+1). For a = 
— 3.5 the order of the excited states is 3/2, 7/2, 1/2, which is the order observed for the negative- 
parity band in Al*, 


rotational states can be obtained for various values of a. In Equation 2, B is 
the vibration-rotation interaction parameter. @ is a parameter which is 
relevant only to the A = 1/2 bands and is a result of the assumption of 
reflection symmetry of the nuclear shape and the coupling of the particle 
motion to the rotation. a, which is often referred to as the decoupling para- 
meter, is given by Equation 8 (Bohr 1954): 


(3) a= >> (-1)*""(j+1/2)C,’. 
J 


In Equation 3, 7 is the total angular-momentum quantum number of the odd 
nucleon in a spheroidal nucleus. K is the projection of the total nuclear angular 
momentum J onto the axis of symmetry of the spheroid. The parity of the 
wave function is a good quantum number, but, because of the absence of 
spherical symmetry, 7 is not a good quantum number. However the projection 
of j along the axis of symmetry of the spheroidal potential is a good quantum 
number (Bohr and Mottelson 1955). This last quantum number is usually 
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called Q. The single-particle wave function in a spheroidal potential x9 can 
then be expressed in terms of a sum of single-particle wave functions for a 
spherical potential y,;: 


(4) xo°-= = Ci;. 


The C, coefficients are normalized so that }>; C;* = 1 and the square of the 
coefficient C,; is then the probability that the single nucleon has total angular- 
momentum quantum number j. 

The system of a spheroidal nucleus plus a single nucleon has a good quantum 
number J together with M, the projection of J upon some axis fixed in space, 
and K, the projection of J upon the axis of symmetry of the spheroid. A is 
equal to 2 in the absence of beta and gamma vibrations, which is usually 
the case for low-lying states. The wave function for such a system has been 
given by Bohr and Mottelson (1953). 

K and Q are, however, good quantum numbers only in the absence of what 
has been called ‘rotation-particle coupling’ or RPC (Kerman 1956). The 
rotations of the spheroidal nucleus perturb the nucleon motion very little 
provided the frequency of rotation of the nucleus is much smaller than the 
frequency associated with the nucleon motion. The latter frequency is related 

to the excitation energy of the single nucleon, and the former frequency is 
related to the energy separation of the rotational states. Thus if the rotational 
bands due to different single-nucleon excitations overlap, then perturbations 
of the bands are to be expected. However, the RPC in first order couples 
together only bands of the same parity and whose A’s differ by one or zero 
(Kerman 1956). It is probably due to this that Mg?> and Al*® show rotational 
bands so clearly amongst their low-lying states since the two principal strongly 
overlapping bands have AK = 2. Other nearby nuclei such as Na?*, which 
might be expected to show rotational spectra, have neighboring bands which 
differ in K by only one unit. Consequently Na?* would be expected to show 
perturbed bands. This point is discussed briefly in Appendix A. 

In the case of the K = 1/2 bands in Mg* and AI*, the effect of the RPC 
is to modify the value of the parameters 4 and a obtained from the bands. 
Until further members of the bands are located the values shown in Fig. 2 
should be treated with some caution, though, as will be discussed later, the 


mutual perturbing effect of the two positive parity K = 1/2 bands is probably 
small. 
The decoupling parameters a obtained from the bands with A = 1/2 are 


intimately related to the state of the odd nucleon in the spheroidal nucleus 
(cf. Equation 3). The case of the anomalous state order of the negative-parity 
AK = 1/2 band in Al*® (Fig. 2) is particularly interesting since it has a value 
of a = —3.16 which is close to the limiting value of a = —4.00 obtained by 
making C7). = 1; Fig. 4 shows the order of the states for these values of the 
decoupling parameter. This suggests that the band can be described fairly 
well by an f7,2 particle in a rotating spheroidal nucleus. Moszkowski (1955) 
has shown that for a nucleus in the vicinity of A = 25 with the shape of a 
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prolate spheroid there would be a low-lying single-nucleon state with Q = 1/2 
which originated from the 1f7,/2 shell. The deviation from spherical symmetry 
of the nucleus results in the If7;2 state being split into states with Q = 1/2, 
3/2, 5/2, and 7/2 with the Q = 1/2 state lying lowest for a prolate nucleus. 
The anomalous order of the low negative-parity states of Al*® can thus be 
well described by the rotational equation 2. From Equation 3 alone no con- 
clusions can be reached concerning the origin of the even-parity A = 1/2 
bands other than that various mixtures of 7 = 1/2, 3/2, 5/2 are allowed by 
the observed values of a. The evaluation by Nilsson (1955) of the single- 
particle wave functions for a spheroidal potential permits the calculation to 
be made of the values of a. This will be discussed later. 


2. The Moments of Inertia Obtained from the Bands 

The moments of inertia of the bands shown in Fig. 2 show that there is 
an appreciable difference between the various bands, with the moments of 
inertia for the K = 1/2 bands being significantly higher than the values for 
the K = 0 and K = 5/2 bands. The difference in moment of inertia between 
even-even nuclei and even—odd nuclei has been discussed in some detail by 
Bohr and Mottelson (1955) and by Kerman (1956). 

The values of the moments of inertia are approximately equal to the values 
for the rotation of a spherical rigid nucleus. For a spherical nucleus with 
mass number 25 and Ro = 4.1X10-' cm., the parameter A = h?/27 shown 
in Fig. 2 is 128 kev., which is quite close to the values for the K = 1/2 bands. 
If the nucleus is considered to be a spheroid of constant density, then for 
mass number 25 the moment of inertia for rigid rotation about one of the 
minor axes is given to first order by the following equation: 


(5) h?/2T neta = 128(1 — 0.36) kev. 


If the nuclear distortion parameter 6, which is approximately equal to AR/R» 
where Ro is the mean nuclear radius and AR is the difference between the semi- 
major and semi-minor axes of the spheroidal nucleus, is obtained from the 
quadrupole moment of the neighboring nucleus Al®’, the value of h?/2J ngia 
is changed very little from the value for a spherical nucleus since 6 ~ 0.3. 
The situation for mass number 25 can be compared with the well-studied 
rotational spectra in heavie- elements. There the moments of inertia deduced 
from the spacing of the nuclear states lie roughly between a third and a 
quarter of the values for rigid rotation. The significance of this fact has been 


discussed at length by Bohr and Mottelson (1955). 


3. Particle Widths in Mg*> and Al*® 

The C;, coefficients in the case of Mg®® and Al** have a readily appreciated 
physical significance. The square of a C, coefficient is proportional to the 
nucleon reduced width for the addition of a nucleon of angular momentum 7 
to the ground state of Mg*4 as shown in Equation 7. 

The general expression for the stripping reduced width for rotational states 
has been given by Yoshida (1954) and by Satchler (1955) and a modified 
form of these expressions is given by Equation 6: 
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» 227, r £ 

(6) 9° = ee @)? DIC TGK AK. 
oJ 4 j 
In Equation 6, (J;,K,) and (J;, K,) are the angular momenta and their 
projections along the axis of symmetry of the nucleus for the initial and final 
states, while (jQ) are the corresponding quantities for the added nucleon. The 
quantity (09)? is dimensionless and expresses how well the nuclear wave 
function for A+1 particles can be split up into the product of a wave function 
for A particles and the wave function for a single nucleon. (4)? is assumed 
to be unity (independent-particle model) for the purposes of this paper. 
Equation 6 reduces to Equation 7 when J; = 0: 


= 2 é ! 2 

(7) iy “zs ICy,|". 

From this equation a sum rule for reduced widths obtained from stripping 
reactions and elastic proton scattering on a zero-spin nucleus follows: 


2 
(8) > ic;,1? ate re. er) ; 
Jf Jf 4 
In equation 8 the sum is taken over all states belonging to one particular 
band. 

In order to compare the nucleon reduced widths with theory, calculations 
of the C, coefficients must be made using some model. Section III.3 of this 
paper contains such a comparison. It is, however, possible to use Equation 8 
without making detailed calculations once the rotational bands have been 
identified. Table | shows some examples of the sum rule. The neutron capture 
probabilities A quoted by Holt and Marsham (1953) and by Hinds et al. 


TABLE | 
EXAMPLES OF THE REDUCED WIDTH SUM RULE 











(1) (2) (3) 
5 15.4 5/2+ 
9 - 17.3 1/2+ 
ui Mg A. 5 1/2+ 
14 >39 1/2— 
9) 14.7 1/2+ 
8} Siz? 16.9 3/2+ 
104 19 7/2— 





(1) Nilsson orbit number. 

(2) The sum of the A in arbitrary units for the 
various bands. 

(3) The Q value and parity of the band. 


(private communication) are proportional to 6,,?(2J;+1). Apart from the 
high value for orbit number 14,* there is good agreement for the nuclei Mg?® 


*The high value for the summation in the case of the negative-parity band is probably due to 
the use of the peak cross section of the angular distribution as a measure of the reduced strip- 
ping width. Any contribution from compound-nucleus formation would tend to make the 
reduced width smaller. Since the / = 3 angular distribution was measured only from 0° to 
60°, this effect cannot be estimated. 
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and Si**,.* The possibility of rotational states in Si?’ has been discussed recently 
by Bromley et al. (1957), who make the orbit number identifications. An 
application of Equation 6 to the higher members of the rotational bands is 
given in Appendix C. In Appendix D the application of Equation 6 to some 
stripping reduced widths in the d-shell is discussed briefly. 


4. Gamma-Ray Transition Widths 

The y-ray decay within a rotational band of states is characterized by the 
enhancement of the £2 radiation and the subsequent effective competition 
with the M1 radiation also present. Enhanced E2 transition widths have 
now been measured for a large number of medium and heavy nuclei and 
their association with rotational groups of states established in many cases 
(Alder et a/. 1956). In the case of Mg?* and Al*> there are no lifetime measure- 
ments of the rotational states and there is no evidence from Coulomb excitation 
so that the enhancement of the £2 transitions within a rotational band cannot 
be established. 

In the case of Mg** however electron scattering measurements by Helm 
(1956) have permitted an estimate of the lifetime of the first excited state 
of Mg** to be made. The results of Helm (1956) for Mg**, Si*’, S%* together 
with those of Devons et al. (1956) for Ne?®, for the lifetimes and ratios of 
the measured lifetime to the Weisskopf £2 unit (Wilkinson 1956), are shown 
in Table II. In all cases there is enhancement of the £2 lifetime over the 
single-particle estimate. Using the measured lifetime for Mg** it is possible 


TABLE II 
E2 TRANSITION WIDTHS FOR SOME EVEN-EVEN NUCLEI WITH 
20 < A < 32 


(1) (2) (3) 








Ne 7.643.3 0.82 11.15 
Mg™ 19 0.34 10.8 
Sis 6 1.1 6.2 
3s” 1.6 4.0 5.7 


(1) Lifetime of the first excited state in units of 10-43 seconds. 

(2) Transition width in millielectron volts. 

(3) Ratio of experimental transition width to the Weisskopf 
estimate (Wilkinson 1956). 


to calculate the intrinsic quadrupole moment of Mg*‘. The value of Qo = 0.43 
barns is very close to the value obtained by Gugelot and Rickey (1956) who 
deduced a value of Qo = 0.45 barns from the elastic scattering of 40 Mev. 
a-particles from Mg‘. It is also close to the values obtainable, assuming 
strong coupling, from the measured quadrupole moments of Na** and Al?? 
of 0.55 and 0.41 barns respectively. 

It seems likely from some of the properties of neighboring nuclei that £2 








* Note added in proof.—The recent studies on the Mg**( p, d) Mg*4 reaction by Bennett (Bull. 
Am. Phys. Soc., Ser. 2, 26, 1958) are relevant here. The reduced neutron widths for the for- 
mation of the first three states in Mg* have been interpreted successfully by Blair (Bull. 
Am. Phys. Soc., Ser. 2, 26, 1958) within the collective model framework. 
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transitions in Mg®®> and Al*® would also be enhanced. Though at present this 
cannot be demonstrated from the absolute transition widths the ratios of the 
reduced M1 to reduced £2 transition widths when expressed in terms of 
Weisskopf units show that either the £2 transition widths are enhanced or 
the M1 transition widths are inhibited. This point is illustrated in Fig. 5 for 
the two identified bands exhibiting 2 crossover and cascade transitions. 
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Fic. 5. The experimental 1/1-E2 ratios within the two lowest bands in Al?> are compared 
with the ratios of the corresponding Weisskopf (W) units. The 1/1 transition widths have been 
taken as unity for purposes of this comparison. 











The prominent £2 crossover transitions to the first excited state from the 
fourth excited state of Mg*> and Al*> observed by Litherland et al. (1956) 
and Gove et al. (1956) were pieces of evidence which led to the rotational 
interpretation of the excited states of Mg?® and AI’, 

The rotational model makes definite predictions concerning the ratios of 
the reduced 2 transition widths (Alaga et al. 1955). There are two cases 
in Al*?* where such a comparison can be made. These two cases are shown 
in Fig. 6 where the results of Litherland et a/. (1956) have been used. Although 
the measured ratios with their errors include the theoretical ratios it is clear 
that more accurate measurements would be worth making. 

For dipole radiation between rotational bands the quantum number A 
must change by 0 or +1. If K changes by more than one unit the transition 
is said to be K-forbidden (Alaga et a/. 1955). The operation of the A-selection 
rule in Al*® was another piece of evidence that led to the rotational interpre- 
tation of the states. It was found by Litherland ef al. (1956) that transitions 
to the ground state of Al?> were in general weak compared with the cascading 
radiations even though nearly all the radiations were dipole in nature.* Fig. 7 














a ‘T he strong inhibition of the ground-state tri insitions from the states at 4.22 Mev. and 4.60 
Mev. is striking. However, it is difficult to attribute this inhibition to a K-selection rule since 
the most likely K value for each state would be 3/2. A possible explanation of the inhibition is 
discussed in Appendix B. 
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Fic. 6. The experimental ratios of the reduced E2 transition widths within the two lowest 
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Fic. 7. Branching ratios in the decay scheme of the state in Al?® at 3.72 Mev. with the 
reduced transition probabilities for the electric dipole and magnetic dipole transitions shown 
in such a manner that K-forbidden transitions have the value unity. 


shows the decay scheme of the 7/2—, K = 1/2 state to the lower states in 
Al*® together with the reduced y-ray transition widths. The operation of the 
K-selection rule is demonstrated quite clearly. The transitions in Mg*® and 
Al®?®> show many examples of the K-selection rule, and these are listed in 
Table III. The numbers in the table are the ratios of the reduced dipole 
transition widths between the excited states, listed horizontally and vertically, 
to the reduced transition width for the ground-state dipole transition. The 
energies are in Mev. The branching ratio of the 3.40 Mev. state in Mg?® is 
from the Mg*4(, y) reaction measured by Campion and Bartholomew (1957). 
Quadrupole mixing has not been allowed for, since in most measured cases 
it has been shown to be small (Litherland et al. 1956). The values shown in 
parentheses are preliminary values (Bromley et a/. 1957). 

Violation of the selection rule is due to K not being a perfectly good quantum 
number, and this is a partial reflection of the presence of a vibration-rotation 
interaction and of the admixture of the bands with another band or bands 
which differ in their K by 0 or +1. Such a perturbing band would be the 
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TABLE III 
EXAMPLES OF K-FORBIDDEN TRANSITIONS IN Al®®> AND Mg?5 
Al 
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expected, but hitherto unobserved, band with AK = 3/2 based upon the d3,2 
single-nucleon orbit. The degree of violation of the K-selection rule varies 
considerably and this may be partly due to the variations of the magnitudes 
of the allowed dipole reduced widths among themselves, as discussed below. 
It is possible using the equations given by Alaga et al. (1955) to give some 
account of the variations of the ratios of the reduced-dipole transition widths. 
For transitions from members of one band of rotational states to another 
band it is possible to form ratios which can be compared with theory. There 
is only one parameter 6 available to fit several ratios. Fig. 8 shows the experi- 
mental ratios compared with the theoretical ratios for the M1 transitions 
between the band based on the 2.51 Mev. state in Al*® and the band based 
on the first excited state of Al*®® at 0.45 Mev. Fig. 9 shows the experimental 
ratios for the £1 transitions between the negative-parity band and the band 
based on the first excited state at 0.45 Mev. Although the fit to the experi- 
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each state is taken to be unity for this comparison. 
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mental ratios is not perfect there is enough qualitative agreement to justify 
the hope that if all the perturbing effects, discussed briefly in Section III.4 
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are compared with the theoretical ratios for bg, = +0.3. The strongest El transition width 
from each state is taken to be unity for this comparison. 


of this paper, are taken into account the agreement would be improved. The 
large variations in the ratios also illustrate one of the reasons why Table III 
of K-forbidden transitions contains such apparent large variations in the 
degree of K-forbiddenness. 


5. Beta Decay of Al*® and Na*® 

The beta decay of the nuclei Al*® and Na®®> has been studied, and several 
points of interest emerge. In the positron decay of Al*® there is, besides a 
transition to the ground state of Mg” with a log ft value of 3.5, a transition 
to the state at 1.61 Mev. with a log ft value lying between 4.7 and 5.6 
(Maeder and Stahelin 1955). Possible values of the spin and parity of the 
1.61 Mev. state in Mg*> and Al*> have been discussed by Gove et al. (1956- 
1957) and by Litherland et a/. (1956) and the conclusion reached that an 
assignment of 7/2+, with a K of 5/2, is the most reasonable. This is in 
agreement with the allowed log ft value. As the ground state of Al*® is 
5/2+, transitions to the state in Mg*®® at 0.98 Mev. (3/2+) would also be 
allowed (AJ = 1), as well as those to the 1.96 Mev. state (5/2+) (AJ = 0). 
These transitions have not been observed and in the case of the transition 
to the 0.98 Mev. state in Mg” there is a lower limit to the log ft value of 4.9 
obtained by Stribel (1956). The measurements, though consistent with the 
K-forbiddenness of the transition to the 0.98 state, are clearly inadequate to 
vield any quantitative idea of the degree of K-forbiddenness. 

An accurate measurement of the positron branching to the 1.61 Mev. state 
in Al*> is of importance because the log ft value can be calculated if the 
rotational character of the states is assumed (Bohr and Mottelson 1953; 
Nilsson 1955). The value of log ft calculated is 4.24, which is smaller than the 
experimental lower limit. The uncertainty of the experimental log ft. value 
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is almost entirely due to the difficulty of estimating the contribution of Al** 
to the positron annihilation quanta in the decay of Al*> studied by Maeder 
and Stahelin (1955). 

The beta decay of Na”® is quite different from that of Al®> as transitions 
with very similar ft values are observed to the ground state (5/2+), second 
excited state (3/2+), and third excited state (7/2+) (Gove et al. 1956-57). 
As the ground state and third excited state have K = 5/2 and the second 
excited state has K = 1/2, the most likely spin, parity, and AK quantum 
number for Na®> are (5/2+ 3/2), which is anomalous, since one might expect 
the ground-state properties of Na?> to resemble those of Na?*. However, as 
discussed in Appendix A, the states of Na**, if they can be described as 
rotational states, are anomalous in that respect and require the introduction 
of RPC to describe their order. It is possible that in the case of Na®> the 
5/2+ first excited state has crossed over the 3/2+, owing to the RPC, to 
become the ground state (cf. Fig. 19). This conjecture implies that the state 
would be of mixed K = 3/2, 5/2, which would make transitions to the states 
of K = 1/2 allowed by the K-selection rule. A low-lying state of J = 3/2+ 
would then be expected for Na?®, 

The foregoing discussion results in the conclusion that a satisfactory 
qualitative description of the spectra of Mg*® and Al*> can be obtained by 
the simple extension of the model of a single nucleon in a spherical potential 
to the model of a single nucleon strongly coupled to a rotating spheroidal 
potential. The eigenvalues and eigenfunctions for a nucleon in a fixed 
spheroidal potential have been calculated by Nilsson (1955), and the applica- 
tion of these calculations to Mg*> and Al*® will now be discussed. 





III. COMPARISON OF THE EXPERIMENTAL DATA FOR Mg*®® AND Al*® WITH 
THE CALCULATIONS OF NILSSON 

1. The Energy States in a Spheroidal Well 

The calculation of the energy states in a static spheroidal potential has 
been carried out by Moszkowski (1955), Nilsson (1955), and Gottfried 
(1956). The calculations of Nilsson are most suitable for our purpose as they 
include the calculations of the eigenfunctions for the region of A = 25. Fig. 10 
shows the energy states of the model in the vicinity of A = 25 obtained by 
extending Nilsson’s calculations of the energy of an individual particle in a 
spheroidal potential to larger values of 7. For this purpose Nilsson’s equation 
(C3) has been rewritten as follows: 


E 
(9) = 
3) 4.Rado 
where the following relationships between e, 7, and 6 were employed: 


(10) Kn = e[1 —1/3.€2—2/27.€3]-'8 
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Fic. 10. The energy of states in a spheroidal potential taken in part from the paper ‘‘Bind- 
ing states of individual nucleons in strongly deformed nuclei’’ by Nilsson (1955). The orbit 
numbers are those given by Nilsson. Orbits 5, 9, 11, 14 are those identified in Mg?® and Al?®. 


The value x, which measures the strength of the spin-orbit coupling, used 
in calculating the curves for Fig. 10 was 0.08. This appears to be a better value 
for light nuclei than 0.05 and, in this connection it is probably worth noting 
that the d5,2-d3,2 splitting of 5.08 Mev. in O'? would require «x = 0.13. In 
view of the uncertainty of the value of « to use for A = 25, all subsequent 
graphs are given as functions of 7 since it is then straightforward to consider 
the variation of the strength of the spin-orbit interaction. 

The value of 7 used in the subsequent comparisons varies from approxi- 
mately +3 to +5. There is no reason a priori for expecting 7 to be a constant 
for each rotational band and the fact that different values of y give better 
agreement between theory and experiment in the subsequent comparisons 
probably implies that neglected effects, such as RPC, are significant. 


2. Sequences of States for Mg*® and Al*® 

Two protons and two neutrons are allowed for each orbit, so that, con- 
sidering Mg*4 to be a prolate spheroid with 0 < 7 < 6 as indicated by the 
magnitude and sign of the neighboring quadrupole moments, orbit numbers 
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1, 2, 3, 4, 6, and 7 (see Fig. 10) are filled. The value of 6 depends upon the 
strength of the spin-orbit interaction through Equations 10 and 11. 

The next four orbits available are numbers 5, 9, 11, and 8, with number 14 
being the lowest negative-parity orbit. Of these orbits all but the Q = 3/2, 
number 8 orbit can be identified amongst the low-lying states of Mg*® and 
Al (cf. Fig. 2). Orbit number 5 can be identified with the ground-state 
band (Q = 5/2) and orbit number 9 with the Q = 1/2 band based on the 
first excited state at 0.58 Mev. in Mg* and 0.45 Mev. in Al*>, The next 
two 2 = 1/2 orbits identified in Mg** and Al*> can be identified with orbits 
number 11 and 14. The latter is based on the If7,2 orbit in a spherical potential. 
There is no clear indication at present for the band based on orbit number 8 
(Q = 3/2). However, the region above 4 Mev. excitation where the lowest 
member of this band might be located is not known quite so well as the region 
below. In the case of Mg*®, magnetic analysis of charged-particle reactions 
indicates a number of unidentified states in the region of 4 Mev. excitation 
(Endt and Braams 1957), and in the case of Al** the region above 4 Mev. 
excitation is complicated by the onset of inelastic scattering in the Mg*4+ p 
reaction. It is therefore possible that the Q = 3/2 band might be found if 
the less well known regions of Mg*® and Al**® were searched.* 

The order of the bands number 5 and 9 in Mg*® and Al** depends sensitively 
on the prolateness of the spheroid, and the close proximity of the first Q = 1/2 
band to the Q = 5/2 ground-state band suggests that the parameter 7 in the 
case of Mg* and Al* is approximately +2.0. However if the individual- 
particle energies of Fig. 10 are summed for the 25 particles the potential- 
energy curves, as a function of , can be calculated and the equilibrium 
distortion estimated. The results of the calculation are shown in Fig. 11. On 
the left-hand side of the figure the normal configurations are plotted in which 
24 particles fill the orbits number 1, 2, 3, 4, 6, and 7 (Nilsson’s notation 
shown in Fig. 10) and the odd particle is placed in orbit 5, 8, 9, 11, or 14. In 
the right half of Fig. 11, potential plots for states of core excitation are shown; 
in each case one of the four particles is removed from orbit 7 and paired with 
the odd particle in the upper state. No contribution for differences in pairing 
energy has been included. The results indicate an equilibrium value of » ~ 4.5 
but do not reproduce the sequence of the observed bands in Al*>.t However 
the potential-energy plots do indicate that orbit 14 is expected to play an 
important role among the low states of Al**, a result which is not immediately 
apparent from Fig. 10. Also the calculations show that the so-called ‘excited- 
core’ states can be expected among the higher excited states; their exact 
position being in doubt because of pairing-energy uncertainties. A tentative 
identification of two members of the band based upon the lowest of the 


*If the 3/2+ state in Al®5 at 4.22 Mev. excitation were the first state of the K = 3/2 band, 
then the low reduced proton width of this state found by Mooring et a/. (1951) would be in 
disagreement with the rotational interpretation of the states of Al®5. The reduced elastic proton 
width of the 3/2+, K = 3/2 state can be calculated from Fig. 13(e) and Equation 6, assuming 
no K admixtures and no inelastic scattering, and is approximately the single-particle width. 
A possible explanation for the low reduced proton width is given in Appendix B. 

{The value of » = +4.5 corresponds to a distortion parameter 6 = +0.23 for x = 0.05 and 
6 = +0.36 for « = 0.08. In Section 3 it was pointed out that for nuclei in the region Na® to 
AI? the electric quadrupole moments indicate 6 ~ 0.3. 
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excited-core configurations of Fig. 11 is given in Appendix B. While these 
calculations are not quite consistent among themselves it is noteworthy that 
this very simple model is in more than qualitative agreement with this very 


complicated nucleus. 
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Fic. 11. The potential-energy plots obtained from Nilsson (1955) for Mg®® and Al? are 
shown for several of the possible configurations. 
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Having assigned orbit numbers to the various bands, it is a straightforward 
matter to calculate the values of the parameters a. Fig. 12 shows the results 
as a function of 7 for the orbit numbers 9, 11, and 14. The experimental values 
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Fic. 12. The decoupling parameter a is shown as a function of the parameter 7 for the 
Nilsson orbits numbers 9, 11, and 14. 
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of a, obtained from the observed order of the energy states, are also shown. 
These comparisons indicate values of 7 in the vicinity of +4, which is approxi- 
mately the same as the value obtained from Fig. 11. Until other members 
of the K = 1/2 bands are identified it would be premature to speculate upon 
the indicated differences of 7 from band to band. 


3. Particle Widths in Al** and Mg*® 

It is also a straightforward matter to calculate the values of C/ for each 
orbit, and these are presented graphically in Fig. 13. The graphs show some 
common pronounced characteristics. At 7 = 0 one of the C;?’s becomes 


10o- 





Fic. 13. Graphs of C,? versus 7. 
(a) Nilsson orbit number 6 based on the 1d5/2 shell. Q = 1/2. 
(b) Orbit number 7 based on the 1Id5,2 shell. Q = 3/2. 
(¢c) Orbit number 9 based on the 2s,/2 shell. Q = 1/2. 
(d) Orbit number 11 based on the Id3,2 shell. Q = 1/2. 
(e) Orbit number 8 based on the 1d3,2 shell. Q = 3/2. 
(f) Orbit number 14 based on the 1f7;2 shell. Q = 1/2. 
The positive and negative signs give the signs of the C;'s. 


unity. This C,? is the one for which the j value is equal to that of the parent 
orbit in a spherical potential. As the prolateness or oblateness of the spheroidal 
potential becomes more pronounced, coefficients with other values of 7 in- 
crease whilst the coefficient whose value was originally unity decreases. The 
relation between the stripping reduced width and the C,? is given by Equation 7. 
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Table IV shows a comparison of the relative reduced widths for the Mg*4 
(d, p)Mg*> reaction compared to the values of C,;* for 7 = 3; Table V shows 
the comparison of the reduced proton-scattering widths* for Mg?*(p, p) obtained 
by Koester (1952) with the C, for » = +3 and +6; the comparisons in 
Tables IV and V show that, although the agreement is fairly good, there 
are some obvious discrepancies. 


TABLE IV 


THE RELATIVE REDUCED STRIPPING WIDTHS FOR 
Mg*4(d, p)Mg?® ARE COMPARED WITH THE VALUES OF 
THE COEFFICIENTS |C;? 





(1) (2) (3) (4) (5) 
0 2 5/2+ 15.4 1 
(0.582 0 1/2+ 5.8 0.44 
0.976 2 (3/2+) 6.0 0.36 
1.612 (7/2+) — 0 
1.957 2 (5/2+) 5.5 0.20 
2.565 0 1/2+ 3.2 0.35 
2.742 (7/2+) 0 
2.806 2 (3/2+) 10.3 0.64 
3.405 1 (3/2—) 11.2 0.16 
3.899 “is (5/2+) ? 0.01 
3.972 3 (7/2—) 25 0.82 


(1) Energy of excited state in Mev. (Endt and Braams 
1957). 

(2) /-value of stripping pattern. 

(3) Spin and parity of excited state. The values in 
parentheses have been obtained by comparison with Al, 
the mirror nucleus. 

(4) A from Holt and Marsham (1953) and Hinds et al. 

(5) |C;|? for » = +3 obtained from Nilsson (1955). 


TABLE V 


‘THE REDUCED WIDTHS FOR THE SCATTERING OF PROTONS 
BY Mg?! ARE COMPARED WITH THE VALUES OF |C;!? 








(1) (2) (3) (4) (5) 


3.09 


3/2— 0.54 0.16 0.30 

3.72 7/2— 0.99 0.82 0.60 
3.85 1/2— 0.21 0.015 0.06 
3.88 5/2+ 0.006 0.01 0.30 
? 5/2— ? 0.005 0.04 


(1) Energy of excited state of Al*®. 

(2) Spin and parity (Koester 1952; Litherland et a’. 
1956). 

(3) Ratio of the reduced width to h?/ya multiplied by 
3(2/+1), where uw is the reduced mass and a is the inter- 
action radius assumed: a = 5.05X10°" cm. (Koester 
1952). 

(4) |C;|? for n 

(5) |C,|? for 


+3. 
+6. 


The agreement between theory and experiment in Tables IV and V can 
be improved by varying the parameter 7 independently for each band. Also 
the agreement could possibly be improved by varying the parameters assumed 


*The reduced width for the 7/2— state is smaller by a factor 3/4 than the value quoted by 
Koester since he assumed 5/2— when calculating the reduced width. 
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by Nilsson for the spheroidal potential and by introducing RPC. Such 
attempts to improve the agreement would scarcely be justified in view of 
the difficulty in obtaining the reduced widths from the stripping angular 
distributions. The neutron capture probabilities A have been used because 
many of the results from the (d, p) and (d, ”) reactions in the region A = 25 
were analyzed using the Born approximation theory of Bhatia et al. (1952). 
If the theory of Butler (1951) is used the reduced widths show a systematic 
deviation from the reduced widths obtained from the theory of Bhatia et al. 
(Hinds et al. private communication). Until the derivation of reduced widths 
from stripping reactions is made more precise it would be premature to 
speculate upon the discrepancies between the calculated reduced widths and 
the reduced widths derived from the experimental results. The use of Equa- 
tion 8 also has surprising success first of all in the case of Table I and later 
in Table V. In Table V the sum of the reduced widths equals 1.74 instead 
of the theoretically expected value of 1. This is reasonable agreement con- 
sidering that the experimental value depends quite sensitively on the assumed 
radius of the interaction. 

The agreement between theory and experiment in Tables I, IV, and V 
indicates that the further study of stripping reactions leading to rotational 
states both in light and heavier nuclei would be profitable. 


4. Gamma-Ray Transition Widths 

The usefulness of the concept of A-forbiddenness in describing the y-ray 
transitions in Mg*® and Al** has been discussed in Section 5, and the results 
for Mg*® and Al*® are given in Table III. 

Before discussing the results of the calculations of the 6 and G? coefficients, 
it is worth estimating to what extent such calculations can be relied upon. 
Using the relations given by Kerman (1956), it is possible to estimate the 
mixing of the various rotational states. The parameters which describe the 
mixing of the Al® states in the first K = 1/2 band (based on the first excited 
state at 0.45 Mev.) with those in the second K = 1/2 band (based on the 
2.51 Mev. state) are given in Table VI. The notation employed in Table VI 
is as follows: the perturbed wave function for the lower band can be ex- 
pressed as 


by = aps”*+bpz", 


TABLE VI 


THE MIXING AMPLITUDE, OF THE BAND (K = 

1/2) BASED ON THE 2.51 MEV. STATE IN Al?5 

IN THE WAVE FUNCTION OF THE BAND (K = 

1/2) BASED ON THE 0.45 MEV. STATE, IS 

SHOWN AS A FUNCTION OF THE J VALUE OF 
THE STATE 


a 
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where a?+6? = 1 and y,” is the unperturbed wave function for the lower 
band and y¥," is the unperturbed wave function for the upper band. J is the 
angular-momentum quantum number of the state. 

The values of a and 8, together with the observed degree of K-forbiddenness, 
indicate that the comparison of theory with experiment in the case of y-ray 
transitions is likely to be much less satisfactory than in the case of reduced 
particle widths unless an extensive calculation is made which takes into 
account the RPC. Such a calculation has been made for W'!** by Kerman 
(1956). 

The calculations of 6 and G? coefficients were carried out for Al*® assuming 
that the RPC was negligible. The formulae for dipole transitions used in 
the calculations have been given by Nilsson (1955). The results of such 
calculations, though encouraging, show that there is much room for im- 
provement. 

In Figs. 8 and 9 the branching ratios of a number of .W/1 and £1 transitions 
have been compared with the theoretical expressions of Alaga et al. (1955) 
to obtain empirical values of by; and bg:. A good fit to the branching ratios 
was obtained by a single value of the parameter } in each case. The Nilsson 
eigenfunctions permit these to be calculated a priori; the results of such 
calculations are shown in Fig. 14. All the £1 transitions considered involve 
transitions between rotational levels based on Nilsson’s orbit number 14 to 
ones based on orbit number 9, and for positive values of the parameter », 


~<-, 14-9 EI TRANSITIONS 1.2 
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Fic. 14. The calculated values of Gg, beiGei, bg: and Gy, byiGai, by are shown as a function 
of ». 


bg: ~ +0.6. This is to be compared with the empirical value of +0.3. Simi- 
larly, for the M1 transitions involving Nilsson’s number 11 to number 9 
orbits, the empirical value of b4;; = —0.5 is to be compared to the theoretical 
value of —1.0 obtained from Fig. 14. 
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An empirical value can also be obtained for Ggi? and Gy? from the 
measured absolute transition probabilities. The values are Gz," (14 — 9) 
= (0.13 and Gy,? (11 > 9) = 2.5. They can be compared to the theoretical 
values given in Fig. 14, about 0.023 and 11.2 respectively, which is good 
agreement for absolute gamma transition probabilities. The 6 and G? co- 
efficients are derived from all the measured transitions between different 
rotational bands which are not K-forbidden. 

Transitions within rotational bands are of considerable interest since it is 
here that enhanced £2 transitions would be expected. It has been shown by 
McManus and Sharp that, in addition to the single-particle contribution to 
the E2 transition probability given by Nilsson, an additional and, for values 
of the parameter 7 > 0.05, considerably greater contribution of a collective 
nature occurs for transitions within a rotational band. This collective E2 
transition is essentially negligible for transitions between different bands. The 
complete expression for £2 transition widths provided by McManus and 
Sharp has been quoted by Bromley et al. (1957) and has been employed in 
the discussion to follow. 

The results for transitions within the rotational band based on Nilsson’s 
level 9 are shown in Fig. 15. The theoretical branching ratio ¥4s/(T iat Tes) 
for the 2.73 and 1.81 Mev. states is plotted as a function of the distortion 
parameter 7. In both cases the experimental ratio is greater than unity. For 
the 1.81 Mev. state, the calculated M1 transition probability vanishes at 
values of » near +4 causing the ratio to be very large at these values. Agree- 
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Fic. 15. The experimental y-ray branching ratios within the band based on the first excited 
state of Al?5 are compared with theory. 

Fic. 16. The experimental E2/1/1 ratio for the transition from the state at 1.81 Mev. to 
the state at 0.95 Mev. in Al® is compared with theory. 
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ment with experiment occurs for n = +3.4 and +4.4. The agreement is not 
as good for the branching of the 2.73 Mev. state since the observed value 
is about 10 times greater than predicted. 

The measured £2—M1 mixture in the transition between the 1.81 and 
0.95 Mev. states is compared with theory in Fig. 16 and again good agreement 
is obtained both for the intensity ratio and phase with a distortion parameter 
n ~ 3.5. 

Similar comparisons for transitions within the rotational band based on 
Nilsson’s orbit 5 are shown in Figs. 17 and 18. Here, however, the measured 
branching for the 3.44 Mev. state is considerably larger than the predicted 
value while the measured phase for the E2-M1 mixture in the 3.44 — 1.61 
Mev. transition is opposite to the predicted value. 
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Fic. 17. The experimental y-ray branching ratio of the 3.44 Mev. state in Al** is compared 
with theory. 

Fic. 18. The experimental E2/M1 ratio for the transition from the state at 3.44 Mev. to 
the state at 1.61 Mev. in Al®* is compared with theory. 


The disagreement between some of the experimental results and the calcula- 
tions based upon Nilsson’s eigenfunctions possibly reflects the inaccuracy of 
the assumption that RPC is negligible (cf. Table VI). It is also possible that 
better agreement could be obtained by a limited variation of the parameters 
assumed by Nilsson for the nuclear potential. Such detailed calculations would 
be of value, particularly if further members of the rotational bands are found. 


IV. CONCLUSIONS 
It has been demonstrated that the strong-coupling wave function is a good 
approximation for the nuclei Mg®® and Al**. The evidence in favor of the 
strong-coupling wave function is as follows: 
(1) The appearance of well-defined rotational bands each of which has a 
sequence of states with the energy separation characteristic of rotational 
states. 
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(2) The nucleon widths for the excited states are in good agreement with 
the sum rule. 

(3) The y-ray decay of the excited states exhibits the phenomenon of 
K-forbiddenness. 

(4) A good fit to the ratios of the reduced dipole widths can be obtained 
using only two parameters 6g; and by71, one for the electric dipole transitions 
and one for the magnetic dipole transitions. 

The single-nucleon part of the strong-coupling wave function has been 
calculated by Nilsson (1955) using a specific spheroidal potential. Using these 
wave functions the decoupling parameters a, the gamma-ray branching ratios, 
and the nucleon reduced widths have been evaluated. Comparison of these 
with the results of experiment led to the following conclusions: 

(1) Although the bands observed for Mg*® and Al*® are predicted by the 
model, the order of the two lowest bands is inverted and one of the expected 
bands is missing. 

(2) The calculated decoupling parameters, for n ~ 4, for the three K = 1/2 
bands are in good agreement with the experimental values. 

(3) The equilibrium distortion of Mg*> and Al* is predicted to lie between 
7 = 4 and 7 = 6. 

(4) No good over-all fit to the reduced nucleon widths is possible. However 
the best agreement is obtained for 7 ~ 3. 

(5) The calculated values for the parameters by; and by; are in fair agree- 
ment with experiment. 

(6) The branching ratios within the lowest bands are only in approximate 
agreement with experiment. 

The evidence presented here would seem to indicate that the rotational 
model has real significance for this region of atomic number. Although the 
Nilsson model of a single nucleon in a deformed potential fails to give quanti- 
tative agreement with all the experimental data, this failure is not surprising 
in view of the relative simplicity of the model. Preliminary calculations by 
Elliott (1957) indicate that there may be a rather close connection between 
the shell model and the collective model and, if this is so, it provides some 
theoretical justification for applying the collective model to nuclei in regions 
between closed shells where calculations using the intermediate coupling shell 
model are prohibitively tedious. 


The authors wish to record their indebtedness to Professor A. Bohr, Dr. 
B. R. Mottelson, Dr. S. G. Nilsson, Dr. A. kX. Kerman, and Dr. K. Gottfried 


for several stimulating discussions and communications. 


APPENDIX A. PERTURBED BANDS IN Na® 
Reference to Fig. 10 indicates that Na?’ is expected to have K = 3/2 and 
K = 5/2 for its first two bands. Such bands, if separated by the order of 
1 Mev., would be expected to interact strongly, and in fact the perturbed 
order of the states of Na**, from the point of view of the rotational model, 
can be explained in this way. Fig. 19 shows a possible fit to the lower states 
of Na** together with the assumed parameters. The nucleus Ne** would, 
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however, be expected to show the same sort of band structure as Mg*® and 
Al** amongst its low-lying states. This is because the odd nucleon in Ne? is 
the fifth neutron after the closed shell at O'*®, which is the case for Mg*® 
(fifth neutron) and Al*> (fifth proton). There is at the moment no evidence 


Mev 
F ~ LOW STATES OF Na 
ise 
3 J 7 E(Mev) 
(%+) 2.398 
. 55455 (+) 2.076 
IL %+% 
%4+ 0.440 
3% +34 ——_—_— ei 


(a) (b) 


Fic. 19. <A possible description of the low-lying states in Na® based upon the two rotational 
bands K = 3/2 and K = 5/2 with the rotation-particle-coupling RPC parameter Ag = 440 
kev. and the parametersh?/2/ = 283 kev. and 340 kev. respectively. (a) shows the unperturbed 
overlapping bands (Ax = 0) and (+) shows the perturbed band (Ag = 440 kev.). 


for rotational-like bands in Ne**, but the measurement of the angular distri- 
butions from the stripping reaction Ne??(d, p)Ne?* by Burrows et al. (1956) 
to the ground and first excited states shows prominent / = 2 and / = 0 pat- 
terns respectively, which is similar to the Mg*‘(d, p)Mg*® reaction studied 
by Holt and Marsham (1953). This is to be contrasted with Ne?°(d, p)Ne?!, 
where the angular distributions show isotropy for the ground-state transition 
and an / = 2 pattern for the first-excited-state transition (Burrows et al. 


1956). 
APPENDIX B. STATES OF CORE EXCITATION IN AI? 

As mentioned in Section 5, the strong inhibition of the ground-state tran- 
sitions from the states in Al®® at 4.22 and 4.60 Mev. requires some discussion, 
which can be based upon Figs. 10 and 11. The first states in Al*> of core 
excitation, contrasted with states arising from the excitation of the odd pro- 
ton, should involve the excitation of a nucleon from orbit number 7 to orbit 
number 5 or 9 which are close together. In either case the unpaired proton 
in orbit number 7 might generate a K = 3/2 band. The possibility gives a 
plausible explanation for the small reduced proton widths for the states at 
4.22 and 4.60 Mev. Also if the two excited nucleons are considered to be in 
the orbit number 9 dipole, y-ray transitions to the ground state are two- 
particle transitions whereas transitions to the band based on the first excited 
state are single-particle transitions which are K-allowed. The potential- 
energy curves of Fig. 11 support this assignment. If the two nucleons are in 
orbit number 5 the ground-state transitions would be expected to be strongest. 
It may not be coincidental that the separation 380 kev. of the 3/2+ and 
5/2+ states is close to the values for the separations of the ground and first 
excited states of Ne*! and Na®* (347 kev. and 440 kev. respectively) which 
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can be interpreted as a perturbed 3/2+, 5/2+ pair (cf. Fig. 19). The first 
state in Mg*‘ which is not a member of the ground-state band is the 2 

state at 4.23 Mev. The prominent / = 0 stripping pattern for this state 
observed in the Na**(d, n)Mg** reaction strongly suggests that the parents of 
the 4.22 Mev. and 4.60 Mev. states in Al*?> are members of a band based upon 
the 4.23 Mev. state in Mg?*. Evidence for such a band has been obtained by 


Newton (1954); the Na?3(d, ~) Mg*‘ reaction is discussed briefly in Appendix D. 


APPENDIX C. PROTON WIDTHS FOR THE UNBOUND 
ROTATIONAL STATES IN AI? 

Equation 6 enables the calculation of reduced widths for the emission of 
protons from rotational states in Al*> to the first excited state of Mg**. Since 
no inelastic scattering from the rotational states belonging to the first four 
observed bands in Al*® has been observed, the calculations can only serve 
as a guide to future experiments to locate the missing members of the bands. 

Consider, for example, the missing 9/2+A = 1/2 state belonging to the 
band based on the first excited state of Al**. If this state occurs near the 
expected bombarding energy of 1.7 Mev. then it can be shown, using the tables 
for barrier penetrability given by Gove (1957), that the inelastic d-wave proton 
width is approximately 25 millivolts. The inelastic proton width increases 
rapidly with energy because of the proximity of the inelastic threshold and 
at 2.0 Mev. bombarding energy becomes 13 ev. As in each case the g-wave 
elastic proton width is expected to be smaller than the inelastic proton width 
the resonance would be difficult to detect. This conclusion probably applies 
to the higher spin members of the rotational bands in Al**. In addition to 
the barrier effect the reduced formation width in each case is expected to be 
very low because the distance of the lg and 1h from the ld and If shells 
results in small admixing. The protons emitted to the first excited state of 
Mg** can however be p, d, or f protons with large reduced widths. 

Equation 6 also shows another way of testing the goodness of the K-quantum 
number since the Clebsch-Gordan coefficient vanishes unless K;+2 = Ky. A 
state whose (J7K) values are (3/2 + 3/2) for example cannot break up into 
an s-wave proton plus Mg” in its first excited state since the above equality 
is violated. If however there was an admixture of K = 1/2 in the state, s-wave 
protons could be emitted. 


APPENDIX D. DEUTERON STRIPPING AND THE 
ROTATIONAL MODEL NEAR A = 25 

Recent studies of the (d,”) and (d, p) reactions by Calvert et al. (1955) 
and Holt and Marsham (1953a, 6) have yielded information concerning 
nucleon widths in the vicinity of A = 25. The Na®*(d, n)Mg* reaction is of 
particular interest because of the rotational spectrum of Mg**. The stripping 
patterns allowed by conservation of angular momentum for the 0+, 2+, 
and 4+ states are / = 2,/ = 0 and/ = 2, and / = 2 respectively. However, 
if the rotational model gives a good description of Mg** Equation 6 and Fig. 13 
can be used to determine the reduced widths as a function of distortion. 
Unfortunately the low yield of neutrons to the low-lying states of Mg?! 
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prevented a determination of the reduced widths (Calvert et a/. 1955). How- 
ever it is possible to make some qualitative statements concerning the relation 
of the results to the rotational model. No marked stripping pattern for the 
first excited state of Mg** was observed in spite of the fact that an / = 0 
transition is allowed from the conservation of angular momentum. Since 
/ = 0 patterns even with small proton capture probabilities can be observed 
by the (d, ) reaction,* this implies that the addition of an s-wave proton 
to Na** to form the first excited state of Mg** is strongly inhibited by some 
selection rule. The rotational model provides such an inhibition because since 
K, = 3/2 and K, = 0, the addition of s-wave protons, 2 = 1/2, is forbidden. 

Fig. 13(6) also shows that the proton capture probability to the ground 
state of My*! is very small for prolate distortions whereas the proton capture 
probability for the first excited state is near the maximum value possible. For 
a distortion 7 = +4 the ratio of the proton capture probabilities to the two 
states is approximately 20. The differential cross section at 0° for the ground- 
state neutron group is a factor of five smaller than that for the first excited 
state, which is possibly a lower limit to the ratio of the proton capture prob- 
abilities in view of the uncertainties introduced by compound-nucleus and 
other effects (cf. footnote to Section 4). The remeasurement with improved 
statistics of the proton-capture probabilities for the two lowest states of 
Meg?! from the Na?3(d, 2) Mg** reactions would be of considerable interest. 

The angular distribution of the sum of the neutron groups to the second 
and third excited states (unresolved) of Mg?! shows a prominent / = 0 
stripping pattern together with a possible / = 2 pattern. The / = 0 pattern 
must be associated with the 2+ state at 4.23 Mev. and the large proton 
capture probability implies that a rather simple description of the state is 
possible.t A description based upon Fig. 10 would be that the unpaired 
proton of Na** in orbit number 7 pairs with a proton in orbit number 9 to 
produce a state with J = 2+ and with A = 2. This state is then the base 
state for a rotational band. Evidence for the second member of this band 
J = 3+ has been given by Newton (1954). A possible connection of this 
band with Al*®> has been considered in Appendix B. 

The Mg**(d, p) Mg?® studied by Holt and Marsham also is of interest since 
the target nucleus is well described by the rotational model. Comparatively 
little is known about the properties of Mg?®. The first excited state is most 
probably 2+ and the excited states do not suggest a rotational structure as 
in the case of Mg?‘. The angular distribution of the protons leading to the 
first excited state of Mg*® shows a mixed / = O and / = 2 pattern. The reduced 
neutron capture probability is 20 times smaller for the addition of an / = 0 
neutron to Mg* to form Mg”*® in its first excited state than for the addition 
of an / = 2 neutron. This observation is quite consistent with the rotational 
model since the first excited state of Mg?* would be expected to be 2+ with 
K = 0. The inhibited / = 0 patterns for the Mg*4 and Mg’® excited states 

*The prominent / = 0 deuteron stripping pattern for neutrons leaving Si? in its first excited 


state has a proton capture probability a factor of 10 times smaller than the observed ground 


state / = 2 pattern (Calvert ef al. 1955). 
tThe / = 2 component could be from either or both of the 2+ or the 4+ member. of the 


doublet. 
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could also be explained on the basis of nearly pure j-j coupling. Both nuclei 
are then in the ds5,2 shell which would imply no / = 0 transitions to states 
which are described by a number of d5,2 particles. On the rotational picture 
the appearance of an / = 0 component to the angular distribution implies 
a small admixture of other K-values to the predominantly K = 0 first excited 
state. It is not possible to correlate further the appearance of forbidden 
stripping patterns with the deviation from a rotational spectrum until more 
experimental evidence is available. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


The Influence of Paramagnetic Resonance Saturation on the Faraday Effect 


The Faraday effect (the rotation of the plane of polarization of light passing through a 
material in a magnetic field) is very large, in general, in paramagnetic salts. The rotation 
of the plane of polarization is proportional to the magnetic moment of the salt, and for a 
single crystal of neodymium-ethylsulphate the maximum rotation corresponding to complete 
magnetization is about 114° per mm. of path length for green mercury light (A = 5461 A) 
propagated along the optic axis (Becquerel et a/. 1938). When the crystal is exposed to micro- 
wave radiation of the frequency which causes paramagnetic resonance, and the microwave 
power is sufficient to cause appreciable saturation of the resonance, the magnetic moment is 
reduced and the Faraday rotation should also be reduced. This effect was predicted by Kastler 
(1951) and a theory of this effect has been developed by Opechowski (1953). 

We have observed this effect in a single crystal of neodymium-ethylsulphate. The crystal 
was mounted with the optic axis parallel to the external d-c. magnetic field, in a TE;o rect- 
angular cavity, with two holes, about 2 mm. in diameter, in the large faces to let the light 
through. Microwave power at 9060 Mc./sec. from a 2K39 klystron was fed into the cavity, 
which was critically coupled to the microwave transmission line and had an unloaded Q of 
about 9000. The cavity was immersed in liquid helium at about 1.5° K. The light, from a 
mercury lamp, was filtered to leave only the green line (A = 5461 A). It was then passed 
through a nicol prism polarizer, through the hollow poles of the electromagnet, through the 
cavity, out through the magnet, and finally through a nicol prism analyzer. The emergent 
light was condensed onto the cathode of a photomultiplier (1P21), and the anode current 
was taken as a measure of the light passing through the system. The light passed along the 
optic axis of the neodymium-ethylsulphate crystal, and the path length was 1.6 mm. The 
ends of the crystal were cut perpendicular to the optic axis and polished. The analyzer was 
rotated until the photomultiplier current was a minimum, and the Faraday rotation was 
read off from a scale on the analyzer housing. 

The results are shown in Fig. 1. The angle of rotation of the plane of polarization is plotted 
as a function of magnetic field for four cases: no microwave power, full microwave power, 
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medium power, and a low value of power. The power levels have not been determined. These 
investigations are being continued and will be reported at a later date. 
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